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Voorwoord 


Door het toenemende gebruik van microcomputers maken de grafische toepas- 
singen een sterke ontwikkeling door. Het beeldscherm is tegenwoordig, veel 
meer dan de printer, het afdrukapparaat bij uitstek. Grafische toepassingen 
(graphics) vinden we bij computerspelletjes, bij computerkunst en in toene- 
mende mate op de werkplek van ingenieurs, ontwerpers, architecten en 
anderen die zich bezighouden met Computer Ondersteund Ontwerpen (Com- 
puter Aided Design). 


De geïnteresseerde leek is vaak verrukt van de mooie 'plaatjes', die bij 
demonstraties getoond worden, en denkt dat hiervoor hele ingewikkelde pro- 
gramma's nodig zijn. Iemand, die iets van wiskunde afweet, weet hoe een 
functie eruitziet en kan doorgaans zonder veel moeite grafische programma's 
maken. 


Ik heb voor het Zwitserse Computertijdschrift 'MIKRO+KLEINCOMPUTER' een 
cursus ‘programmeren met hoge resolutie! geschreven, die in een aantal afle- 
veringen van het blad is verschenen. Al direct na het verschijnen van het 
eerste artikel werd ik werkelijk overstelpt met enthousiaste reacties van de 
lezers. Dit heeft ons aangezet tot het maken van dit boek. 


Veel van de grafische toepassingsprogramma's worden geschreven voor een 
bepaald merk en type microcomputer of voor een bepaalde plotter. Wie een 
dergelijk programma wil herschrijven voor een ander merk microcomputer zal 
veel moeilijkheden ondervinden en vaak zelfs zijn pogingen staken. 


De veertig programma's in dit boek zijn geschreven in standaard BASIC en 
gebruiken geen POKE- en PEEK-opdrachten. Bovendien worden slechts twee 
grafische opdrachten gebruikt, te weten het inschakelen van de hoge-resolu- 
tiestand en het verbinden van twee punten door een rechte lijn. De illustra- 
ties die als voorbeeld dienen van hetgeen de programma's tekenen zijn 
gemaakt op de miniplotter van de Sharp PC-1500. 


Wij hopen dat de lezers veel plezier aan deze programma's zullen beleven en 
dat zij voor velen een aanmoediging zullen zijn om het terrein van de compu- 
tergraphics verder te verkennen. Mijn dank gaat uit naar de heer H.J. 
Ottenbacher van Micro+Kleincomputer; zonder hem zou dit boek nooit het 
daglicht hebben gezien. 


Marcel Sutter 
voorjaar 1984 


Voorwoord bij de Nederlandse 
uitgave 


De auteur, Marcel Sutter, heeft dit boek oorspronkelijk geschreven 
in standaard Mierosoft-BASIC. In de programma's maakte hij alleen 
gebruik van de opdracht om twee punten met elkaar te verbinden. 
In de MSX-versie die wij van zijn boek hebben gemaakt geven we 
daarentegen suggesties voor het werken met kleuren en laten we 
bij een groot aantal programma's zien hoe u met de CIRCLE- en 
PAINT-opdrachten de programma's kunt uitbreiden en verfraaien. 
In de appendix vindt u een aantal uitgewerkte voorbeelden. Hierin 
laten we ook zien hoe u tekst en getallen op het grafische scherm 
kunt afdrukken. 


Alle programma's zijn ingetoetst op een MSX-computer en direct van- 
uit het geheugen van de computer op een printer afgedrukt. Als u 
de programma's precies zo intoetst als ze in het boek staan afge- 
drukt, moeten ze allemaal goed werken. De programma's zijn inge- 
toetst met het toetsenbord in de 'kleine-letterstand'. Hierdoor wordt 
de tekst in de commentaarregels en in INPUT- en PRINT-opdrachten 
in kleine letters afgedrukt, terwijl de computer van alle BASIC- 
opdrachten (zoals rem, print, for, input, line, enz.) en van alle 
variabelen (zoals xl, u, h, k, yy, enz.), die we met kleine letters 
intoetsen, zelf hoofdletters maakt. U kunt, als u dat wilt natuurlijk, 
ook alles in de 'hoofdletterstand' intoetsen. 


Wij raden u aan reeds bij de eerste programma's te experimenteren 
met de COLOR-opdracht (verschillende kleuren voor de voorgrond, 
de achtergrond en de boven- en onderrand), met de LINE-opdracht 
(verschillende afdrukkleuren), met de PAINT-opdracht (voor het 
kleuren van vlakken) en met de DRAW-opdracht (voor het afdruk- 
ken van tekst op het grafische scherm). U kunt uw ervaringen dan 
direct bij volgende programma's gebruiken. In de tekst vindt u 
suggesties hoe u de bovengenoemde opdrachten kunt gebruiken. 


De auteur heeft voor de tekeningen gebruik gemaakt van een Sharp- 
plotter met een resolutie van 220 bij 220 puntjes. Op een aantal 
plaatsen hebben wij zijn tekeningen vervangen door afdrukken, die 


wij zelf van het grafische scherm van de computer gemaakt hebben. 
Deze afdrukken geven een beter beeld van hoe u de tekeningen en 
grafieken op uw eigen MSX-scherm ziet. 


Alle programma's zijn kort. Verfraai ze, verbeter ze en breid ze 
uit. De mogelijkheden zijn onbegrensd. 


Het is een boek voor elke MSX-computer. Diegenen die de veelal 
wiskundige uitleg bij de programma's teveel van het goede vinden, 
zullen alleen al door het intikken en draaien van de programma's 
verrukt zijn van het resultaat. Anderen, die zich nog iets van de 
sinus en cosinus uit hun schooltijd herinneren of die nu op school 
zitten, zullen weinig of geen moeite hebben met de theorie in dit 
boek. Leraren en leerlingen kunnen wellicht hun voordeel doen met 
de programma's voor het tekenen van functies. Ook het driedimen- 
sionaal tekenen wordt behandeld. Erg leuk zijn de vijf programma's 
waarmee in BASIC getekend kan worden zoals dat met LOGO ook 
kan. Dit geeft zeer fraaie grafieken, vooral de turtle-graphics. De 
laatste vijf programma's zijn toepassingsprogramma's, waarin 
‘graphies! de hoofdrol spelen. 


Hopelijk hebt u net zoveel plezier met het draaien en verfraaien van 
deze programma's als wij hebben gehad bij de bewerking van dit 
boek. 


Nok van Veen 
februari 1985 
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In de Hoge-Resolutiestand (screen 2) bestaat het beeldscherm van 
een MSX-computer uit 256 beeldpunten horizontaal en 192 beeldpun- 
ten verticaal. Vanuit een BASIC-programma kunnen we elk van deze 
49152 puntjes aan- en uitzetten (PSET en PRESET); we kunnen 
lijnen trekken tussen twee beeldpunten (LINE); we kunnen cirkels 
tekenen (CIRCLE); we kunnen vlakken kleuren (PAINT) en we 
kunnen van een speciale grafische macrotaal gebruikmaken (DRAW). 
De hoge-resolutiestand kiezen we met de SCREEN-opdracht, terwijl 
we de afdrukkleur, de achtergrondkleur en de kaderkleur kunnen 
instellen met de COLOR-opdracht. 


In de programma's gebruiken we in principe de volgende drie gra- 
fische opdrachten: 


COLOR 1,15,4 : afdrukkleur is kleur 1 (zwart); 
achtergrondkleur is kleur 15 (wit); 
boven- en onderranden (kader) 
zijn kleur 4 (blauw) 


SCREEN 2 : stelt het hoge-resolutiescherm in 


LINE (X1,Y1)-(X2,Y2) : verbindt twee beeldpunten met 
beeldschermcoördinaten (x1,y1) en 
(X2,y2) door een rechte lijn 


Ook erg mooi is een roodbruine achtergrond (kleur 6) en een gele 
afdrukkleur (kleur 10 of 11). Dit stellen we in met COLOR 11,6,6. 
Hierbij is er geen boven- en onderrand maar een achtergrond in 
één kleur. 


Met de SCREEN 2-opdracht krijgen we de beeldschermindeling van 
p-.2 bovenaan. Linksboven in het beeldscherm ligt de oorsprong van 
het beeldschermcoördinatenstelsel. De bovenrand van het scherm 
noemen we de X-as. De linkerrand noemen we de Y-as. De Y-as 
wijst dus naar beneden. Een punt met een grotere Y-coördinaat ligt 
dus onder (voor ons) een punt met een kleinere Y-coördinaat. Als 
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(0,191) (255,191) 


middelpunt van het scherm nemen we het punt (U,V) met U=128 en 
V=96. Wellicht moet u even wennen aan de, naar beneden wijzende, 
Y-as. 


(256,0) 


Ya 


P‚(xa sy) 


(0,192) 


De LINE-opdracht trekt een lijn tussen twee beeldschermpunten in 
de afdrukkleur die met COLOR is ingesteld. We moeten echter 
oppassen. Het volgende programmaatje zal op het scherm een recht- 
hoek tekenen en géén vierkant! 


10: COLOR 1,154 £-- SCREEN. 2 
20 LINE (50,50)-(150,50) 
30 LINE (150,50)-(150,150) 
40 LINE (150,150)-(50,150) 
50 LINE (50,150)-(50,50) 
60 AS=INPUT$(1) 

70 END 


MSX-graphics 3 


De horizontale zijden zijn 1,4 maal zo lang als de verticale zijden. 
Als we een vierkant willen tekenen, zullen we alle acht de X-coördi- 
naten in de regels 20-50 door 1,4 moeten delen. Omdat we in vrijwel 
alle programma's een verhouding van 1 op 1 voor de x- en y-rich- 
ting willen zien, passen we de X-coördinaten in de LINE-opdrachten 
(en ook in de CIRCLE- en PAINT-opdrachten) als volgt aan: 


oude X 


nieuwe X = 37 + 


Hierdoor komt de rechterbovenhoek (X=255) terecht op het beeld- 
punt met X = 37 + 255/1,4, dus X=219,143 en de oorsprong (X=0) 
komt op X=37. Omdat beeldschermecoördinaten altijd gehele getallen 
zijn, nemen we de INT-functie bij het omrekenen van de X-coördi- 
naten. Regel 20 van het vorige programma zou er dan zo uitzien: 


20 LINECINT(37+X1/1.4) „Y1)-CINT(357+X2/1.4),Y2) 


In de programma's doen we het iets anders. We maken van deze 
omrekening een functie met de DEF FN-opdracht, die er zo uitziet: 


DEF FNXOO=INT (37+X/1.4+.5) 


De X tussen haakjes is het argument van de functie. De functie heet 
FNX. Dit argument, dat in een programma elke variabele mag zijn, 
wordt door 1,4 gedeeld; er wordt 37 bij opgeteld en door er nog } 
bij op te tellen voordat de INT-functie genomen wordt, ronden we 
de nieuwe X-coördinaat af. In het vorige programma zou dit er zo 
uitzien : 


10 COLOR 1,15,4 : SCREEN 2 

15 DEF FNXOX=INT (374X/1.4+.5) 

20 LINE (FNX(O),0)- (FNX (100) ,0) 
30 LINE (FNX (100) „0)- (FNX (100) „100) 


“ enzovoorts 


In vrijwel elk programma komt u deze functie tegen. De algemene 
LINE-opdracht voor het trekken van een lijn tussen de punten 
Pi(x1,y1) en P2(x2,yz2) ziet er dus zo uit: 


LINE(FNX(X1),Y1)-(FNX(X2),Y2) 
Deze LINE-opdracht trekt een lijn in de kleur die met COLOR als 
afdrukkleur gekozen is. Kijk in uw MSX-handleiding wat u allemaal 
nog meer met LINE kunt doen. 


Bij hoge-resolutie-graphics onderscheiden we puntgraphics en lijn- 
(of vector-) graphics. 
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Puntgraphics 


Bij puntgraphies worden de coördinaten x en y van een bepaald 
punt berekend en wordt dit punt door een bepaalde graphic- 
opdracht (PSET(x,y)) op het scherm zichtbaar. Willen we een der- 
gelijk punt weer onzichtbaar maken, dan is daar een andere 
opdracht (PRESET (x,y)) voor nodig. In dit boek zullen we echter 
voor het tekenen van figuren niet van deze puntgraphics-techniek 
gebruik maken. 


Lijngraphics 


Bij lijngraphies worden door een bepaalde opdracht, die als 
machineroutine in de microcomputer aanwezig is, razendsnel twee 
berekende punten Pi(X:,Y1) en P2(X2,Y2) door een rechte lijn met 
elkaar verbonden. Zo'n rechte verbindingslijn wordt opgebouwd uit 
een groot aantal kleine horizontale en verticale lijnstukjes. Op een 
afstand van zo'n 50 em van het scherm lijken deze verbindingslijnen 
inderdaad bijna recht. 


In onze programma's maken we gebruik van deze lijngraphics. Alles 
wat we gaan tekenen ontstaat door het steeds verbinden van twee 
punten door een rechte lijn (of door een recht lijntje). De program- 
ma's zijn kort, zodat duidelijk wordt op welke wijze de tekeningen 
ontstaan. U kunt de programma's zelf uitbreiden en verfraaien. U 
kunt kleuren gebruiken of regels toevoegen waarmee de door u in 
te toetsen waarden gecontroleerd worden, zodat dit niet leidt tot 
het afbreken van het programma. In de appendix vindt u een aantal 
uitbreidingen van enkele programma's. Hierin ziet u hoe u met kleu- 
ren kunt werken, hoe u vlakken vult en hoe u tekst op het gra- 
fische scherm kunt afdrukken. 


Programmastructuur 

In alle programma's gebruiken we steeds dezelfde variabelen en 
dezelfde programmastructuur. De verschillende variabelen beteke- 
nen : 


X1,Y1l coördinaten van het beginpunt 
X2,Y2 coördinaten van het eindpunt 


U,‚V oorsprong van het wiskundige coördinatensysteem; 
dikwijls is dit het midden van het beeldscherm (128,96) 

H de waarde 0,5 voor het afronden op helen 

K de vermenigvuldigingsfactor voor functiewaarden 

W‚W1 hoeken bij trigoniometrische functies 


RD het getal 1/180 voor het omrekenen van graden in radialen 
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Een grof structuurdiagram voor de programma's ziet er zo uit: 


doe zolang nodig 


nieuw punt P2(X2,Y2) berekenen 
verbindt Pl met P2 


punt P2 wordt punt Pl (X1=X2 


Y1=Y2) 


In alle programma's gebruiken we twee belangrijke transformatie- 
formules. Voor het schermcoördinatensysteem van de MSX-computer 
is de linkerbovenhoek van het scherm de oorsprong. Wijzelf gebrui- 
ken echter bijna altijd een coördinatensysteem met de oorsprong in 
het midden van het scherm. Om de beeldschermcoördinaten X,Y 
(zoals de computer ze gebruikt) te berekenen uit de coördinaten x,y 
zoals wij ze gebruiken hanteren we de transformatieformules: 


X=INT(U+*x+H) : Y=INT(V=y*+*H) 
Bekijk de hiernaast- ï 
staande figuur eens. \ 
Duidelijk is dat elk (0,0) X X 
punt p(x,‚y) in ons 


coördinatenstelsel 
(gestippeld assen- 
kruis) door boven- 
staande formules 
wordt getransfor- 
meerd (overge- 
bracht) naar een 
punt P(X,Y) in het 
coördinatenstelsel 
van de computer. 
Hierbij hoeft onze 
oorsprong (U,‚V) 
natuurlijk niet 
altijd precies in het midden van het beeld te liggen. 


Y 


Genoeg theorie, nu de programma's. De ideeën voor de programma's 
hebben wij gekregen bij het doorbladeren van Amerikaanse, Duitse en 
Franse computertijdschriften. Alle programma's zijn echter eigen 
creaties of bewerkingen en vereenvoudigingen van reeds gepubli 
ceerde programma's. 
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Programma 1 tekent een 'diagonaalweb'. Elk van de acht bovenste 
punten wordt door een rechte lijn verbonden met elk van de onder- 
ste acht punten. 


100 REM programma 1 diagonaalweb 
110 CSS 

120 COLOR 1,15,4 : SCREEN 2 

130 Y1=0 : Y2=192 


NEET INZET DTT 
NANDA 
Nakheebedbedpe 00e sgf 

„ 


AN 
He 


Arsen ens df 
…. ves ed 


/ 


140 FOR X1= 16 TO 240 STEP 32 gdsdnesn 

150 FOR X2=16 TO 240 STEP 32 bi RA AAAAA 
160 CENE (XY 1D CHEN 2I BVs eie \ A 
170 NEXT X2 | 9, 


180 NEXT X1 
190 AS=INPUTS$(1) 
200 COLOR 15,4,4 : END 


verse 

ere ereaih 
AAS 
BENEN 


Als u na regel 120 opneemt: 125 DEF FNXOX)=INT(374+X/1.4+.5) 
en regel 160 verandert in: 160 LINE CFNXCX1),Y1)=(FNX(X2) ,„Y2) 
hebben de x- en de y-richting dezelfde orde van grootte. 

Het bovenstaande programma geeft een nogal brede figuur. 


Programma 2 iseen fraai voorbeeld van het moiree-effect. Geniet ervan. 


100 REM programma 2 moireeveffect 
110 CES 

120 COLOR 1,15,4 : SCREEN 2 

130 FOR J=0 TO 136 STEP 4 


140 LINE (60,0)-(J+60,192) 

150 LINE (60,0) (196, 19214) 
160 NEXT J 

170 FOR J=0 TO 136 STEP 4 

180 LINE (1967192160 192104) 
190 LINE (196,192)-(60+J,0) 

200 NEXT J 


210 AS=INPUT$S(1) 
eeU COLOR 157444 5 END 


EE 


mr 


| 
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Programma 3 tekent een reeks driehoeken in perspectief. Het 'cen- 
trum van vermenigvuldiging' heeft de coördinaten S(0,72). De drie 
hoekpunten van de driehoek die vermenigvuldigd wordt zijn 
P1(6,84), P2(20,81) en P3(12,66). Alle coördinaten zijn beeld- 
schermcoördinaten (oorsprong links bovenaan). De drie richtingen 
van vermenigvuldiging (richtingsvectoren) zijn SP,, SP, en SP3. 
Deze worden in het programma als SX(1);SY(1), SX(2);SY(2) en 
SX(3);SY(3) vastgelegd. De vermenigvuldigingsfactor K is de lus- 
variabele van de buitenste FOR-NEXT lus. K loopt van 0 tot 10 in 
stappen van 0.4. Hieronder zien we de uitgangssituatie. het resul- 
taat van het programma en het programma zelf. 


(0,0) 6 12 20 „ sl 
KN 


100 REM programma 3 driehoeken in perspectief 
TĲD GES 

120 COLOR: 1,15,4 : SCREEN 2 

130 DIM KCS ESKCH) SN CSI 

140 FOR J=1 TO 3 

150 READ SXCJ),SYCI) 

160 NEXT J 

170 X(0)=0 : Y(O)=72 

180 FOR K=0 TO 10 STEP .4 


190 FOR J=1 TO 3 

200 XIX CONHKESK CI) 
210 YEII=YCONHKESY (JI) 
220 NEXT J 


230 LEINE (ECN ANISECR ELD eN CZ 
240 LINE (NCR AN CIN CHI NCSI 
250 LINE (X(3),Y(3))I- CCD, YCHDD 
260 NEXT K 

270 AS=INPUTS$(1) 

280 COLOR 15,4,4 : END 

290 DATA. '6,12,20,9,12,"6 
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Programma 4 tekent een reeks ingeschreven regelmatige zeshoeken. 
Behalve bij de eerste zeshoek, liggen alle zes de hoeken van elke 
zeshoek precies in het midden van een zijde van de omgeschreven 
zeshoek. De coördinaten van de hoekpunten van de eerste (buiten- 
ste) zeshoek worden met trigoniometrische functies (sinus- en 
cosinusfunctie) berekend. De hoek die hierbij als argument van de 
functies gebruikt wordt doorloopt de waarden 0°, 60°, 120°, 180°, 
2409, 300° en 360°. De middelpunten van de zijden van een zeshoek 
worden aangegeven met MX(0);MY(0) t/m MX(5);MY(5). Dit worden 
de hoekpunten van de volgende zeshoek. 


100 REM programma & ingeschreven 
zeshoeken 

110 CES 

120 COLOR 1,15,4 7 SCREEN 2 

130 DEF FNX(X)=INT(37+X/1.4+.5) 

140 DIM X(6),YC6) ,MX(6) „MY(6) 

150 PI=4&kATN(I) 

160 U=128:V=96:R=100: HE. 5 

170 W=60*PI/180 

180 FOR J=0 TO 6 


190 Wi=JtW EN 
200 XUJDZINTCUFR*COS (WI) +H) Van 
210 VOID ZINT (V=RESINCWI) +H) (en 
220 NEXT J SEN 
230 FOR N=1 TO 20 es 
240 FOR J=0 TO 5 Sf 
250 LINE CENXCX(JDD,Y(IDDT 
CENXCKCJHID DY CJHID) 
260 NEXT J 
270 FOR K=0 TO 5 
280 MXCKDZINTCOKCKD+X K+) /2+H) 
290 MYCKD=INTCCYCK)+Y K+) /2+H) 
300 NEXT K 
310 MXC6)=MXCO) : MYC6)=MYCO) 
320 FOR J=0 TO 6 
330 XOJDEMXOJD 2 YI EMY ID 
340 NEXT J 
350 NEXT N 


360 AS=INPUT$S(1) 
370 COLOR 15,4,4 : END 


Zie de appendix voor een speciaal effect met dit programma. 
Het gebruik van de FNX-functie (regel 130) is hier essentieel, 
anders wordt het een in horizontale richting uitgerekte figuur. 
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Programma 5 tekent alle diagonalen in een n-hoek. De n hoekpun- 
ten liggen op de omtrek van een ellips of van een cirkel afhankelijk 
van de waarden die we in het begin van het programma voor de 
halve lange as a en halve kleine as b kiezen. Voor de berekening 
van de coördinaten van de punten op de omtrek van de ellips ge- 
bruiken we niet de (cartesische) vergelijking (x?/a?) + (y?/b?) = 1, 
maar de parametervorm x = a.cosp en y = b.sinw. 


Als de parameter w loopt van 0° tot 360°, dan beschrijft het daar- 
aan toegevoegde punt P(x,y), met x = acospen y = bsinw, precies 
de omtrek van een ellips. Hieronder zien we een bepaald punt 
P(x1,y1) op de omtrek van de ellips met de daarbij horende waarde 
van de parameter wp). Voor dat punt geldt: xj = acosp; en 

Yi = bsinp:. 


Als xj =acosp;, en yi =bsinp, dan geldt ook 
x? =a?cos? pi: en y} =b?sin?p, en ook 


xî yi 
+ =COos? Ei en br =sin?p,, dus dan is 


u 


bz eos? pi + sin? wi 
nn A 
a 7 b Ì 


en hebben we de cartesische relatie tussen x; en y: afgeleid uit de 
parametervorm. 


en omdat cos?; + sin?p: gelijk is aan 1 geldt 


Voor het berekenen van de coördinaten van het j-de hoekpunt 
(X(J),Y(J)) verdelen we de maximale middelpuntshoek van 360° in 
n gelijke hoeken van elk 360/n graden. 
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Als A de lengte van de halve lange as is èn 
B de lengte van de halve korte as èn 
N het aantal hoekpunten van de n-hoek èn 


W gelijk is aan 360/N èn 

Wi de hoek die hoort bij het j-de hoekpunt èn 

X(J) de x-coördinaat van het j-de hoekpunt t.o.v. het 
middelpunt van de ellips èn 

Y(J) de y-coördinaat van het j-de hoekpunt t.o.v. het 
middelpunt van de ellips 


dan zou je in een BASIC-programma de coördinaten X(J) en Y(J) 
als volgt kunnen berekenen: 


w = 360/N : WI = Jx*W 
X(J) = INT(A*COS(WI)) 
Y(J) = INT(B#SIN(WIJ)) 


Als we dit doen maken we twee fouten. De eerste fout is dat de 
hoek W1, als argument van de COSinus- en de SINusfunctie, in 
graden is uitgedrukt, terwijl BASIC vereist dat het argument van 
deze functies in radialen wordt uitgedrukt. De tweede fout is dat 
het HRG-systeem*) van de MSX-computer de oorsprong van het 
coördinatenstelsel voor het scherm in de linkerbovenhoek van het 
scherm verwacht en niet in het middelpunt van de ellips. Wat we 
dus nog moeten doen is: 


8: bereken Wl in radialen èn 
b. transformeer X(J) en Y(J) naar beeldschermcoördinaten. 


We weten dat een hoek van 360° overeenkomt met 2n radialen, 
waarin n = 3,14159..... . Dit betekent dat een hoek van 1 graad 
overeenkomt met een hoek van 2n/360 of n/180 radialen. De trans- 
formatieformules voor de transformatie van ‘onze! coördinaten naar 
beeldschermcoördinaten hebben we op p.5 gegeven. De drie BASIC- 
opdrachten voor het berekenen van de beeldschermcoördinaten van 
het j-de hoekpunt van de n-hoek worden nu: 


W = (360/N)*n/180 : W1 = J*W 
X(J) = INT(U + A*COS(WI) + H) 
Y(J) = INT(V - BxSIN(WI) + H) 


Voor n gebruiken we in het programma de variabele PI die als 
waarde krijgt 4.ATN(I). ATN is de inverse tangensfunctie. Er 
geldt dat tan(45°)=1, dus tan(n/4)=1. Dit betekent dat de inverse 
tangens van 1 gelijk is aan n/4, dus atan(1) = n/4, hetgeen tot 
gevolg heeft dat n = 4.atan(1); in BASIC: PI=44ATN(1). We had- 
den ook PI=3.14159 kunnen nemen! 


*) HRG = Hoge Resolutie Graphics 
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In het navolgende programma komen we bovenstaande opdrachten 
tegen. Kiezen we voor A en B dezelfde waarden, dan krijgen we 
een regelmatige n-hoek met al zijn diagonalen. Bovendien is de 
ellips dan een cirkel geworden. 


ev 
AVS 


pe 


pe 
N 


100 REM programma 5 diagonalen in een 
n=hoek 

1:10: CS 

120 INPUT "halve grote as A<=128"; A 

130 INPUT "halve kleine as B<=96"; B 

140 INPUT "aantal hoekpunten N<25'';N 

150, CLS 

160 COLOR 1,15,4 : SCREEN 2 

170 DEF FNX(X)=INT(37+X/1.4+.5) 

180 DIM XCN),YON) z PI=Z4KATN(I) 

190 U=128:V=96:H=. 5: W=(360/N)*PI/180 

200 REM punten op de ellipsomtrek 
berekenen 

210 FOR J=0 TO N=1 


220 WI=J*W 

230 XO) =INTCUHARCOSCWIJ +H) 
240 YOJDZINT (V-BESINCWI)+H) 
250 NEXT J 


260 REM punten op de ellipsomtrek 
met elkaar verbinden 
270 FOR I=0 TO N=-2 


280 FOR J=I+1 TO N=-1 

290 LINE CFNXCXCIJ),YCIJD= 
CENXCXCJHD DY) 

300 NEXT J 

310 NEXT I 


320 AS=INPUT$(1) 
530: COLOR: 154,4 ‘= END 
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Programma 6 is wat veeleisender in die zin dat het wat voorberei- 
ding nodig heeft. 


We willen een reeks ingeschreven vierkanten zo tekenen dat de 
hoekpunten van een ingeschreven vierkant uit de reeks op de zij- 
den liggen van zijn voorganger. We zien deze situatie voor het eer- 
ste, tweede en derde vierkant van de reeks in de onderstaande 
tekening. De afstand tussen een hoekpunt van een vierkant en een 
hoekpunt van het volgende vierkant (S, en S>) is steeds een vast 
deel van de zijde waarop de hoekpunten liggen. Zo is Sj = ai/k en 
S2 = a2/k. In het algemeen is Sn = an/k, waarbij an de zijde is van 
het n-de vierkant in de reeks en k de verkleiningsfactor. De waar- 
de voor k kunnen we zelf kiezen. k=16 geeft een fraaie tekening. 


XX(3), YY (3) 
X(4),Y(4) a, X(3),Y(3) 
XX(4), YY(4) fi 
a, a 
a, a, 
Az Ay 
XX(2), YY(2) 

XD, Ya) a, X(2), Y(2) 


XX, YY) 


CKC ,YCIDD); (X(2),Y(2)); (X(3),Y(3)) en (X(4),Y(4)) zijn de 
coördinaten van de hoekpunten van een bepaald vierkant uit de 
reeks. Hoe berekenen we nu de coördinaten van de hoekpunten van 
het volgende ingeschreven vierkant ((XX(1),YY(1)); (XX(2),...)? 


Voor het eerste hoekpunt van het eerste ingeschreven vierkant geldt: 


Xx(I) =X(1) + EDE an 


YY) = YI) + DD 


Ga na dat dit klopt in bovenstaande tekening. 
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Voor het tweede hoekpunt van het ingeschreven vierkant geldt: 


vy(a) =v(2) «+ KOE 
In het algemeen geldt: 
KRIEK) EEDE) èn 


dod CD Bi dl DN. EDE) 


voor J = 1,2,3,4, waarbij X(5) = X(1) en 
Y(5) = YI). 


Als we het tweede vierkant getekend hebben, veranderen we 
XX(I) in X(1), YY(I) in YI), XX(2) in X(2), YY(2) in Y(2), 
enzovoorts, en tenslotte X(5) in X(1) en Y(5) in Y(1) en we gaan 
met dezelfde formules opnieuw XX(1),YY(1),... enz. berekenen, 
maar dan voor de hoekpunten van het volgende vierkant. We zien 
dit gebeuren in de regels 270 t/m 300 van het volgende programma. 
De regels 200 t/m 220 tekenen het vierkant, terwijl de regels 230 
t/m 260 de hoekpunten van het volgende vierkant berekenen. Een 
structuurdiagram voor het programma ziet er zo uit: 


begin programma ingeschreven vierkanten 
HRG-voorbereidingen 


lees waarde voor K in 
bepaal coördinaten voor het eerste vierkant 


voor het eerste t/m 40-ste vierkant doe 


teken dit vierkant 


bepaal coördinaten voor het volgende vierkant 


einde programma 
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Hier komt het programma: 


100 REM programma 6 ingeschreven 
vierkanten 

110 ELS 

120 INPUT "toets K in 1<K<20"'; K 

130 COLOR 1,15,4 : SCREEN 2 

140 DEM XC5I AV CHIEKK CDI 5 VVCS 

150 FOR J=1 TO 5 

160 READ X(J),Y(J) 

170 NEXT J 

180 GES 2 Hess 

190 FOR N=1 TO 40 


200 FOR J=1 TO 4 
210 LINE (XJ), YJ) = 
OKCIHIDD YH) 
220 NEXT J 
230 FOR J=1 TO 4 
240 XK EKD + 
INT (AX AID =X (JI) /K+H) 
250 YY EVD + 
INT (CAY J#1D =Y JD) /K+H) 
260 NEXT J 
270 FOR J=1 TO 4 
280 XCIDEKKAI) 2 YI =EYYCID 
290 NEXT J 
300 XC5IEXCI) 2 YS) =H) 
310 NEXT N 


320 AS=INPUT$(1) 

330 COLOR 15,4,4 : END 

340 DATA 60,192,196,192,196,0,60,0 
„60,192 


Re 


U kunt dit programma als uitgangspunt voor een uitgebreider pro- 
gramma gebruiken. Verdeel het beeldscherm in, bijvoorbeeld, twee 
rijen van vijf gelijke vierkanten (elk vierkant 50x50 puntjes). In elk 
van deze tien vierkanten tekent u, met bovenstaand programma, een 
reeks van ingeschreven vierkanten. Teken tien identieke reeksen, 
of kies steeds een andere waarde voor K of probeer ze ten opzichte 
van elkaar te laten draaien. 


2 Grafieken van functies in 
cartesische vorm 


In het eerste hoofdstuk hebben we alleen rechtlijnige patronen gete- 
kend, bijvoorbeeld het programma voor het tekenen van alle diago- 
nalen in een regelmatige n-hoek. 


In dit en het volgende hoofdstuk zullen we ons uitvoerig bezig hou- 
den met het tekenen van de grafiek van een aantal continue en niet- 
continue functies. In de wiskunde noemen we de grafiek van een 
niet-lineaire functie een kromme. De grafiek van een lineaire func- 
tie (bijvoorbeeld de functie y = 4x + 3) noemen we een rechte. We 
kunnen de vergelijking van zo'n kromme op drie manieren formule- 
ren: 


A ; met cartesische coördinaten : y = f(x) 
B ; met poolcoördinaten :r = f(@) 
C ; of in parametervorm :x=f(t), y =g(t) 


Niet-wiskundigen kennen vaak alleen de gewone (cartesische) vorm 
y=f(x). Een voorbeeld van een niet-lineaire functie in cartesische 
vorm is de functie y = 2x? - 3x + 4. De grafiek van deze functie 

is een parabool. 


De functie r = 110.cos(4wp) stelt een kromme in poolcoördinaten 
voor. Als we de hoek wp laten lopen van 1 tot 360 graden en we 
tekenen bij elke hoek p onder die hoek een punt op afstand r van 
de oorsprong, dan krijgen we een soort bloem als op pagina 33. 
Deze vorm r = 110.cos(4p) heet de poolcoördinatenvorm. Een punt 
in het platte vlak wordt nu niet gekenmerkt door de afstand van 
dat punt tot de x- en y-as (cartesisch), maar door de afstand tot 
de oorsprong (r) en de hoek p die de x-as maakt met de lijn die dat 
punt met de oorsprong verbindt. Poolcoördinaten komen in hoofd- 
stuk 3 aan de orde (zie tekening op p.16). 


Er is nog een manier om de vergelijking van een kromme weer te 
geven en dat is de parametervorm. Hierbij worden de x- en y- 
coördinaat van een punt op de kromme afhankelijk gemaakt van een 
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P(x1,y1) P(r,p) 


pool 


Hetzelfde punt P in een 
cartesisch coördinatenstelsel en een poolcoördinatenstelsel 


bepaalde parameter. Een voorbeeld zijn de vergelijkingen 
Xx =a.cos(t) en y =a.sin(t). 


Als we t laten lopen van 0 tot 360°, beschrijven de daarbij horende 

punten (x,y) precies de omtrek van een cirkel met straal a. De car- 
tesische vorm van deze cirkelvergelijking is x? + y? = a?, die velen 

direct zullen herkennen als de 'cirkelvergelijking'. 


In dit hoofdstuk houden we ons dus alleen bezig met het tekenen 
van grafieken van functies en relaties die door middel van cartesi- 
sche coördinaten beschreven worden. De andere twee vormen komen 
in hoofdstuk 3 aan de orde. 


Continue functies 


Een continue functie is een functie waarvan de grafiek in één vloei- 
ende beweging van ons 'potlood', dat wil zeggen zonder het potlood 
van het papier te hoeven halen, getekend kan worden. Er mogen 

'knikken' in voorkomen maar geen onderbrekingen (zie tekening op 


Del). 


De linkergrafiek kunnen we in één beweging, zonder het potlood 
van het papier te halen, tekenen; bij de rechter grafiek kan dat 
niet. De linker grafiek is de grafiek van een continue functie: de 
rechter van een niet-continue functie. Als we alleen links of alleen 
rechts van het punt xo kijken dan is de functie op die intervallen 
natuurlijk wel continu! De functie is alleen niet-continu in het punt 
X0. 
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X2 
Knikpunten 
+ 


Functie is niet-continu 
in het punt xo 


We willen nu een algemeen programma ontwikkelen dat, gegeven de 
grenzen a en b van een bepaald interval (asxsb) de grafiek tekent 
van een willekeurige continue functie y = f(x). Mochten de x-as en 
de y-as in het gebied liggen waarin de grafiek van de functie wordt 
getekend, dan willen we dat ook beide assen door het programma 
getekend worden. Om het programma kort, maar toch algemeen, te 
houden brengen we de volgende vereenvoudigingen aan. 


1. De vergelijkingen van de functie, die getekend moet worden, 
wordt in een subroutine, vanaf regel 1000, opgenomen. De 
functiewaarde y = f(x), bij een bepaalde x-waarde, wordt 
berekend door op dat moment met de gewenste x-waarde de 
subroutine aan te roepen (GOSUB 1000). Er wordt dus niet 
gebruik gemaakt van de methode waarin de functie als tekst 
(INPUT-opdracht met stringvariabele) wordt ingelezen waarna 
het programma deze tekst zelf omzet in een goede BASIC 
functiedefinitie. Een andere mogelijkheid zou zijn om de functie 
met een DEF FN-opdracht in het programma te definiëren. 


2. Op de coördinaatassen wordt niet automatisch een schaalverde- 
ling aangebracht en ook wordt er geen tekst bij afgedrukt. 
Via een omweg kan uw MSX-computer wel tekst op het hoge- 
resolutiescherm afdrukken. We geven hiervan in de appendix 
twee voorbeelden. Kijk vast in de appendix als u naast teke- 
ningen ook tekst en getalwaarden wilt afdrukken. Wilt u toch 
graag weten wat bij een bepaalde waarde van x de functie- 
waarde (y) is, dan kunt u bijvoorbeeld na het tekenen van de 
grafiek de computer (op de printer) een lijstje met x- en y- 
waarden laten afdrukken. 
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Terug naar het programma-ontwerp. In het eerste deel van het pro- 
gramma berekent de computer, na het inlezen van de waarden van 
de intervalgrenzen a en b, de grootste en de kleinste functiewaarde 
(y-waarde) in het opgegeven interval [a,b]. De waarde van HP 
(Hoogste Punt) en LP (Laagste Punt) worden op het beeldscherm 
afgedrukt. We weten zo in welk gebied de computer gaat tekenen. 
Hierna vraagt het programma of we HP nog groter en of we LP nog 
kleiner willen maken om daarmee een fraaiere tekening te krijgen. 
Zo niet, dan toetsen we voor HP en LP dezelfde waarden in als het 
programma heeft berekend; anders toetsen we de door ons gewenste 
maximale en minimale functiewaarden in. 


Het volgende programmadeel (regels 290-340) bevat de lus 

(FOR X=A TO B STEP DX) voor het tekenen van de grafiek. De 
coördinaten x,y van een punt op de grafiek moeten met de juiste 
transformatieformules omgezet worden in beeldschermcoördinaten 
Arke 


Om het interval [a,b] voor te kiezen x-waarden (asxsb) om te 
zetten in het HRG-interval [0,255] (Os Xs 255) gebruiken we de 
volgende evenredigheid: 


(x-a) : (b-a) = X : 255 ofwel Ei = En 


in BASIC: 
XX = INT(KX*(X-A)+0.5), met KX = 255/(B-A). 


Omdat BASIC in hoofdletters 'werkt' hebben we in deze formule 
voor x de variabele X genomen, voor X de variabele XX, voor a de 
variabele A en voor b de variabele B. 


We doen dit ook voor het afbeelden van het functiebereik LPS ys HP 
op het HRG-beeldbereik Os Ys 191: 


(HP-y) : (HP-LP) =Y : 191 ofwel Hedp= 11 


Dit geeft in BASIC: 
YY =INT(KY*(HP-Y)+0.5) met KY = 191/(HP-LP) 
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We zien deze berekeningen van XX en YY in de regels 310 en 320. 

In de volgende regels worden twee, op deze wijze berekende, pun- 
ten (X1,Y1) en (X2,Y2) door een recht lijntje met elkaar verbonden. 
Vervolgens (regels 350-380) wordt de ligging van de x- en y-as 
bepaald. Kiezen we in bovenstaande transformatievergelijkingen voor 
X en Y de waarde nul, dan vinden we respectievelijk de ligging van 
de y-as en van de x-as (regels 350-370). 


Met deze uitleg is de werking van het programma hopelijk te volgen. 


Test u het programma 7 met willekeurige, maar continue, functies. 
In dit voorbeeldprogramma hebben we de functie y = e-0. Ìx, cos(x) 
gekozen, hetgeen de grafiek van een gedempte trilling oplevert. 
Hieronder zijn enkele ideeën voor minder moeilijke functies, 


functie regel 1000 waarde voor a waarde voor b 
y = sin(x) Y = SIN(X) 0 2n (= 6,2832) 
v dinhe. 4 Y =XAZ2 -4 +4 

y =eX Y = EXP(X) -3 +3 
y=x3-2x?-x Y = XA3-2*XA 2-X -1 +3 


hf P 
Li y=e 0,1 
(zie regel 1010 programma 7 op p.20) 


"X cos(x) 
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100 REM programma 7 grafiek van een 
continue functie 

118 CIS 

120 INPUT "Linker intervalgrens "; A 

130 INPUT "rechter intervalgrens"; B 

140 IF A>B THEN SWAP A‚,B 

150 HP==100000!:LP=100000!: 

DX=(B=A)/64 
160 FOR X=A TO B STEP DX 


170 GOSUB 1000 'functiewaarde 
berekenen 

180 IF Y>HP THEN HP=Y 

190 IF Y<LP THEN LP=Y 

200 NEXT X 


210 PRINT "grootste Y=waarde:"; HP 

220 PRINT "kleinste Y=waarde:"; LP 

230 INPUT "bovengrens voor Y'"; HP 

240 INPUT "ondergrens voor Y"; LP 

250 CLS 

260 COLOR 1,15,4 7 SCREEN 2 

270 DEF FNX(X)=INT(37+X/1.4+.5) 

280 KX=255/(B-A):KY=191/(HP=LP) :H=. 5 

290 FOR X=A TO B STEP DX 

300 GOSUB 1000 'functiewaarde 
berekenen 


310 XXZINT(KX*CX=A) +H) 
320 YYZINTCKY*CHP=Y) +H) 
350 IF X=A 


THEN X1=XX : Y1=YY 
ELSE X2=XX : Y2=YY: 
ENE. GENKERID NY 19 
CENX(CX2) ,„ Ye): 
Xi=X2 © Yi=Yye 
340 NEXT X 
350 X1=0:Y1=INT(KY*HP+H) :X2=255:Y2=Y1 
360 IF Y1>=0 AND Y1<=191 THEN 
LINE CEN%CATDG YI CENKENZD, VL) 
370 X1=INT(KX*C=A)+H)EY1=0: 
x2=X1:Y2=191 
380 IF X1>=0 AND X1<=255 THEN 
LINE (ENXCX1) YI) =CFNX(X2) ,Y2) 
390 AS=INPUTS$S(1) 
400 COLOR 15,4,4 1 END 
1000 REM subroutine functiewaarde 
berekenen 
1010 Y=EXP(-. TX) COS CX) 
1020 RETURN 
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Nu volgen drie korte programma's. 


Programma 8 is de eerste van deze drie. Dit programma tekent tien 
in fase verschoven sinuskrommen, allemaal in hetzelfde coördinaten- 
stelsel. De vergelijking van het stelsel is 


y =sin(xtnp), met als fase p = 5 (20°). 


We vinden de tien vergelijkingen door voor n de waarden 0,1,2 t/m 
9 te kiezen. De lusvariabele J (regel 160) komt overeen met de 
beeldschermcoördinaat X2. Omdat we in dit programma de beeld- 
schermcoördinaat X2, dat wil zeggen J, in stappen van 5 naar 255 
laten lopen (regel 160), moeten we eerst de bij X2 behorende waar- 
de voor x berekenen. Dan pas kunnen we de functiewaarde y = f(x) 
berekenen. We kunnen de evenredigheid 


Cara): © (Dad KL 7 205 
natuurlijk ook gebruiken om x uit X2 te berekenen; hieruit volgt: 


„ (b-a) 
xXx = 55 2 + a. 


In het onderstaande programma zijn de gekozen intervalgrenzen 
a =0enb = 2r, zodat bovenstaande vergelijking wordt: 


_ 2m 
x= 755 XL. 


Als we C = JEE nemen, wordt in het programma het argument van 
de sinusfunctie dus C*X2+N*P, ofwel C+J+N*P. 


In het programma wordt de waarde C+#J als X berekend, dus zien we 
in regel 170 X+N*P als argument van de SINusfunctie. Natuurlijk 
moeten we de functiewaarde y = SIN(X+N*P) nog omzetten naar de 
beeldschermcoördinaat Y zodat we krijgen 


Y2 = INT(V-KxSIN(X+N*P)+H). 


Door voor de schaalconstante K de waarde 90 te kiezen (regel 140) 
krijgen we een mooie 'grote' tekening. 


Wie met kleuren wil werken kan dit programma uitbreiden door de 
diverse krommen andere kleuren te geven. Iets moeilijker zal het 
zijn de 'banen' tussen de sinuskrommen in te kleuren. Probeer het 
eens. 
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100 
110 
120 
130 
140 
150 
160 
170 


180 


190 
200 
210 
220 
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REM programma 8 sinuskrommen 
CLS z PI= 4XATN(1) 
COLOR: "1,154 SCREEN 2 
DEF FNX(X)=INT(37+X/1.4+.5) 
V=96:K=90O:H=. S:P=PI/I: C=2KPI/255 
FOR N=0 TO 9 
FOR J=0 TO 255 STEP 5 
X=JkC 2 
YZINT(V=K*SIN(X4AN*P) +H) 
IF J=0 
THEN X1=J:Y1=Y 
ELSE 
xèsJiYe2=Y: 
LINECFNX(XI) YI) -CENX(X2) ,„Y2) 
sxi=X2:Y1=Y2 
NEXT J 
NEXT N 
AS=INPUT$S(1) 
COLOR 15,4,4 : END 


MSX kleurcodes 


transparant 
zwart 
middelgroen 


lichtgroen 
donkerblauw 
lichtblauw 
donkerrood 
hemelsblauw 


middelrood 
lichtrood 
donkergeel 
lichtgeel 
donkergroen 
magenta 
grijs 

wit 


4 
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Programma 9 tekent een stelsel parabolen met de vergelijking 
vain et 


Alle parabolen (bepaalde t-waarden) hebben dezelfde snijpunten 
met de x-as (x = len x =-1). 


100 REM programma 9 paraboolstelsel 
110: GLS 

120: COLOR 11ste e= SCREEN 2 

130 DEF FNXCX)=INT(374X/1.4t.5) 

140 U=128:V=96:H=.5 

150 FOR K== 90 TO 90 STEP 10 


160 FoR X=e110 TO- 110 STEP 5 

170 XX=INTCU+X+H) 

180 Y==K&XX/6400FK: Y=INT(V-Y+H) 
190 EF X= =110 


THEN X1=XX:Y1=Y 
ELSE X2=XX:Y2=Y: 
LINE (FNXOI),Y1)= 
CENXCX2) ,Y2): 
Xl=X2:Y1=Y2 
200 NEXT X 
210 NEXT K 
220 AS=INPUTS(1) 
230 COLOR 15,4,4 : END 


Voeg toe: 165 IF X<=0 THEN KLEUR=10 
ELSE KLEUR=4 
en verander de LINE-opdracht in regel 190 in: 


LINECENX(X1),Y1)- 
CENX(X2) ,Y2) „KLEUR: 


Ook heel mooi wordt het met: 
165 IF INT((K/10)/2)*2=K/10 
THEN KLEUR=10 
ELSE KLEUR=5 


Soms willen we de oppervlakte tussen de grafiek van een functie 
en de x-as berekenen (de integraal bepalen) of laten tekenen. 
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Programma 10 kleurt zo'n oppervlak tussen de x-as en de grafiek 
van de functie 


y =cos(x) - coal) „eos(5x) cos(7x) 


Ook nu is de lusvariabele J de HRG-coördinaat X2. De waarde voor 
a is -r en die voor b is tr, waaruit volgt dat 
„J.2n 
“209 
2n 


Ook nu nemen we in het programma C = 55 


xX Eee 


100 REM programma 10 oppervlak onder 
een kromme 

110 GUUS 

120 COLOR 1,154 2 SEREEN 

130 DEF FNX(X)=INT(37+X/1.4+.5) 

140 V=96:H=.S:PI=4KATN(I) :K=80: 

C=2PI/255 

150 FOR J=0 TO 255 

160 X=JkC-PI 

170 GOSUB 1000 'functiewaarde 

berekenen 

180 YZINT (V=-K«*Y+H) 

190 LEENE CENX CIN IE GENN ESD 5 VI 

200 NEXT J 

210 AS=INPUTS(1) 

220 COLOR 157454 s END 

1000 REM subroutine functewaarde 

berekenen 
1010 Y=COS(XI-COS(E4XI/3+COS(SHKI/S5= 
COS(7*X)/7 
1020 RETURN 


Probeer dit programma zo te maken dat de drie oppervlakken in bij- 
voorbeeld rood, wit en blauw gekleurd worden. Gebruik hiertoe de 

mogelijkheid om in de LINE-opdracht de afdrukkleur van de lijn op 

te geven (zie p.23). 
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Niet-overal-continue functies 


2 
De functie y = en kan niet met behulp van programma 7 gete- 


kend worden. Zouden we bijvoorbeeld voor x het interval -5<x£5 
kiezen, dan liggen hierin de punten x = 3 en x = -2. Dit zijn de 
twee punten waarvoor de noemer van bovenstaande functie nul is. 
Als de computer bij het tekenen bij één van deze twee punten 
belandt, zal op het scherm een foutmelding (DIVISION BY ZERO) 
‘delen door nul' verschijnen en zal het programma afgebroken wor- 


den. 


Dit zal ook gebeuren bij logaritmische functies met een niet-positief 
argument of bij wortelfuncties met een negatief argument. Zelfs bij 
het tekenen van continue (nette) functies kunnen problemen optre- 
den. We komen hierbij in de problemen als de functiewaarden heel 
groot of heel klein worden. Transformatie van dergelijke waarden 
naar HRG-beeldschermeoördinaten (0-255 en 0-191) heeft dan geen 
enkele zin meer, omdat de aard van het verloop van de functie in 
het geheel niet meer tot uitdrukking komt. Het is daarom beter om 
de functiewaarden van een continue functie naar boven en beneden 
te begrenzen. Dit heeft tot gevolg dat de grafiek van een continue 
functie er uit zou kunnen zien als de grafiek van een niet-continue 
functie. 


Het zal duidelijk zijn dat je een algemeen programma voor het teke- 
nen van een willekeurige continue of niet-continue functie niet 
zomaar even opschrijft. Dergelijke programma's kom je in de vak- 
literatuur dan ook niet of nauwelijks tegen, en mocht je er wel een 
tegenkomen, dan betreft het òf een heel groot programma òf een 
programma dat voor een bepaalde functie, waarvan men van te 
voren weet waar de discontinuiteiten zitten, geschreven is. 


We geven nu een verbazend kort programma voor het tekenen van 
de grafiek van een willekeurige functie y = f(x). Veel wiskunde- 
leraren, scholieren en studenten zullen nu hun ‘oren! spitsen. U 
begrijpt dat, gezien het bovenstaande, hierbij wel enkele beperkin- 
gen gelden: 


l. Degene die het programma gaat gebruiken moet een beetje kun- 
nen programmeren. De functie moet namelijk in gedeelten in 
een subroutine (vanaf regel 1000) beschreven worden. 

2. De programmagebruiker moet voldoende wiskundige kennis 
bezitten om te kunnen bepalen voor welke x-waarden functies 
moeilijkheden kunnen opleveren. 


We zullen zien dat, normaal gesproken, slechts een paar BASIC- 

regels nodig zijn om de functiebeschrijving te programmeren. Het 
berekenen van nulpunten met behulp van een of ander numeriek- 
wiskundig algoritme is in het geheel niet nodig. 
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In het onderstaande programma 11 gebruiken we twee variabelen 
FZ en FA als zogeheten vlaggen (Flags). Een vlag is een variabele 
die slechts twee waarden (vaak 0 en 1) kan aannemen. 


Om de werking van de vlag FZ duidelijk te maken geven we hieron- 
der de subroutine 1000 met de functiebeschrijving van de functie 


x2+3 


VT Ex 


REM subroutine functiewaarde 
berekenen 
N=XkX= erb 
IF N=0 THEN FZ=1 
ELSE Y=(X*X+3) IN: 
IF Y<LP OR Y> HP 
THEN FZ=1 
ELSE FZ=0 
RETURN 


De vlag FZ wordt alleen op 1 gezet als het punt (x,y) niet op het 
scherm getekend kan worden. Dit is het geval als de functie niet- 
eontinu is in het punt (x,y) (N = 0; x = 3en x =-2) of als de 
functiewaarde buiten het opgegeven functiebereik (LP<y<HP) valt. 


Waarvoor dient nu de tweede vlag FA? Bekijk het volgende stukje 
programma eens: 


FA=1 
FOR X=A TO B STEP DX 
X2=INT(KX*CX=-A) +H) 
GOSUB 1000 'functiewaarde 
berekenen 
IF FZ=1 THEN FA=1 
ELSE IF FA=Í 

THEN X1=X2:FA=0: 
YAZINT(KY*CHP=Y) +H) 

ELSE Y2=INT(KY*CHP=-Y)+H) 2 
LINECENXCX1, YI) = 
CEFNX(X2) „Yeal ext=Xea:Yi=Y2 

NEXT X 


Als FZ = 1 dan wordt FA ook gelijk aan 1 gemaakt en wordt de vol- 
gende X-waarde bekeken (FOR-lus) zonder dat het nieuw berekende 
punt met het laatstgetekende punt verbonden wordt. Is FZ echter 

0 (berekende punt ligt in het beeldvlak) dan wordt onderzocht of 
FA daarvoor soms op 1 gezet is. Is dit zo (tweede THEN) dan 

moet het nieuw berekende punt als nieuw beginpunt (X1,Y1) geko- 
zen worden en dit mag dan ook niet met het vorige getekende punt 
verbonden worden. FA wordt in dit geval nul gemaakt. 
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Dus als zowel FZ als FA nul zijn, wordt het nieuw berekende punt 
verbonden met het laatstgetekende punt. 


Om te zorgen dat het programma probleemloos werkt moet voor het bere- 
kenen van het eerste punt van de grafiek, dus voor de FOR-lus, de 
vlag FA op 1 gezet worden, waardoor we er zeker van zijn dat de waar- 
den Xl en Y1 berekend worden. Hier volgt het volledige programma: 


100 REM programma 11 grafiek van een 
willekeurige functie 


110 ELS 

120 INPUT "linkergrens voor X "; A 
130 INPUT "rechtergrens voor X "; B 
140 INPUT "ondergrens voor Y me LP 
150 INPUT "bovengrens voor Y EH 


160 IF A>B THEN SWAP A,B 
170 KX=255/(B-A):KY=191/ (HP=LP) :H=.5 
180 DX=(B-A)/255 
190 CLS 
200 COLOR 1,15,4 : SCREEN 2 
210 DEF FNX(X)=INT(37+4+X/1.4+.5) 
220 FA=1 
230 FOR X=A TO B STEP DX 
240 X2=INT(KX*(X-A)+H) 
250 GOSUB 1000 'functiewaarde 
berekenen 
260 IF FZ=1 THEN FA=1 
ELSE IF FA=1 
THEN X1=X2:FA=0: 
YAZINT(KY*CHP=Y)+H) 
ELSE Y2=INT(KY*(HP=Y)+H): 
LINECENXCX1) „Y1)= 
CENX(X2) „Y2)sxi=X2:yYi=Y2 
270 NEXT X 
280 X1=0:Y1=ZINT(KY*HP+H) :X2=255:Y2=Y1 
290 IF Y1>0 AND Y1<=191 THEN 
LINE CFNX(X1), YI) -CENX(X2) ,Y2) 
300 X1=INT(KX*(=A)+H):Y1=01X2=X1: 
y2=191 
310 IF X1>0 AND X1<=255 THEN 
LENE GENKCKTI YTS CENKCR2I YAI 
320 AS=INPUTS$S(1) 
530 COLOR 19,44 © END 
1000 REM subroutine functiewaarde 
berekenen 
1010 N=XkXeX=6 
1020 IF N=0 THEN FZ=1 
ELSE Y=(X&X+3) /N: 
IF Y<LP OR Y> HP 
THEN FZ=1 
ELSE FZ=0 
1030 RETURN 


28 Grafieken van functies in cartesische vorm 


Dit programma tekent de grafiek van de functie 
2 
y= EN ; kies voor A, B, LP en HP respectievelijk 
nm -5, 5, -10 en 10 


Wilt u een andere functie tekenen, bijvoorbeeld y = In(x?-2), her- 
schrijf dan subroutine 1000 als volgt: 


1000 REM subroutine functiewaarde 
berekenen 
1010 N=X*X=-2 
1020 IF N<=0 THEN FZ=1 
ELSE Y=LOG(N): 
ER YEP OR YS KR 
THEN FZ=1 
ELSE FZ=0 
1030 RETURN 


U krijgt dan deze grafiek: 


Kies voor A, B, LP en HP de waarden -25, 25, -5 en 8. 


Niet-overal-continue functies 


Hier volgt een hele fraaie: 


y=3-xt+In | 2 


Subroutine 1000 wordt nu: 


1000 REM subroutine functiewaarde 
berekenen 
1010 N=X=3 
1020 IF N=0 THEN FZ=1 
ELSE Y=3-X+LOG(ABS((X-1) 


PND 
LF YSLP OR TO HP 
THEN FZ=1 
ELSE FZ=0 


1030 RETURN 


Kies voor A, B, LP en HP de waarden -5, 10, -4 en 8. 


29 


3 Krommen in poolcoördinaten 
en in parametervorm 


In het vorige hoofdstuk hebben we ons beziggehouden met het teke- 
nen van de grafiek van een continue of niet-continue functie, waar- 
van de vergelijking als y = f(x) geschreven kon worden; dus in 
cartesische vorm. Nu gaan we grafieken tekenen van functies waar- 
van de vergelijking in poolcoördinaten geschreven wordt of in de 
zogeheten parametervorm. Dit levert zeer fraaie ‘beelden! op. 


Krommen met poolcoördinaten 


Functies, waarvan de grafiek een gesloten kromme laat zien, zijn 
vaak eenvoudiger met poolcoördinaten te beschrijven dan in carte- 
sische vorm. Als voorbeeld noemen we de 'hartkromme' (kardioide), 
waarvan de vergelijking in poolcoördinaten eenvoudig is, namelijk: 


r = k(1l + ecosp), 
maar waarvan de cartesische vorm nogal ingewikkeld is, namelijk 


(x?-y?-kx)? í 
k?(x7+y7) ’ 


Om de volgende programma's te doorgronden (niet om ze te draaien!) 
is een beetje theorie nodig. Het poolcoördinatenstelsel wordt in de 
wiskundelessen op school niet of nauwelijks behandeld. Eigenlijk is 
dit jammer, want hiermee (en met de parametervorm) kunnen juist 

de mooiste functies (en daarmee de fraaiste grafieken) beschreven 
en getekend worden. 


In een poolcoördinatenstelsel wordt elk punt in het platte vlak met 
twee coördinaten, te weten r en pw, bepaald. r is de afstand tussen 
het punt en de oorsprong en w (phi, spreek uit: fie) is de hoek 
tussen de horizontale lijn door de oorsprong en de lijn door de oor- 
sprong en het punt (zie tekening op p.31). 
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—- nul-richting 


oorsprong 


De oorsprong van het poolcoördinatenstelsel leggen we in het 
middelpunt van het beeldscherm (U = 128, V = 96). Zoals we weten 
ligt de oorsprong van het HRG-beeldschermcoördinatenstelsel in de 
linkerbovenhoek van het beeldscherm (HOME-positie). De nul-rich- 
ting in ons poolcoördinatenstelsel is horizontaal en wijst naar rechts. 
p=90° (nr /2 radialen) is verticaal naar boven; p=180° (mn radialen) is 
horizontaal naar links en w=270° (3n radialen) verticaal naar beneden. 


We gebruiken de volgende twee transformatieformules om het, uit 
de functievergelijking berekende, punt P(r‚p) om te zetten in het 
beeldschermpunt P(X,Y): 


INT(U + RxCOS(P) + 0.5) èn 
INT(V - R*SIN(P) + 0.5) 


X 
x 


In onderstaande figuur zien we hiervoor de verklaring. 


(0,0) 128 255 % 


Y 
P(X,Y) 

1 

1 

nr sing 

' 


_ a mn en 


ze 
rcosp 


De hoek wp, in de programma's voorgesteld door de variabele P, 
doorloopt steeds in positieve richting (tegen de klok in) het inter- 
val van 0 tot 2m. Dit is een interval in radialen (0-360 in graden). 


In alle programma's is als lusvariabele de hoek W in graden gekozen. 
Met de formule P = W.(7/180) (in de programma's P = W*RD) wordt 
de hoek in radialen berekend. Het voordeel van een lusvariabele 
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die het interval 0-360° doorloopt is dat de stapgrootte waarmee de 
lusvariabele steeds wordt opgehoogd een geheel getal kan zijn (1 

bijvoorbeeld) en dat hierdoor de programma's beter leesbaar zijn. 
Een neveneffect is dat alle programma's bijna woordelijk in Pascal 
zijn over te zetten; Pascal kent immers geen gebroken stapgrootte 
in een lus. 


Meer theorie hebben we niet nodig. De volgende programma's spre- 
ken grotendeels voor zichzelf. 


Deze figuur ontstaat als we in programma 12 (zie pagina hiernaast) 
voor N de waarde 14 nemen. 
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Programma 12 tekent de grafiek van de functie 
r = keos(no). 


In het programma kiezen we k = 95 (K) en n = 4 (N). 


100 REM programma 12 grafiek van de 
functie R=K*cos(4*Phij 

110: BLR 

120: COLOR 17154: E SCREEN (2 

130 DEF FNX(X)=INT(37+4+X/1.4+.5) 

140 PI=4XATN(I) iN=4 

150 U=128:V=96:H=. 5:RD=PI/180:K=I5 

160 FOR W=0 TO 360 STEP 3 


170 P=W*RD 

180 GOSUB 1000 'functiewaarde 
berekenen 

190 X=INT (U+R*COS(P)+H) 

200 Y=INT (V=-R*SIN(P)+H) 


210 IF P=0 
THEN X1=X:Y1=Y 
ELSE Xa Yez: 
LINE GENNER NIS 
CFNRUKEIN VES 
Xi=X2:Yl=Y2 
220 NEXT W 
230 AS=INPUTS(1) 
240 COLOR 15,4,4 : END 
1000 REM subroutine functiewaarde 
berekenen 
1010 R=EK*COS(N*P) 
1020 RETURN 


Het programma tekent een achtdelige draaisymmetrische figuur; een 
bloem met acht blaadjes. Als u met dit programma gaat experimen- 
teren zult u snel ontdekken dat voor even waarden van n een 2n- 
bladerige en voor oneven waarden van n een n-bladerige bloem ont- 
staat. Kunt u een tweekleurige bloem maken? 
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Wilt u dat uw computer sneller tekent? Neem dan een stapgrootte 
van 3° of van 5° in plaats van 1°. Bedenk dan wel dat de blaadjes 
wat hoekiger worden. 


Als u regel 1010 R=K4COS(N*P) vervangt door één van de 
onderstaande functies, krijgt u steeds een andere mooie figuur: 


R = COS(4xSIN(2*P)) : R = ABS(R) 
R = COS(4xSIN(3*P)) : R = ABS(R) 
R = SIN(3*SIN(24P)) : R = ABS(R) 
R = SIN(5*COS(2*P)) : R = ABS(R) 


U kunt naar hartelust trigonometrische functies met verschillende 
parameters 'mengen'. Laat uw creativiteit en nieuwsgierigheid de 
vrije loop. Probeer eens: 


1010 R=EK*SIN(S*COS(2*SIN(3*COS(4*P)))) :R=ABS CR) 


Omdat de poolcoördinaat r per definitie niet negatief mag zijn, moe- 
ten we in regel 1010 ook de opdracht R = ABS(R) opnemen. 


Nog veel mooiere figuren krijgen we als we niet één grafiek, maar 
een stelsel van gelijkvormige grafieken tekenen. 
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Programma 13 tekent het stelsel krommen: 


R = K*COS(4#SIN(2*P)). 


De parameter K loopt van 20 tot 100 in stappen van 10. 


100 


110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 


230 
240 
250 
260 


REM proramma 13 grafiek van de 
functie R=K*cos(4*sin(2*Phi)) 
CLS 
COLOR: 1,154 £ SCREEN 2 
DEF FNX(X)=INT(37+X/1.4+.5) 
PI=4&kATN(I) 
U=128:V=96:H=. S:RD=PI/180 
FOR K=20 TO 100 STEP 10 
FOR W=0 TO 360 STEP 2 
P=W*RD 
GOSUB 1000 ' R berekenen 
X=ZINTCU+RKCOS(P) +H) 
Y=INT(V=-R&SIN(P)+H) 
IF P=0 THEN X1=X:Y1=Y 
ELSE Hes EY A= Es 
GENE GENWEKTI NIS 
CENX(X2) ,Y2): 
xi=X2:Y1=Y2 
NEXT W 
NEXT K 
AS=INPUT$(1) 
COLOR 15,4,4 : END 


1000 REM subroutine R berekenen 


1010 


REKKCOS(4*SIN(2*P)) 


1020 RETURN 
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Programma 14 tekent sinuskrommen die cirkelvormig gekromd zijn. 


100 REM programma 14 sinuskrommen in 
cirkels 
1140 CLS 
120 COLOR 1,15,4 7 SCREEN è 
130 DEF FNX(X)=INT(37+X/1.4+.5) 
140 PI=Z4&KATN(I) 
150 U=128:V=96:H=. S:RD=PI/180 
160 FOR K==40 TO 40 STEP 10 
170 FOR W=0 TO 360 STEP 2 
180 P=W«&RD 
190 GOSUB 1000 ' R berekenen 
200 X=INT(U+R*COS(P) +H) 
210 Y=INT(V=-R&SIN(P)+H) 
220 IF P=0 THEN X1=X:Y1=Y 
ELSE X2=X:Y2=Y: 
LINECENX(X1) ,Y1)= 
CFNX(X2) „Y2): 
xi=X2:Y1=Y2 
230 NEXT W 
240 NEXT K 
250 AS=INPUT$S(1) 
260 COLOR 15,4,4 : END 
1000 REM subroutine R berekenen 
1010 R=60 + K*SIN(4*P) 
1020 RETURN 


Probeert u eens 
1010 R = 2#TAN(2*P)+KKSINC2KCOSC(SINC(6P))) 


voor een schitterend spinneweb. De variaties zijn echt onbegrensd! 
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Programma 15 tekent een bloemvormige figuur met blaadjes die in 
het midden aan steeltjes vastzitten. In het programma is N = 4 geko- 
zen. U kunt gerust andere waarden voor N proberen. Het effect 

is hetzelfde als bij programma 12. 


Ne 
GN 


100 REM programma 15 bloemen WE 


110 GES 
120 COLOR 1,15,4 2 SCREEN 2 
130 DEF FNX(X)=INT(37+X/1.4+.5) 
140 PI=4*ATN(I) iN=4:C=.25 
150 U=128:V=96:H=. S:RD=PI/180 
160 FOR K=30 TO 80 STEP 10 
1:20 FOR W=0 TO 360 STEP 3 
180 P=W&RD 
190 GOSUB 1000 ' R berekenen 
200 X=INTCU4+RKCOS(P)+H) 
210 Y=zINT (V=-R&SIN(P)+H) 
220 IF P=0 THEN X1=X:Y1=Y 
ELSE Xe2=X:Y2=Y: 
LINECFNX(X1) ,Y1)= 
GENKCK2N ) VDE 
xt=X2:Y1=Y2 
230 NEXT W 
240 NEXT K 
250 R=30:P1=180/NERD 
260 FOR J=1 TO N 


270 X1=ZINTCUER*COS(P) +H) 

280 Y1=INT (V=R&SINCP)+H) 

290 X2=INTCU+R*COS(P) +H) 

300 Y2=INT (V=R&SINCP)+H) 

310 LINE (FNX(X1) „YTD =(FNX(X2),Y2) 
320 NEXT J 


330 AS=INPUT$S(1) 

340 COLOR 15,4,4 : END 

1000 REM subroutine R berekenen 
1010 R=EK*(T+CHABS(SIN(N*P))) 
1020 RETURN 
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U kunt na elke doorloop van de K-lus de afdrukkleur veranderen. 
Ook de stelen van de bloem kunnen een eigen kleur krijgen. Tekent 
u een aantal van dergelijke bloemen naast of onder elkaar dan hebt 


u een begin gemaakt met 'computer-art', 
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Spiralen zijn altijd geliefde figuren. Programma 16 tekent logarit- 
mische spiralen of spiralen van Archimedes. Deze laatste soort spi- 
ralen hebben als vergelijking 

r=c.p 


In het programma op p.40 is voor C de waarde 3 gekozen. Kies 
gerust andere waarden, maar wel tussen 0,5 en 20. 


De logaritmische spiraal heeft (naar Bernoulli) de vergelijking 
r = ke°®, 


In het programma dat de onderste spiraal getekend heeft hebben 
we voor K de waarde 110 en voor C de waarde -0,2 gekozen. 


r=2p 


r=110e 029 


Omdat de logaritmische spiraal zich oneindig vaak rond de oor- 
sprong zal winden moet ervoor gezorgd worden dat het programma 
na bepaalde 'tijd' wordt afgebroken. 
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100 REM programma 16 spiralen 
110 CLS 
120 COLOR 1,15,4 : SCREEN 2 
130 DEF FNX(X)=INT(37+X/1.4+.5) 
140 PI=4XATN(I) :C=3 
150 U=128:V=96:H=. S:RD=PI/180 
160 FOR wW=0 TO 10000 STEP 3 
170 P=W&RD 
180 GOSUB 1000 ' R berekenen 
190 X=INTCU4R*COS(P)+H) 
200 Y=INTCV=R&*SIN(P)+H) 
210 IF X<0 OR X>255 OR Y<O OR Y>191 
THEN AS=INPUT$S(1) : STOP 
220 IF P=0 THEN X1=X:Y1=Y 
ELSE X2=X:Y2=Y: 
LINECENX(X1) ,Y1)=- 
CFNXCX2) „Y2) z 
xi=X2:Y1=Y2 
230 NEXT wW 
240 AS=INPUTS(1) 
250 COLOR 154,4,4 : END 
1000 REM subroutine R berekenen 
1010 R=C*P 
1020 RETURN 


Voeg toe regel 235 U = 133 : GOTO 160 
en zet in regel 210 achter THEN GOTO 235 in plaats van STOP 
en zie hoe dit de figuur verfraait. 


Voor het tekenen van een logaritmische spiraal maken we C gelijk 
aan -0,2 en herschrijven subroutine 1000 als volgt: 


1010 R=110*EXP(C*P) 
1020 IF R<5 

THEN AS=INPUT$(1) : STOP 
1030 RETURN 
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Als laatste programma met poolcoördinaten geven we het programma 
17 voor het tekenen van een willekeurige, in poolcoördinaten gefor- 
muleerde, continue of niet-continue functie. De programmastructuur 
komt overeen met de structuur van programma 11 (zie p.27). U kunt 
daar de werking van de vlaggen FZ en FA nog eens bestuderen. 
Ook de andere variabelen hebben dezelfde betekenis als hun naam- 
genoten in programma 11. Als voorbeeld in het programma hebben 
we een vrij ingewikkelde functie, die niet overal continu is, geko- 
zen: 


t = sin(1,5.p) 
1-2.cosp 


De figuur is getekend met A =-2, B = 2, LP =-2, HP = 2, WO = 0° 
en WN = 7209, 


Wilt u een andere functie proberen, dan hoeft u alleen de regels 
1010 en 1020 te veranderen door hierin de nieuwe functie R te bere- 
kenen. Mocht er iets fout gaan, dan kan het zijn dat u de volgende 
regel moet toevoegen: 1025 IF R<O THEN FZ=1 


ea 


3 ER 4.sin(1,5p+2) 
Deze zeldzame vlinder heeft als vergelijking R = en 
cos(p)(1+ SSEP) 


Voor A, B, LP, HP, WO en WN nemen we in programma 17: 
4, +4, -4, +4, O en 720. 
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100 REM programma 17 grafiek van de 
functie R=f(PHI) 

110: GES 

120 INPUT "linkergrens voor X 

130 INPUT "rechtergrens voor X ''; 

140 INPUT "ondergrens voor Y Ae 


ER 

id 

id 
150 INPUT "bovengrens voor Y We NP 

, 

‚ 


160 INPUT "startwaarde voor PHI'' 
170 INPUT "eindwaarde voor PHI "'; 
180 PI=4XATN(1) 
190 KX=255/(B-A):KY=191/ CHP=-LP) :H=. 5 
200 RD=PI/180 
210 LS 
eel COLOR 1,15 74 1 SCREEN 2 
230 DEF FNX(X)=INT(37+X/1.4+.5) 
240 FA=1 
250 FOR W=WO TO WN 
260 P=W*RD 
270 X2=INT(KX*(X=A)+H) 
280 GOSUB 1000 'functiewaarde 
berekenen 
290 IF FZ=1 THEN FA=1 
ELSE IF FA=1 

THEN X1=X2:FA=0: 
YIZINT(KY*CHP=Y)+H) 

ELSE Y2=INT(KY&*CHP=-Y)D+H): 
LINE(FNX(X1),Y1)= 
CENX(X2) „Vads t=2as vise 

300 NEXT W 
310 AS=INPUT$(1) 
320 COLOR 15,4,4 s END 
1000 REM subroutine functiewaarde 
berekenen 
1010 N=1-2*COS(P) 
1020 IF N=0 
THEN FZ=1 
ELSE R=SIN(3*P/2) /N 
:XERKCOS(P) : YERASIN(P): 
IF X<A OR X>B OR Y<LP OR 
Y>HP THEN FZ=1 
ELSE FZ=0 
1030 RETURN 


De parametervorm 43 


De parametervorm 


De meest interessante krommen krijgen we bij functies waarvan de 
vergelijking in parametervorm wordt opgesteld. Dit houdt in dat 
zowel de X- als de Y-coördinaat als functie van dezelfde, derde, 
parameter t worden uitgedrukt. In de natuurkunde zien we vaak 
de tijd als parameter (vandaar de t). In de wiskunde is de parame- 
ter t bijna altijd een hoek in radialen. 


Zo is de parametervorm van een cirkel met straal r en het middel- 
punt in de oorsprong: 


Xx = rcost en y = rsint 


Uit deze twee vergelijkingen kan gemakkelijk de cartesische vorm 
x?+y?=r? gedistilleerd worden. De vergelijking van een ellips met 
het middelpunt in de oorsprong en met een halve lange as a en een 
halve korte as b is: 


Xx = acost en y = b sint 


Het is eenvoudig programma's te ontwikkelen voor het tekenen van 
stelsels concentrische of excentrische cirkels. Ook stelsels ellipsen 
zijn, dankzij de parametervorm, eenvoudig te tekenen. We doen dit 
echter niet, omdat de figuren niet zo bijzonder zijn en vrij bekend. 
In hoofdstuk 1 hebben we trouwens al een ellips met behulp van de 
parametervorm geprogrammeerd. 


De volgende programma's tekenen figuren die u mogelijk nog niet 
kent en die zich voor computer-graphics bijzonder goed lenen. 
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Programma 18 tekent een Lissajousfiguur. In het algemeen is de 
parametervorm voor een dergelijke figuur: 


x= k‚-sin(f,-t+p‚) + kj-cos(f‚.t) 


y= kg-sin(fg.t+p3) E kg-cos(f‚-t) 


Hierin is p de fase, f de frequentie en k de amplitude. Om een mooie 
figuur te krijgen hebben we de volgende waarden gekozen: 


k,=k3=100, £,=16, p‚=0, fo=17, pP3=35 en ko=k, =f=f,=0 


3 24 2 4 


100 REM programma 18 Lissajousfiguren 
110 CLS 

120 PRINT "toets K1,F1,P1,K2,F2 in"; 
130 ENPUT K1,F1,P1,K2,Fè 

140 PRINT:PRINT 

150 PRINT "toets K3,F3,P3,K4,F4 in"; 
160 INPUT K3,F3,P3,K4, F4 

170 PI=4XATN(1) 

180 U=128:V=96:H=. S:RD=PI/180 

190 ELS 

200- COLOR 1,,15,4 2: SCREEN -2 

210 DEF FNX(X)=INT(37+X/1.4+.5) 

220 FOR W=0 TO 360 


230 T=WeERD 

240 GOSUB 1000 'x en y berekenen 
250 X=INTCU+X+H) 

260 Y=INT(V=Y+H) 


270 IF W=0 
THEN X1=X:Y1=Y 
ELSE X2=X:Y2=Y: 
LINECFNX(X1) ,„Y1)= 
CENX(X2) ,Y2): 
Xi=X2:Y1=Y2 
280 NEXT W 
290 AS=INPUT$S(1) 
300 COLOR 15,4,4 : END 
1000 REM subroutine x en y berekenen 
1010 XEKT&SINCFIAT+PI)+K2KCOS(F24T) 
1020 Y=K3*SIN(F3T+P3I)+K4KCOS (FAT) 
1030 RETURN 


Experimenteer met het programma! Zoek zelf de juiste waarden voor 
de constanten zodat u mooie figuren krijgt. Natuurkundigen zullen 
u hierbij graag helpen; zij weten welke waarden u moet nemen! 
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Programma 19 presenteert een niet alledaagse figuur. Deze figuur 
wordt dikwijls de 'vliegekopkromme' genoemd. Om de vliegekop hori- 
zontaal op het beeldscherm te krijgen moet de parameter t het inter- 
val 90°-450° doorlopen. 


100 REM programma 19 vliegekopfiguur 
110 CLS 

120 COLOR 1,15,4 : SCREEN 2 

130 DEF FNX(X)=INT(37+X/1.4+.5) 

140 PI=4XATN(I) 

150 U=128:V=96:K=30:H=. 5: RD=EPI/180 
160 FOR W=90 TO 450 STEP 3 


170 T=W&RD 

180 GOSUB 1000 'x en y berekenen 
190 X=INTCU+*X+H) 

200 Y=INT(V=Y+H) 


210 IF W=90 
THEN X1=X:Y1=Y 
ELSE X2=X:Y2=Y:t 
LINECFNX(X1) ,Y1)= 
CENX(X2) „Y2): 
xi=X2:Y1=Y2 
220 NEXT wW 
230 AS=INPUT$(1) 
240 COLOR 15,4,4 : END 
1000 REM subroutine x en y berekenen 
1010 X=EK&SIN(24*T)*(2.5+COS(34T)) 
1020 Y=Kk*2xCOS(3*T) 
1030 RETURN 


Het tekenen van de x- en y-as is achterwege gelaten om de 'kop' 
beter te laten uitkomen. 


Voor een echte vliegekop met 'ogen' moet u in de appendix kijken. 


Hoe geraffineerder u de functies x = f(t) en y = f(t) kiest, des te 
ongewoner worden de figuren. Probeer hier hoe creatief u kunt zijn. 


De 'huiswiskundigen' bij computerfabrikanten hebben vaak hun 
eigen lievelingsfuncties. Die zie je dan ook vaak in een demonstra- 
tieprogramma optreden. 
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Programma 20 tekent een stelsel krommen. Hewlett-Packard publi- 
ceerde deze tekening enige tijd geleden. Onder het programma staat 
een tabel met waarden voor A en B, die heel mooie figuren opleve- 
ren. Probeer zelf andere combinaties. U zult verrukt zijn over wat 
u ziet! 


100 REM programma 20 vlinders 
110 CLS 
120 PI=Z4XATN(I) 
130 PI=Z4&*ATN(I) 
140 U=128:V=92:H=.5:RD=PI/180 
150 KX=10:KY=10 
160 INPUT "toets A en B in"; A,„B 
10 GES 
180 COLOR 1,15,4 : SCREEN 2 toets in: -6,1 
190 DEF FNX(X)=INT(37+X/1.4+.5) 
200 FOR N==3 TO 3 
210 FOR W=0 TO 360 STEP 3 
220 T=W&RD 
230 X=(A+B)*COS(T)= 
N*B*COS ((A+B) /B*T) 
240 Y=(A+B)&SIN(T)= 
N*B*SIN((A+B) /B*T) 


250 XX=INTCUÊKX*X+H) 
260 YYZINT(V=KY*Y+H) 
270 IF W=0 


THEN X1=XX:Y1=YY 
ELSE X2=XX:Y2=YY: 
LINECFNX(X1) „Y1)=- 
CFNX(X2) ,Y2): 
xi=X2:Y1=Y2 
280 NEXT W 
290 NEXT N 
300 AS=INPUTS$(1) 
310 COLOR 15,4,4 : END 


Tabel voor mooie figuren 


Kies in deze programma's de stapgrootte voor hoek W (FOR W=..... 
gerust 1; dit duurt wat langer maar geeft mooiere tekeningen. 
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Het laatste programma uit dit hoofdstuk tekent 'supersymmetrische' 
figuren. Ook nu geven we een tabel met waarden voor de parame- 
ters A, B en C. Het is ongelooflijk hoe de figuur totaal verandert 
als we andere waarden voor één of meer parameters kiezen. Het idee 
voor programma 21 is afkomstig uit het Amerikaanse tijdschrift 
Creative Computing. Het kan soms even duren voor het tekenen 
wordt hervat. Druk niet direct op een toets, want dan wordt na het 
tekenen direct teruggeschakeld naar het gewone scherm en bent u 
uw tekening kwijt. Druk pas op een toets als u zeker weet dat de 
tekening af is. 


A=6/B=6/C=4 
toets in: 6,6,4 


A=2/B=71/C=3 
toets in: 2,7,3 


100 REM programma 21 symmetrische 
figuren 
110 CLS 
120 PI=4KATN(I) 
130 U=128:V=92:H=. 5:RD=EPI/180:K=100 
140 INPUT "toets A‚,B en C in"; A,B,C 
450 CLS 
160 COLOR 1,15,4 £ SCREEN 2 
170 DEF FNX(X)=INT(37+X/1.4+.5) 
180 FOR W=0 TO 360 
190 TEWERDER=EKKSINCCHT) 
200 X2=INTCU4RKCOS CAT) +H) 
210 Y2=INT(V=-R*SIN(B&T)+H) 
220 IF W=0 
THEN X1=X2:Y1=Y2 
ELSE LINECFNXCX1),Y1)= 
CENX(X2) ,Y2): 
x1=X2:Y1=Y2 
230 NEXT W 
240 AS=INPUTS(1) 
250 COLOR 19744 £ END 


Kijk in de appendix hoe u de waarden van A, B en C links naast de 
figuur op het graphic-scherm kunt afdrukken. 
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Tabel voor mooie figuren 


Het zou jammer zijn als u bij het zien van een mooie figuur de waar- 
den voor A, B en C niet hebt genoteerd. Verander het programma 

zo, dat na afloop van het tekenen de waarden van A‚ B en C op het 
beeldscherm of op de printer worden afgedrukt. Ook heel leuk zijn 


de combinaties: 


A 20 -40 100 
B sl -40 -0,5 
C 3 10 3 


A=4/B=6/C=1 


ee ee 
Te 
A=-1/B=-40/C=0,5 


4 Tekenen van driedimensionale 
figuren 


In dit en het volgende hoofdstuk gaan we ons bezighouden met het 
tekenen van driedimensionale lichamen (zo heet dat in de wiskunde) 
zoals kubussen, prisma's, pyramiden, kegels en bollen, en met het 
tekenen van driedimensionale grafieken van functies. Zo'n drie- 
dimensionale grafiek is een vlak in de ruimte met een vergelijking 
van de vorm z = f(x,y). Zo,krijgen we oe hel es ipv Mexicaanse 
hoed' (zie p.73) als we de functie z = e” Ess ) tekenen (e is het 
grondtal van de natuurlijke logaritme; e = 2,71828...). 


Er bestaan veel programma's waarmee driedimensionale figuren 
getekend kunnen worden. Waarom we toch ook in dit boek derge- 
lijke programma's laten zien komt, omdat er aan de meeste van deze 
programma's drie nadelen kleven: 


1. De meeste programma's voor driedimensionaal-tekenen zijn 
voor een bepaald graphics-systeem ontworpen en zijn slechts 
met veel inspanning geschikt te maken voor andere systemen. 


2. De wiskundige theorie die aan het driedimensionaal-tekenen 
ten grondslag ligt wordt zelden of zeer summier behandeld. 


3. Veel programma's zijn uiterst ingenieus ontworpen en draaien 
op de kleine microcomputers, in BASIC, heel langzaam. 


Met de volgende programma's en stukjes theorie hopen we deze drie 
nadelen uit de weg te ruimen. 


Er bestaan verschillende methoden om een driedimensionaal lichaam, 
bijvoorbeeld een kubus, in een plat vlak te projecteren. Wij hante- 
ren de projectiemethode, die u wellicht op school gebruikt hebt (of 
nog gebruikt). Stelt u zich een doorzichtige kubus voor met ribben 
van draadijzer. U staat voor deze kubus een beetje rechts van het 
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midden en kijkt er zo'n beetje schuin bovenop. Uw gezichtsstralen 
projecteren elke hoek, met de daaraan verbonden ribben, van de 
kubus in een, achter de kubus liggend, projectievlak. Zo ontstaat 
de bekende ‘schuine! afbeelding van een kubus. 


De verticale en horizontale ribben worden op ware grootte getekend. 
De ribben die naar achteren lopen worden verkort weergegeven. 


a' = k.a k = verkleiningsfactor 


De rechte hoek tussen de ribben BA en AD in het grondvlak lijkt 
eveneens kleiner te zijn geworden (zie a in de bovenstaande 
figuur). U kunt gemakkelijk nagaan dat door het vastleggen van 
alpha (bijvoorbeeld 45°) en k (bijvoorbeeld 0,5) de projectierich- 
ting eenduidig bepaald is. 


Wat we nodig hebben zijn transformatievergelijkingen die voor een 
bepaalde a en k de coördinaten x,y,z van een punt op de ruimte- 
lijke kubus projecteren op de coördinaten x' en y' van het gepro- 
jecteerde punt in het platte (projectie-) vlak. Deze formules vor- 
men de kern van alle volgende programma's. 


Afleiding van de transformatievergelijkingen ol 


Afleiding van de transformatievergelijkingen 


Het projectievlak is ons beeldscherm. We gebruiken de MSX- 
schermresolutie van 256 bij 192 punten. De oorsprong van het ruim- 
telijke coördinatensysteem (xyz) moet precies in het midden van het 
beeldscherm geprojecteerd worden. De oorsprong van het beeld- 
schermcoördinatenstelsel ligt weer in de linkerbovenhoek. De x-as 
wijst naar rechts; de y-as naar beneden. 


De afbeelding van het ruimtelijke punt P(x,y,z) op het punt 
P'(x',y') in het platte vlak komt tot stand met behulp van de vol- 
gende transformatievergelijkingen : 


U +x + k.y.cosa èn 
V = (k.y.sina + z) 


x! 
y' 


Hoe we aan deze vergelijkingen komen is uit onderstaande figuur 
direct (nou ja, direct... ) duidelijk. 


y 
P(x,y,z) =p'(x',y) 


| 


Nemen we C = KxCOS(A), S = KxSIN(A) en H = 0,5, dan zijn de 
transformatieformules in BASIC: 


X2 = INT(U+X+C*Y+H) 
Y2 = INT(V-S*Y-Z+H) 


Meer wiskunde hebben we voor dit hoofdstuk niet nodig. 
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Programma 22 laat zien hoe je van elk, door rechte lijnen begrensd 
lichaam (polyeder) een projectie kan maken. Hiervoor is het nood- 
zakelijk dat de waarden van de coördinaten x,y,z in alle hoekpun- 
ten bekend zijn. Als voorbeeld van het gebruik van dit programma 
tekenen we een kubus met een ingeschreven achtvlak (oktaëder). 
De gevolgde programmeertechniek leidt ook bij een viervlak, prisma 
of pyramide tot het gewenste resultaat. 


We nemen aan dat het middelpunt van de kubus samenvalt met de 
oorsprong van het ruimtelijke coördinatensysteem. Om het program- 
ma ‘sneller! te maken nemen we de coördinaten van de acht hoekpun- 
ten ABCDEFGH van de kubus en de coördinaten van de zes hoek- 
punten IJKLMN van de achthoek op in DATAregels. Het tekenen 
van de ribben wordt door de letterstring Z$ bestuurd. Zo betekent 
ZS = "ABBCCDDA" dat de computer van A naar B, van B naar C, 
van C naar D en van D naar A rechte lijnen moet trekken. Met de 
MIDS-functie halen we steeds één letter uit de string Z$ en maken 
er vervolgens met de ASC-functie een getal van (A=1, B=2, C=3, 
EEE ). Hiermee moet de werking van het programma duidelijk zijn. 
Kies voor alpha 45° en voor k 0,5; dat geeft de mooiste projectie. 


100 


110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
s1:0 


320 
530 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
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REM programma 22 kubus 
met achtvlak 
CLS 
INPUT "alfa in graden (45) MieÂ 


INPUT “verkleiningsfactor (.5)";K 
PI=4&xATN(1) 
U=128:V=96:H=.5:RD=EPI/180 
WEA*RD z: CEK#COSCW) 7 S=K*SINCW) 
DIM X(14),Y(14) „ZS (2) 
FOR J=1 TO 14 
READ X,Y,Z 
XCJDZINTCUIXHCHY+H) 
YOJ)=ZINT(V=S*Y=Z+H) 
NEXT J 
CLS 
COLOR 1,15,4 : SCREEN 2 
DEF FNX(X)=INT(37+X/1.4+.5) 
FOR N=1 TO 2 
READ Z$(N) : L=LEN(ZS(N)) 
FOR M=1 TO L=-1 STEP 2 
I=ASC(MIDS(ZS(N) „M,1))=64 
J=ASCC(MIDS(ZE(N) ,M+1,1)) =64 
GENE: CENXERCIIN ZY CID 
CENX CX) DY) 
NEXT M 
NEXT N 
AS=INPUTS$S(1) 
COLOR 15,4,4 : END 
1 


DATA =60,-60,-60, 60,-60,-60 
DATA 60, 60,-60,-60, 60,-60 
DATA =60,-60, 60, 60,-60, 60 
DATA 60, 60, 60,-60, 60, 60 
DATA 0,0,-60,0,-60,0,60,0,0 
DATA 0,60,0,-60,0,0,0,0,60 


DATA "ABBCCDDAAEBFCGDHEFFGGHHE" 
DATA "IJIKILIMJKKLLMMJJNKNLNMN" 
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In programma 23 wordt in een kubus een zadelvlak getekend. Het 
geeft het idee van een gekromd vlak in een kubus. 


Bekijk de bovenstaande figuur. De diagonalen BG en ED zijn elk 
in acht (n) gelijke stukken verdeeld en de zo ontstane punten zijn 
paarsgewijs door een rechte lijn met elkaar verbonden. Kies in het 
programma voor n de waarde 18 of 15. 


“ 


100 


110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 


320 
330 
340 
350 
360 
370 
580 
390 
400 


410 


420 
430 
440 
450 
460 
470 
480 
490 
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REM programma 23 kubus 
met zadelvlak 


CLS 

INPUT "alfa in graden (45) VEA 
INPUT "verkleiningsfactor (.5)"';K 
INPUT "hoeveel lijnen (32) “IeN 


PI=Z4&KATN(I) 

U=128:V=96:H=.S:RD=PI/180 

WEAERD 3 C=K*COSCW) 7 S=KESINCW) 

DIM X(8),Y(8) 

FOR J=1 TO 8 

READ X,Y,Z 
X(JD=INTCU-X+CHY+H) 
YOJD=ZINT(V-S*Y=-Z+H) 

NEXT J 

GES 

COLOR 1,15,4 : SCREEN 2 

DEF FNX(X)=INT(37+X/1.4+.5) 

READ Z$ : L=LEN(ZS$) 

FOR M=1 TO L=1 STEP 2 
I=ASC(MIDS(Z$S,M,1)) -64 
J=ASC(MIDS(ZS,M+1,1)) =64 
LINE KPNWERCIJI EN EIJ De 

CENX CKD D ,YCIID 

NEXT M 

FOR J=0 TO N 
XAZINT(XC2)IN EIJK ZD -XC2)) NH) 
YIZINT(Y (2) + IKCY TD -YC2I ID N+H) 
X2=INT(X (5) AIC (4D XSI DID IN+H) 
Y2zINTC(Y(S5)+IJKCY 4D =YC5I ID IN+H) 
LGENEGENKGERTD) VAINSCENGEENEN „WZI 

NEXT J 

LINE GENXCXG2I) sY C2) Je 

GENK CCL) INCA 

LINE CENKENCHI JY CSI AS 

CENXCX (4D) ,„YC4)) 

AS=INPUTS$(1) 

COLOR 15,4,4 : END 

1 


DATA =60,-60,-60, 60,-60,-60 
DATA 60, 60,-60,-60, 60,-60 
DATA -=60,-60, 60, 60,-60, 60 
DATA 60, 60, 60,-60, 60, 60 
DATA "ABBCCDDAAEBFCGDHEFFGGHHE" 
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Met programma 24 kunnen we, naar keuze, cilinders, kegels of 
afgeknotte kegels tekenen. Als we voor Rl en R2 dezelfde waarde 
invoeren, krijgen we een cilinder. Als R2 kleiner is dan R1 krijgen 
we een afgeknotte kegel en voor R2 = 0 krijgen we een kegel. 


We plaatsen het lichaam zo dat het grondvlak bij z = -60 ligt en het 
bovenvlak (de top) bij z = +60. Het programma tekent zeven 
breedtecirkels met een onderlinge afstand van z = 20 en 16 lijnen 
op de kegel- (of cilinder-) mantel loodrecht op de breedtecirkels. 
Tot slot wordt verticaal de z-as getekend. 


100 REM programma 24 cylinders en 
kegels 
A0 GS 
120 INPUT "alfa in graden (45) EA 
130 INPUT "verkleiningsfactor(.5)'';K 
140 INPUT "stralen RI en R2 "s;RI,R2 
150 PI=4&kATN(1) 
160 U=128:V=96:H=. 5: RD=EPI/180 
170 DR=E(RI=-R2)/6:N=0 
180 WEAERD 7 C=K&COS(W) 7: S=K*SINCW) 
190 CLS 
200: COLOR 115,4 5 SCREEN 2 
210 DER FNXOOSINT (GBTEN/ 1 ebt. 5) 
220 FOR Z==60 TO 60 STEP 20 
230 R=ERI=-N*DR 
240 FOR W=0 TO 360 STEP 3 
250 WI=W&RD 
260 XERKCOS CWI) : YERKSINCWI) 
270 XX=INTCU4X+CHY HH) 
280 YYzZINT(V=S«Y=Z+H) 
290 IF W=0 THEN X1=XX:Y1=YY 
ELSE Xe va=YYe 
LINECENXCXTDI YI 
GENK: CK2I ; VL) 3 


xi=x2:vi=v2 
300 NEXT W 
310 N=N+1 
320 NEXT Z 
330 FOR W=0 TO 360 STEP 23 
340 WI=W*RD 
350 XERIKCOS CWI) 2 YERIKSINCWI) 
360 XAI=INT(U4X+C#Y+H) 
370 YA=ZINT(V=S*Y+60+H) 
380 KER2KCOSCWIJ : YERZESINCWI) 
390 X2=INT(U4X+C*Y+H) 
400 Y2=INT(V=S*Y=60+H) 
410 LINECENX (XI) YI) =CENX(X2) ,Y2) 
4620 NEXT wW 


430 LINE (FNXCU),0)=(FNX(U), 192) 
440 AS=INPUTS(1) 
4650 COLOR 15,4,4 : END 
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Erg leuk is het tekenen van een bol met breedtelijnen en meridianen. 
Zoals bekend is kunnen we elk punt op het boloppervlak eenduidig 
vastleggen met de straal van de bol (r), de geografische breedte 
(@) en de geografische lengte (}) (zie figuur). Het omzetten van 
deze bolcoördinaten r‚p,‚) in cartesische coördinaten (x,y,z) 
geschiedt met de volgende drie vergelijkingen: 


rcospcos 
reospsin) 
rsinp 


Xx 
y 
z 


De bolvergelijking in cartesische vorm is overigens x°+y?+z?=r?. 
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In programma 25 wordt de hoek p door de W en de hoek à door P 
vertegenwoordigd. Kies voor alpha 90° en voor k 0,5. Andere waar- 
den vervormen de bol. Zie de appendix voor mooiere bollen. 


100 REM programma 25 bol 


110: CLS 

120 INPUT "alfa in graden (90) sA 
130 INPUT "verkleiningsfactor(.5)'";K 
140 INPUT "straal,maximaal 75 WER 


150 PI=4kATNCI) 
160 U=128:V=96:H=.5:RD=PI/180 
170 W=A*RD 7 C=KKCOS(W) 7 S=K#SINCW) 
180 CLS 
190 COLOR 1,15,4 1 SCREEN 2 
200 DEF FNX(X)=INT(37+X/1.4+.5) 
210 FOR W==90 TO 90 STEP 15 
220 WI=WERDERIERKCOS (WI) 
230 FOR P=0 TO 360 STEP 3 
240  P1=P4RD 
250  X=1.154RI4COS(PI) :YERIKSINCPI): 
ZER*SINCWI) 
260 _ XX=INT(U+X+SEN(P) *C*Y+H) 
270 _ YY=INT(V-SGN(P)*S*Y-Z+H) 
280 IFP = 0 
THEN X1=XX:Y1=YY 
ELSE X2=XX:Y2=YY: 
LINECENXCX1) „Y1)=- 
CENXCX2) „Y2): 
xi=X2:Yi=Y2 
290 NEXT P 
300 NEXT W 
310 FOR P=0 TO 180 STEP 15 
320 P1=P&RD 


rd TERREN 
ESD 
GATEN 


330 FOR W=0 TO 360 STEP 3 Pons EKD 
340  WI=WERDERTERKCOS (WI) (enen 
350  X=1.154RIKCOS(PI) Veenen) 
360  Y=RIKSIN(P1) & CUD, / 


HN 1) Jb 


370  Z=RASINCWI) Nr 


380 _ XX=INT(U4X+C*Y+H) 
390 _ YY=INT(V=S*Y-Z+H) 
400 IF W=O THEN X1=XX:Y1=YY 
ELSE X2=XX:Y2=YY: 
LINECENX (XI) YI) = 
CENX(X2) ,Y2): 
Xi=X2:Y1=Y2 


410 NEXT W 

420 NEXT P 

430 AS$S=INPUT$S(1) 

440 COLOR 15,4,4 : END 
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Een geliefd demonstratieprogramma van computerleveranciers is een 
om zijn as draaiend driedimensionaal lichaam. Meestal wordt als 
'draai-as' de z-as gekozen, die dan samenvalt met de lengte-as van 
het lichaam. 


Microcomputers die alleen een BASIC-vertolker bezitten zijn te lang- 
zaam om zo'n draaiing real-time uit te voeren. Er zijn namelijk nogal 
ingewikkelde trigonometrische berekeningen voor nodig, die tame- 
lijk veel tijd in beslag nemen. Daarnaast duurt het tekenen van het 
lichaam in een bepaalde stand in BASIC zo'n 1 à 2 seconden. Zouden 
we in een BASIC-programma steeds het lichaam tekenen, uitwissen, 
draaien, tekenen, uitwissen, draaien, ... enz., dan zien we het 
draaiende lichaam in plaats van een vloeiende beweging een schok- 
kerige beweging maken. De beste oplossing voor dit probleem is de 
volgende. 


Deel de totale middelpuntshoek van 360° in n even grote hoeken. 
Voor elk van deze hoeken 


ä. = 360° pn 
1 nr Se 

berekenen we van te voren de coördinaten van de hoekpunten van 
het lichaam. Voor elke stand slaan we deze hoekpuntscoördinaten 
in een array op. Al deze berekeningen voeren we in BASIC uit bij 
een leeg beeldscherm. Als de berekeningen klaar zijn wordt een 
machinetaalprogramma uitgevoerd dat de eerste groep coördinaten 
uit de array haalt, het lichaam tekent, na een bepaalde tijd het 
beeldscherm wist, de volgende groep coördinaten ophaalt, .………. 
enz. Deze oplossing heeft een redelijk vloeiende draaiing tot gevolg. 
Dit programma is geschreven op een CBM 3016 die een 6502-micro- 
processor bezit. Hier doen we echter anders (zie volgende pagina). 


ww =n.w, = 3609 
n 1 
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Programma 26 laat een driezijdig prisma om de z-as draaien. De z- 
as is zowel de lengte-as als de symmetrie-as van het prisma. Het 
programma tekent een prisma en wacht vervolgens met het draaien 
en opnieuw tekenen tot we een toets indrukken. Elk volgend prisma 
wordt over zijn voorgangers heen getekend. We kunnen het draaien 
dus vertraagd gadeslaan. 


100 REM programma 26 draaiend prisma 
118-EES 

120 INPUT "alfa in graden (45) PIAN 
130 INPUT "verkleiningsfactor(.5)';K 
140 INPUT "draaihoek in gr. (45)'"';OM 
150 PI=4KATN(1) 

160 U=128:V=96:H=. 5:RD=PI/180:R=80 
170 WEAERD : CEK*COSCW) 7 S=K&SINCW) 
180: CES 

190 COLOR 1,15,4 s SCREEN 2 

200 DEF FNX(X)=INT(37+4+X/1.4+.5) 

210 LINE (FNXCU),0)=(FNX(U),„192) 

220 DIM X(7),Y(7) 

230 FOR W=0 TO 360 STEP OM 

240 FOR J=0 TO 3 

250 WI=CWEJKT20) *RD 

260 XERKCOS CWI) 

270 YERKSINCWI) 

280 XCJDZINTCUFXHCHY HH) 

290 YOJDZINT(V=-S«Y+604H) 

300 X(JT4) EX (JI) 

310 YAJ+4)D=ZINT(V=S«Y=-60+H) 

520 NEXT 


350 "* 

340 FOR J=0 TO 2 

350 LINE CENXCXCJII,YCII D= 
CENX(XCJHIDD ,YCIHID) 

360 LINE CFNXCXCJ)I,YCIDD= 
CENXCXCJH4DD,YJH4)) 

370 LINE (FNX(CX(J+4)) ,Y(IJH4DD= 


CENX(XCJ+5)I DI ,YCI+H5)) 
380 NEXT J 
390 AS$S=INPUTS(1) 
400 NEXT W 
410 COLOR 15,4,4 : END 
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Programma 27 tekent de projectie van een regelmatig twintigvlak 
(ikosaëder). Het is bekend dat er slechts vijf 'echt-regelmatige' 
veelvlakken (polyeders) zijn, namelijk een tetraëder (regelmatig 
drievlak), een kubus (regelmatig zesvlak), een octaëder (regelmatig 
achtvlak), een dodekaëder (regelmatig twaalfvlak) en een ikosaëder 
(regelmatig twintigvlak). Een ikosaëder heeft 12 hoeken, 30 ribben 
en 20 vlakken. Een vlak is een gelijkzijdige driehoek. Als we pro- 
gramma 22 willen gebruiken om zo'n ikosaëder te tekenen dan heb- 
ben we de coördinaten van alle 12 hoekpunten nodig. Als we het 
ikosaëder in een speciale stand zetten kunnen we de hoekpuntcoör- 
dinaten met behulp van de techniek van de 'gulden snede! relatief 
eenvoudig berekenen. Wie zich voor de hierachterliggende wiskun- 
dige theorie interesseert moet er maar eens een meetkundeboek op 
naslaan. De 'gulden-snede'-techniek deelt een lijnstuk AB in twee 
stukken AT en TB op zo'n manier dat de volgende evenredigheid 
geldt: 


AB. AR AT £ TB 
el 
EN sb B 


Kiezen we voor AB de lengte 1, dan kunnen we de lengte van AT 
berekenen; immers 


dus dan moet gelden t? = 1-t ofwel t?+t-1=0. 
Met de alombekende abe-formule berekenen we nu 


t ml 1 +tvVl*+4 
En tn 
De positieve wortel geeft: toe À 5 vs 


Dit is de lengte van het grootste van de twee stukken die we krij- 
gen als we een lijnstuk, ter lengte 1, volgens de 'gulden snede! in 
twee stukken verdelen. (t = 0,618; 1-t = 0,382). 

Deze waarde wordt in regel 170 berekend en in de regels 220-339 
gebruikt om de hoekpuntcoördinaten van het ikosaëder te bereke- 
nen. Om te zorgen dat het twintigvlak voldoende groot getekend 
wordt zijn alle coördinaten met 80 (F) vermenigvuldigd. Kies voor 
alpha (A) 90° en voor k (K) de waarde 0,4. Andere waarden ver- 
tekenen het beeld. Alpha = 180° en K=0,6 geeft trouwens iets 
onverwachts! 
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100 REM programma 27 ikosaeder 

119 ELS 

120 INPUT "alfa in graden (90) MEN 
130 INPUT "verkleiningsfactor(.4)";K 
140 PI=Z4XATN(I) 

150 U=128:V=96:H=. 5:RD=PI/180 

160 WEAERD z£ CEK#COSCW) 7 S=KESINCW) 
170 T=(SQAR(5)=-1)/2 : F=80 

180 CLS 

190 COLOR 1,15,4 : SCREEN 2 

200 DEF FNX(X)=INT(37+X/1.4+.5) 

210 DIM XET2IE CTI EZ CT LIZE CS) 


220 x(1)= Oe NC FAT ef 
230 X(2)= Os YC2)= =FeT: ZIE sf 
240 X(3)= FP: YC5)= 0: 2Z(3)==FAT 
250 X(4)= =Fs Y(4)= 0: Z4)=-F4T 
260 XC5)= FATs YC5)= F: Z(5)= 0 
270 X(6)=-FT: YC6)= F: Z(6)= 0 
ar MOTIE FRTE YOT =p ZE 0 
280 X(8)=-FAT: Y(B)= Rs 23E 0 
290 X(9)= F: Y(9)= Ös 209) BET 
300 X(10)= =F: Y(10)= 0:Z(10)= F*T 
510 XC11)= Oz YLTTI= FAT:ZC I= F 
se0 KCI2I= Os VAI PREZ NS F 


330 FOR N=1 TO 3 

340 READ ZS(N) : L=LEN(ZS(N)) 
350 FOR M=1 TO L-1 STEP 2 

360 _ I=ASC(MIDS(ZS(N),M,1))=64 
370 _ J=ASC(MIDS(ZS(N) ,M+1,1))-64 
380 _ X1=INTCU4X(IJ+CHY(I)+H) 

390 _ YAZINT(V-S*Y(I)=Z(I)+H) 

400 _ X2=INTCUFX(II+CHY (IJ) +H) 

410  Y2=INT(V=S*Y(I) =D +H) 

420 GENE (CENWCRTD TI CENK CNS) YZ) 
430 NEXT M 

440 NEXT N 

450 AS=INPUTS(1) 

460 COLOR 15,4,4 : END 

4670 * 

480 DATA "BCCEEFFDDBABACAEAFAD" 
490 _ DATA "BHHDDJJFFKKEEIICCGGB" 
500 DATA "GHHJJKKIIGLGLHLJLKLI" 
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ruimte 


In het vorige hoofdstuk hebben we ons uitvoerig beziggehouden met 
het tekenen van driedimensionale lichamen. We herhalen hiervan nog 
eens de belangrijkste punten: 


1. We gebruiken de parallelprojectie om een driedimensionaal 
lichaam met n hoekpunten P(x,y,z) in een plat vlak te tekenen. 
Voor de projectiehoek alpha (a) gebruiken we bij voorkeur de 
waarde 45° of 60°. De schuin-naar-achteren lopende ribben 
worden korter getekend dan ze in werkelijkheid zijn. Als ver- 
kleiningsfactor kiezen we vaak 1/2 of 1/3. 


2. Elk punt P(x,y‚z) van het ruimtelijke lichaam wordt op een 
punt P'(x',y') in het projectievlak geprojecteerd. De hierbij 
gebruikte projectieformules zijn : 


x'=U + Xx + k.y-cosa 
y'= V = k.y.sina = z 


Zoals gebruikelijk zijn U en V de coördinaten van het midden 
van het beeldscherm. Tot nu toe hebben we U = 128 en V = 96 
verondersteld. 


In de volgende BASIC-programma's zien we de bovenstaande trans- 
formatievergelijkingen in de vorm: 


XG 
YG 


INT(U+XX+C*YY+H) 
INT(V-SxYY-Z+H) 


De betekenis van de gebruikte variabelen is: 


XG,YG : beeldschermecoördinaten x'‚y' 
RNL Y : de lopende coördinaten x,y 
® : KxCOS(We&RD) 

S : K*SIN(W*RD) 


w : projectiehoek in graden 
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RD : factor 7/180 voor omrekenen van graden in radialen 
K : verkleiningsfactor 


Tekenen van driedimensionale functies 


Als paradepaardje van menige demonstratie van Hoge-Resolutie- 
Graphics wordt vaak een programma gebruikt dat een of andere 
fraaie grafiek van een driedimensionale functie tekent. Als leek 
vraag je je af hoe ze weten welke functies ze moeten nemen, hoe het 
tekenen van de grafiek van zo'n functie geprogrammeerd wordt, en 
hoe ze het programmatechnisch voor elkaar krijgen dat de 'onzicht- 
bare lijnen' ook inderdaad niet getekend worden. Bij dergelijke 
demonstraties wordt doorgaans geen programmalisting getoond. 
Mocht dit wel het geval zijn dan is het programma vaak zo machine- 
afhankelijk dat het omzetten van het programma naar een andere 
machine lastig, zo niet ondoenlijk is. 


In dit hoofdstuk hopen we u antwoord te kunnen geven op de vol- 
gende vragen: 


: @ Hoe vind ik 'mooie' driedimensionale functies? 


2. Hoe ziet een algemeen programma voor het tekenen van een 
dergelijke functie eruit? 


3. Hoe onderdruk je het tekenen van de onzichtbare lijnen 
(hidden lines)? 


We beginnen met de definitie van een driedimensionale functie. 


Elke functie van de vorm z = f(x,y) heet een driedimensionale func- 
tie en vormt een, meestal, gekromd vlak in een driedimensionaal 
(ruimtelijk) coördinatenstelsel. Speciaal voor niet-wiskundigen volgt 
nu een voorbeeld van een functie z = f(x,y) en het gekromde vlak 
dat als 'grafiek' bij de functie hoort. 


Stel dat we voor z = f(x,y) de volgende vergelijking nemen: 


-(x?+v? 
=e Kragt) (e is het grondtal van de 
natuurlijke logaritme; 


e = 2,718284183...) 


z 


Voor elk punt (x,y) in het (platte) X-Y-vlak kunnen we nu de 
daarbij horende waarde van z uitrekenen. Als x = 0,5 en y = 0,1 
dan is 
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- & (0,25*0,01) 
-0,26 
e 


knd 


ZZ = nd 


z = 0,771, want 2,71828... tot de macht -0,26 is 0,771. 


Zet nu (in gedachten) in het voetpunt P(0,5;0,1) in het X-Y-vlak 
een staafje met een lengte van 0,771 rechtop. Doe ditzelfde voor 
nog een groot aantal punten (x,y). 


Leg nu over al deze staafjes een elastische folie. Dit (golvende) 

stuk folie vormt een goed model van het vlak met vergelijking 
Ee 

z=e ee R, 

Hoe dit eruit zou zien kunt u zien op p.73. 
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Zoek met behulp van programma 7 uit hoofdstuk 2 een willekeurige 
econtine functie die symmetrisch ten opzichte van de y-as is. De 
grafiek van de functie moet 'bergen en dalen! vertonen (de functie 
moet maxima en minima hebben). Erg geschikt zijn trigonometrische 
functies (sin, cos) en combinaties van deze functies. Ook geschikt 
zijn functies waarin alleen termen voorkomen met 'even-machten' 
van x en exponentiële functies. 


Enkele voorbeelden : 


2 
y=e 8 4 y = sin(x) 
n _ COsS(3X)  cos(5X) cOS(TX) 
NENDE Tp Ma 
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Bekijk nu de grafiek van zo'n functie in het interval -asx<a. Dit 
‘stuk grafiek gaan we draaien rond de y-as. Zo ontstaat een, ten 
opzichte van de y-as, draaisymmetrisch driedimensionaal vlak. Als 
we nu de y-as als z-as kiezen, dan wordt de vergelijking van een 
dergelijk draaisymmetrisch ruimtelijk vlak (rond de z-as) van de 
vorm 


z=f(r) met r = Vx?+y? 


We zullen dit aan een eenvoudig voorbeeld proberen te verklaren. 
We kiezen als voorbeeld de eenvoudige paraboolvergelijking y=x?. 


Stel we bekijken de grafiek voor -2sxs2: 


Als we nu de x-as rond de y-as laten draaien (de x-as komt als het 
ware loodrecht het papier uit), dan ontstaat een draaisymmetrisch 
vlak rond de y-as. Dit is een kegel: 
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Elk punt P op de kegelmantel heeft de eigenschap y=r?, waarbij r 
de afstand is van het geprojecteerde punt P' tot het punt O. We 
kunnen het vlak V dat door de ronddraaiende x-as wordt gevormd 
als volgt tekenen: 


Welnu, als we nu de z-as als y-as nemen en we nemen de y-as in 
vlak V loodrecht op de x-as, dan zien we dat de afstand r geschre- 
ven kan worden als 


r = VX?+y? (stelling van Pythagoras, zie rechts in de 
volgende tekening) 
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en nu wordt y=r? geschreven als z=r? met r? = x?+y?, dus de ver- 
gelijking van de paraboloïde (kegel) wordt dan 


Z= x?+y? 


Dus y=x? wordt bij draaiing om de y-as en bij een z-as die de plaats 
van de y-as inneemt z = x?+y?. 
Voor de bovengenoemde functies geldt het volgende: 


seer wordt PE a di 
y= wordt z = mr) 
y = cos(x) - cos(3x) ,cos(5x) _ cos(7Xx) ora 


met r = Vx?+y? 


Voor deze klasse van draaisymmetrische figuren ontwikkelen we een 
algemeen tekenprogramma. 
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Antwoord op de tweede vraag 


Bekijk de onderstaande tekening. U kijkt recht bovenop het ruimte- 
lijke vlak van een driedimensionale functie. We willen dat het op het 
grondvlak geprojecteerde vlak een vierkant is met zijden van 2a, 
Straks zullen we zien dat, als we een mooie tekening van zo'n ruim- 
telijk vlak op het beeldscherm willen maken, a=95 een heel goede 
keus is. 


Als ‘doorgewinterde! programmeur ziet u natuurlijk direct dat het 
programma voor het tekenen van het vlak een geneste FOR-NEXT 
(dx‚dy) zal bevatten! Als we een vlak als grafiek van z = f(x,y) 
willen tekenen, moeten we namelijk voor een aantal combinaties van 
x en y de functiewaarde z = f(x,y) berekenen. We beginnen met 
y=-95 en laten x met stapjes dx (bijvoorbeeld dx=3) van -95 naar 
+95 lopen. Voor elke combinatie x,y met y=-95 berekenen we de 
functiewaarde z = f(x,y). Zo ontstaan de punten (x,y,z) van het 
ruimtelijke vlak waarvan de geprojecteerde punten (x,y) in de 
bovenstaande tekening op de onderste zijde van het vierkant als 
bolletjes getekend zijn. Nu verhogen we y met dy en laten x weer 
van -95 tot +95 lopen. Dit geeft de tweede reep van het vlak. Op 
deze wijze ontstaat het hele vlak, waarvan de punten (x,y,z) 
natuurlijk nog getransformeerd moeten worden naar beeldscherm- 
coördinaten (x',y'). 


Een eerste globale opzet voor het tekenen van een 'vierkant'-stuk 
vlak is: 
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FOR Y=-95 TO 95 STEP DY 
FOR X=-95 TO 95 STEP DX 
GOSUB 1000 
NEXT X 
NEXT y 


In de subroutine 1000 wordt voor een punt (x,y) de waarde van 

z = f(x,y) berekend. Met de bekende transformaties wordt vervol- 
gens het punt P(x,y,‚z) overgebracht naar een punt P'(x',y') op 
het beeldscherm. Nu kan het berekende punt P(x,y,z) dat in het 
vlak ligt, als het punt P'(x',y') op het scherm worden getekend. 


Jammer genoeg lenen niet alle driedimensionale functies zich voor de 
intervallen -95s xs 95 en -95sys 95; vandaar dat de waarde voor a 

met een INPUT-opdracht ingelezen wordt. Als we in het programma 

toch -95s XX5s 95 en -I5s YY= 95 kiezen, moet de ingetoetste waarde 

A als volgt gebruikt worden: 


= \ En A 
X = XX. 55 en Y = YY.55 
dat wil zeggen er geldt: -AsXsA en -AsYsA 


De z-waarden van met name de trigonometrische driedimensionale 
functies zal tussen -1 en 1 liggen. Om deze waarden wat 'op te bla- 
zen' kan een vermenigvuldigingsfactor (Kl) worden ingetoetst. 
Goede waarden voor Kl liggen tussen 30 en 80. 


We kunnen ons programma nu als volgt opschrijven: 


opzet voor een programma voor 
het tekenen van de functie 
Z=f XY) 


variabelen waarde geven en 
hoge resolutie instellen 
FOR YY==95 TO 95 STEP DY 
Y=YYxaAf 
FOR XX==95 TO 95 STEP DX 
XEXXAF 
GOSUB 1000 ‘functiewaarde 2 
XG=INT(U4XX+C+HYY+H) 
YG=INT(V=S*YY=Z +H) 


teken punt P(XG,YG) 


NEXT XX 
NEXT YY 


programma afsluiten 


REMbeschrijving van de functie 
Zef lX,Y) 
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Zoals beloofd laten we nu zien waarom de intervallen -95< XX < 95 

en -95sYYs95 zo geschikt zijn om driedimensionale functies te teke- 
nen met een hoog oplossend vermogen. We willen graag de fraaie 
projectie met a=45° en k=0,5 gebruiken. Wat worden dan de beeld- 
schermcoördinaten XG en YG voor de hoekpunten 'linksonder' en 
‘rechtsboven! van het vierkant op p.69? Deze punten bepalen immers 
of de hele grafiek op ons beeldscherm van 256 bij 192 puntjes gete- 
kend kan worden. Voor Kl nemen we 50, dat ligt zo'n beetje tussen 
30 en 80! Nemen we voor de z-waarde ook 50, dan gaat het punt 
(x,y,z) met coördinaten (-95,-95,50) over in het beeldschermpunt 
(XG,YG) met 


XG 
YG 


INT (128-95-0,5.95.cos(45)+0,5) > 
INT(96+0,5.95.sin(45)-50+0,5) = 


XG 
YG 


== 
80 


uu 
uu 
un 


De YG-waarde ligt tussen 0 en 192, maar de XG-waarde ligt net bui- 
ten het scherm. Voor dit hoekpunt is dit echter geen bezwaar, het 
ligt immers in de uiterste hoek van onze tekening! De 'breedte' (256) 
van het scherm wordt heel goed benut, kijk maar eens naar de coör- 
dinaten (XG,YG) voor het punt (95,95,50) rechtsboven: 


XG = INT(128+95+0,5. 95.cos(45)+0,5) » 257 
YG = INT(96-0,5.95.sin(45)-50+0,5) = 12 


We benutten dus ruimschoots de hele schermbreedte (-1sXG5< 257) 
voor het tekenen van de driedimensionale figuur. In programma 28 
(regel 300) en in programma 29 (regels 330 en 340) is rekening 
gehouden met eventuele negatieve XG-waarden en XG-waarden gro- 
ter dan 255. Als dit voorkomt, worden deze waarden op respectieve- 
lijk O en 255 gezet, zodat het programma niet door een foutmelding 
afbreekt. 


Op de onder- en bovenrand (YY=-95 en YY=+95) gelden vaak heel 
kleine z-waarden, zodat bij a=45°, k=0,5 en z=0 de grafiek op het 
scherm zal liggen tussen ‚z-waar de 
YG = INT(96+0,5.95.sin(45)-0+0,5 > 130 


YG = INT(96-0,5.95.sin(45)-0+0,5 » 62 


Ons tekenvlak is dus ongeveer -1sXG5s257 èn 
62s YG< 130 


In het bovenstaande programmavoorstel wordt gebruik gemaakt van 
puntgraphies'. Hierbij worden de punten puntje voor puntje gete- 
kend. Om een enigszins acceptabele tekening te krijgen moeten erg 
veel 'puntjes' berekend worden (DX=1,DY=1). Dit duurt in BASIC 
(geneste lussen en vele trigonometrische berekeningen) erg lang. 
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Voor de onderstaande tekeningen betekent dit al snel een tekentijd 
van meer dan 30 minuten, soms wel een uur. Dit is weinig bevredi- 
gend! Veel sneller gaat het als we de 'vectorgraphics'-methode 
gebruiken. Deze techniek hebben we in alle voorgaande programma's 
gebruikt; we verbinden steeds twee naast elkaar liggende punten 
door een recht lijntje. Hierbij hoeven we lang niet zoveel punten te 
berekenen als bij de 'puntgraphics'-techniek. Bovendien ziet de 
tekening er beter uit. 


Als antwoord op de tweede vraag komen we dan met het volgende 
algemene programma voor het tekenen van draaisymmetrische ruim- 
telijke vlakken. 


100 REM programma 28 grafiek van de 
functie Z=f(X,Y) 
110 CLS 
120 INPUT „alfa var grad: 45155): Ket 
130 INPUT -"'verklesfactorts5.75) — "sk 
140 INPUT "rechtergrens x (>0) a \ 
150 INPUT "vergr.factor functie Z'"';K1 
160 PI=4XATN(1) 
170 U=128:V=96:H=.S:RD=PI/180 
180 C=K*COSCWERD) z: S=KKSINCWERD) 
190 DX=3 : DY=5 £ AF=A/95 
200 CLS 
€10 GOLOR 41554 ft SCREEN 2 
220 FOR YY==95 TO 95 STEP DY 
230 Y=YYxAF 
240 FOR XX==95 TO 95 STEP DX 
250 XEXXkAF 
260 GOSUB 1000 'functiewaarde Z 
270 XG=INTCUIXX+HCHYY+H) 
280 YG=INT(V=S*YY=Z +H) 
290 IF XG<O THEN XG=0 
300 IF XG>255 THEN XG=255 
310 IF YG<O OR YG>191 
THEN PRINT "foute K1": STOP 
320 IF XX==95 
THEN X1=XG:Y1=YG 
ELSE X2=XG:Y2=YG: 
GENECKT SY 1 ICR SY AI 8 
xi=X2:Y1=Y2 
330 NEXT XX 
340 NEXT YY 
350 AS$S=INPUT$(1) 
360 COLOR 15,4,4 : END 
st 
1000 REM subroutine functiewaarde Z 
1010 ZEKTKEXP(- XXY) 
1020 RETURN 
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Met dit programma maakt u de onderstaande tekening. Kies op uw 
MSX-computer de volgende waarden: W=45°®, K=0,5, A=3 en K1=70. 


— (2 +v2 
zn en) 


U zult over de bovenstaande tekening niet tevreden zijn. De voor- 
grond is goed, maar in het achterste stuk zijn ook alle lijnen die 
voor ons onzichtbaar zijn getrokken. In onderstaande tekening zien 
we het resultaat met dezelfde waarden voor W‚, K,‚, A en Kl, maar 
waarin de niet-zichtbare lijnen ook niet getekend zijn. Hoe krijgen 
we dit voor elkaar? 
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Bekijk de onderstaande tekening. Een punt P(x,y,z) van het te 
tekenen vlak is alleen dan zichtbaar wanneer: 


a. het punt 'hoger' ligt dan alle voorgaande punten met dezelfde 
x-coördinaat of 

b. het punt 'lager' ligt dan alle voorgaande punten met dezelfde 
x-coördinaat. 


Vanzelfsprekend geldt dit ook voor de op het beeldscherm gepro- 
jecteerde punten P(XG,YG). 
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Bedenk bij de navolgende uiteenzetting dat de coördinaten (XG,YG) 
van een beeldschermpunt 'gemeten' worden ten opzichte van de 
linkerbovenhoek van het beeldscherm, de oorsprong van ons HRG- 
coördinatenstelsel. 


Een punt, in beeldschermcoördinaten uitgedrukt, ligt dus 'hoger!' 
dan eenander punt (bij dezelfde XG-waarde) als de YG-waarde 
groter is en ligt 'lager' als de YG-waarde kleiner is. 


We gaan als volgt te werk. We gebruiken twee arrays (HI(XG) en 
H2(XG)) waarin we voor alle XG-coördinaten in ons vlak de respec- 
tievelijk kleinste en grootste YG-waarde bewaren. Vinden we voor 
een punt met een bepaalde XG-waarde een bijbehorende YG-waarde 
die òf kleiner is dan HI(XG) òf groter is dan H2(XG), dan is het 
punt zichtbaar en wordt HI(XG) of H2(XG) gelijk gemaakt aan deze 
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nieuwe kleinste of grootste YG-waarde. Geldt voor een bepaalde XG 
dat HI(XG)<YG<H2(XG), dan is het punt (XG,YG) niet zichtbaar. 
HI(XG) en H2(XG) zijn te zien als twee horizonnen. 


Als voor een bepaald punt Pi(x,‚y,‚z) uit het vlak het beeldpunt 
PI(XG,YG) berekend is, gaan we kijken of YG kleiner dan of gelijk 
aan HI(XG) is (YG<=HI(XG)). Is dit zo dan is het punt Pl op het 
beeldscherm zichtbaar. Nu zetten we de vlag Fl op 1 en we maken 
HI(XG) gelijk aan YG (HI(XG)=YG). Is dit niet het geval (YG> 
HI(XG)), dan is Pl onzichtbaar en blijft de vlag Fl op 0 staan. 


Nu gaan we naar het volgende punt P:(x,y,z) uit het vlak (we ver- 
hogen x met dx). Opnieuw berekenen we het beeldpunt P2(XG,YG). 
Weer wordt gekeken of YG<=HI(XG) en zonodig wordt de vlag F2 
op 1 gezet. De punten Pl en P2 worden alleen dan door een rechte 
lijn (lijntje) verbonden als beide vlaggen F1 en F2 de waarde 1 heb- 
ben, dus als FlxF2=1. 


Ditzelfde doen we voor de tweede ‘horizon! H2(XG). We kijken of 
LGSEHACKG) sn vase ‚ enzovoorts. 


Omdat we bij een vaste y-waarde de x-waarde steeds met dx opho- 
gen (bijvoorbeeld dx=3) en hiermee steeds twee buurpunten 
PI(XG,YG) en P2(XG,YG) berekenen slaan we als het ware een aan- 
tal punten over waarvoor geen HI(XG)- en H2(XG)-waarden bere- 
kend worden. Door lineaire interpolatie tussen de horizonwaarden 
(HI(XG) en H2(XG)) van twee buurpunten zouden we de horizon- 
waarden voor de hiertussenliggende punten kunnen berekenen. We 
hebben deze waarden wel nodig, omdat het kan voorkomen dat, als 
we y met dy verhogen, we XG-waarden vinden waarvoor nog geen 
HI(XG)- en H2(XG)-waarden berekend zijn. 


Een programma dat de hierboven geschetste werkwijze zou volgen 
zou, in BASIC geschreven, veel te langzaam zijn. Om de tekensnel- 
heid acceptabel te houden, dat wil zeggen niet langer dan 10 minu- 
ten tekenen voor één figuur, passen we de volgende vereenvoudi- 
gingen toe: 


Ls We berekenen alleen de onderste (voor de kijker de bovenste) 
horizon HI(XG). 


Ld 


We passen geen lineaire interpolatie toe bij het berekenen van 
de array-waarden in HI(XG). Alle beeldpunten in een interval 
ter lengte dx krijgen dezelfde H1(XG)-waarde. 


Als we de stapgrootte dx maar klein genoeg kiezen (een goede waar- 
de is dx=3), krijgen we ondanks deze twee simplificaties toch accep- 
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tabele tekeningen. (Bekijk de volgende tekeningen en oordeelt u 
zelf.) Het niet-tekenen van 'voor de kijker! onzichtbare punten 
komt tot stand door het volgende stukje programma : 


F1=0:L=INT(XG/DX) 
IF YG<=H(L) THEN Fl=1:H(L)=YG 


Laten we dit met een voorbeeld illustreren. Stel we berekenen van 
een punt P(x,y,‚z) in het ruimtelijke vlak de projectie op het beeld- 
scherm. Dit gebeurt door eerst met x en y de z-waarde te bereke- 
nen. (In programma 29 krijgt Y in regel 270, X in regel 290 en Z in 
regel 1010 een waarde.) Als de x-, y- en z-waarde van een punt in 
het vlak berekend zijn, wordt dit punt op het beeldscherm gepro- 
jecteerd door het berekenen van de twee coördinaten XG en YG. 

(In programma 29 gebeurt dit in regel 310-320.) Stel nu dat als 
beeldpunt het punt (XG=40 en YG=126) berekend is. Het programma 
zal de vlag Fl op nul zetten (regel 2010) en de waarde L=INT(XG/DX) 
berekenen. Voor XG=40 en DX=3 wordt deze waarde INT(40/3), dus 
L=13. Dit betekent in feite dat de drie punten met XG=39, XG=40 en 
XG=41 dezelfde horizonwaarde L=13 opleveren. Mocht nu eerder in 
H(13) als horizonwaarde de waarde 132 opgeslagen zijn, dan is 
YG<=H(L) (regel 2020) waar, want 126<=132 is waar, en is het punt 
(40,126) zichtbaar. Op dit moment wordt de vlag F1 gezet en wordt 
YG=126 de nieuwe horizonwaarde bij L=13 (regel 2020). 


De regels 


200 DIM H(255) 

210 FOR L=0 TO 255 
220 H(L)=1000 
230 NEXT L 


zorgen ervoor dat, voordat het tekenen begint, het hele scherm 
‘zichtbaar! is. Dit is zo als de onderste horizon (voor de kijker de 
bovenste horizon) de onderrand van het beeldscherm is. De waarde 
1000 geeft in feite een horizon die 'ver onder' het beeldscherm ligt. 
Op p.77 staat het programma dat alle voor ons onzichtbare punten 
ook niet tekent. 


Met W=45°, K=0,5, A=3 en K1=80 krijgt u de grafiek van 


= 2 2 
z=e GRT zoals die een paar pagina's eerder is getekend; de niet- 
zichtbare punten die achter de hoed liggen zijn inderdaad niet te 
zien! 
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REM programma 29 grafiek van de 
functie Z=f(X,Y), hidden Lines 
CES 
INPUT ‘alfa An grade 665135) “el 
ENPUT "verkletactortso=/5) ek 
INPUT "rechtergrens x (>0) er 
INPUT “vergr.tactor- fUAcE Te ZK 
PI=4XkATN(1) 
U=128:V=96:H=.5:RD=PI/180 
CEKKCOSCWERD) 7 S=KKSINCWERD) 
DX=3 : DY=5 z: AF=A/95 
DIM H(255) 
FOR L=0 TO 255 
HCL)=1000 
NEXT L 
CLS 
COLOR Tetaskh. © SCHEEN 2 
FOR Vya=95 TO 95 STEP DY 
Y=YYxAF 5 
FOR XX==95 TO 95 STEP DX 
XEXXAF 
GOSUB 1000 'functiewaarde Z 
XG=INTCUÊXX+HCHYY+H) 
YG=INT(V=S*YY=Z +H) 
IF XG<O THEN XG=0 
IF XG>255 THEN XG=255 
IF YG<O OR YG>191 
THEN PRINT "foute K1'": STOP 
IF XX==95 
THEN GOSUB 2000 
ELSE GOSUB 35000 
NEXT XX 
NEXT YY 
AS=INPUT$(1) 
COLOR 15,4,4 z END 
1 
REM subroutine functiewaarde Z 
ZEKTKEXP (XXY) 
RETURN 
REM subroutine eerste xx=waarde 
F1=0 : L=INT(XG/DX) 
IF YG<=H(L) THEN F1=1:H(L)=YG 
X1=XG:Y1=YG 
RETURN 
REM subroutine volgende xx'en 
F2=0 : L=INT(XG/DX) 
IF YG<=H(L) THEN F2=1:H(L)=YG 
X2=XG:Y2=YG 
IF F1*F2=1 
THEN CUCINE CX1,Y1I CRZ, V2) 
xi=X2:Y1=Y2 
RETURN 
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Hoe kleiner we DX en DY maken, des te gedetailleerder zal de figuur 
worden. Maar ook 'des te langer zal het tekenen duren'. De combi- 
natie DX=3 en DY=5 is een goed compromis tussen de tekensnelheid 
en de mate van gedetailleerdheid. 


Dit programma is een algemeen programma voor het tekenen van 
draaisymmetrische ruimtelijke vlakken. Als u andere figuren wilt 
maken, verander dan alleen de functie in de subroutine 1000 en kies 
mogelijk andere waarden voor de variabelen. 


Nu volgt een aantal voorbeelden. 


1000 REM subroutine functiewaarde Z 

1010 RESQR(XAXAYAY) ARD 

1020 Z=KIK(COS(RI-COS(FER)/3+ 
COS(S*R)/5-COS(7ER)/7) 

1030 RETURN 


U kiest: W=45, K=0,5, A=180°, K1=45. 


5 _ cos(3r) 
ECO) er + 
„eos(5r) _ cos(7r) 

5 ú 


1000 REM subroutine functiewaarde Z 
1010 RESQR(KEKHYAYD ERD 
1020 IF R=0 THEN Z=K1 

ELSE Z=K1*SINCR)/R 


1030 RETURN 


U kiest: W=45, K=0,5, A=1080, K1=60. 
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1000 REM subroutine functiewaarde Z 
1010 R=SQR(XAXHYAY) 

1020 Z=KIKEXP(-COS(R/16)) 

1030 RETURN 


U kiest: W=45, K=0,5, A=80, K1=30. 


zz exp(-cos(75)) 


Voor de volgende twee tekeningen geldt: 


1000 REM subroutine functiewaarde Z 
1010 R=ESQR(KAXHYAY) 

1020 Z=KICOSCR/16) EN ZEKI*SIN(R/16) 
1030 RETURN 


U kiest: W=45, K=0,5, A=90, K1=50. 


Met dit programma hebben we talrijke andere grafieken gemaakt. 
Elke keer dat u een ‘leuke functie' tegenkomt, kunt u meteen de 
daarbij horende draaisymmetrische figuur maken. 
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Tot slot volgt nog een programma voor het tekenen van een fraaie 
vier-toppige grafiek. U bent hem misschien wel eens in een compu- 
tertijdschrift tegengekomen. Om de snelheid op te voeren is de sub- 
routine 1000 in het hoofdprogramma opgenomen. U kunt dit natuur- 
lijk ook in alle andere programma's doen. 


100 REM programma 30 mooie functie 
110 GES 
120 INPUT "alfa in grad.(45=135) "sw 
130 INPUT "verkletactpta 79 eK 
140 PI=4XATN(I) 
150 U=128:V=96:H=. S:RD=PI/180 
160 C=K&COSCWERD) z: S=KSINCWERD) 
170 DASS s DY=5 © Kil=t5 
180 DIM H(255) 
190 FOR L=0 TO 255 
200 H(L)=1000 
210 NEXT L 
220 CLS 
230 COLOR 1,15,4 : SCREEN 2 
240 FOR YY=-95 TO 95 STEP DY 
250 MI=COS(YY*2*PI/95-PIJ+1 
260 FOR XX==95 TO 95 STEP DX 
270 M2=COS(XXK24PI/GS-PIJ+] 
280 ZaK1&MI*M2 
290 XG=INTCUIXXHCHYY+H) 
300 YG=INT(V=-S*YY=Z +H) 
310 IF XG<O THEN XG=0 
320 IF XG>255 THEN XG=255 
330 IF XX==95 
THEN GOSUB 1000 
ELSE GOSUB 2000 
340 NEXT XX 
350 NEXT YY 
360 AS=INPUT$S(1) 
370 COLOR 15,4,4 : END 
380 * 
1000 REM subroutine eerste xx=waarde 
1010 F1=0 : L=INT(XG/DX) 
1020 IF YGS<=H(L) THEN F1=1:H(L)=YG 
1030 X1=XG:Y1=YG 
1040 RETURN 
2000 REM subroutine volgende xx'en 
2010 F2=0 : L=INT(XG/DX) 
2020 IF YG<=H(L) THEN F2=1:H(L)=YG 
2030 X2=XG:Y2=YG 
2040 IF F1*F2=1 
THEN LINE (X1,Y1)=(X2,Y2) 
2050 x1=X2:Yi=Y2:Fl=F2 
2060 RETURN 
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Voor de onderstaande tekening kiest u W=45°®, K=0,5 en K1=15 of 
20. 


EL EE 2ny _ 
z = (cos(a 1) +Meos Te m)+1) 
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6 Turtle-graphics en 
LOGO-simulatie 


LOGO is vooral door de turtle-graphies beroemd geworden. Er is 
geen andere programmeertaal waarmee zo eenvoudig zeer moeilijke 
grafische figuren gemaakt kunnen worden dan met LOGO, In dit 
hoofdstuk zullen we vijf LOGO-programma's in BASIC vertalen. Deze 
BASIC-programma's bevatten wel 25 tot 30 regels, terwijl LOGO 
hiervoor slechts 3 of 4 regels nodig heeft. 


"Eerst LOGO, daarna andere programmeertalen" 


is het motto van vele informatici en pedagogen die zich met de 
invoering van de informatica op scholen bezighouden. 


Wanneer men deze stelling wil begrijpen, moet men niet alleen de 
taal LOGO en haar mogelijkheden, maar ook de omgeving waarin 
LOGO ontwikkeld werd goed kennen. Een korte historische terug- 
blik lijkt daarom op zijn plaats. 


Het schrijven van computerprogramma's is een intellectuele presta- 
tie en een creatief proces. Aan de wijze waarop wij programma's 
ontwikkelen, herkennen wij direct onze manier van denken. Aan 
het Massachusetts Institute of Technology (MIT) heeft men al in het 
begin van de zestiger jaren een speciale programmeertaal ontwik- 
keld, met behulp waarvan men kunstmatige intelligentie ging bestu- 
deren. Men noemde deze taal LISP, een afkorting van LISt- 
Processing. 


In LISP werden programma's geschreven die een met kunstmatige 
zintuigen uitgeruste elektromechanische muis in staat stelden een 
uitweg uit ieder willekeurig doolhof te vinden. In LISP worden pro- 
gramma's ontwikkeld die de computer tot een medespeler met leer- 
vermogen maakt. Hoe meer spelletjes de computer tegen een mense- 
lijke tegenspeler speelt, des te beter wordt hij. Goede zetten van 
de tegenstander onthoudt hij en legt hij vast in zijn geheugen, ter- 
wijl hij de eigen slechte zetten voor toekomstige spelletjes uit zijn 
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geheugen wegveegt. Een 'afvalprodukt' van deze studies aan het 
MIT zijn de moderne schaakprogramma's. 


LOGO is een hoog ontwikkeld dialect van de LISP-taal. Seymour 
Papert, leerling van de bekende onderzoeker Jean Piaget en mede- 
werker aan het MIT, heeft in twaalf jaar tijd LOGO ontwikkeld. Jean 
Piaget heeft in zijn bekende boek Hoe kinderen leren de kinderlijke 
denkstructuren laten zien. Hieruit blijkt dat tekenen en knutselen 
tot de eerste creatieve handelingen van kinderen behoren. Wij zul- 
len zien hoe LOGO deze bezigheden ondersteunt. 


Amerikanen gaan door voor een volk dat graag experimenteert. Het 
was hun echter duidelijk, dat jonge kinderen, wij denken dan aan 
zes- tot dertienjarigen, niet in staat zijn een programmeertaal 
zoals BASIC, laat staan Pascal, te leren en algorithmen voor nume- 
rieke, niet-numerieke of grafische problemen te schrijven. Daarom 
heeft men een 'kinderlijke' taal gemaakt (waarachter echter een 
imponerende software schuilgaat), waarmee het kind met gemak 
tekeningen kan maken. 


Wordt het LOGO-systeem ingeschakeld, tegenwoordig meestal een 
64K-microcomputer met bijbehorende software, dan verschijnt in 
het midden van het beeldscherm een kleine driekhoek, waarvan de 
top naar boven wijst. De kinderen noemen die driehoek 'turtle', het 
Engelse woord voor schildpad. Met behulp van eenvoudige bevelen 
kan het kind de schildpad willekeurig over het beeldscherm laten 
rondlopen. En al naar gelang het wenselijk is laat de turtle wel of 
geen spoor na. 


De volgende LOGO-opdrachten zijn beschikbaar : 


FORWARD 100 = 100 passen vooruit 
BACK 50 = 50 passen achteruit 


RIGHT 90 = draaiing van 90° met de klok mee 
LEFT 45 = draaiing van 45° tegen de klok in 
PENUP = haal pen van papier 

PENDOWN = zet pen op papier 

HIDETURTLE = haal Turtle van het beeldscherm 


Normaal gesproken geldt de toestand "PENDOWN", Als we bijvoor- 
beeld de hoofdletter F op het scherm willen tekenen, kan dat in 
LOGO als volgt: 


FORWARD 100 RIGHT 90 FORWARD 50 
RIGHT 90 PENUP FORWARD 50 PENDOWN 
RIGHT 90 FORWARD 50 

HIDETURTLE 


Voor 100, 50, 90 en 45 kunnen ook andere waarden gekozen worden. 
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Omdat LOGO een vertolkend systeem is, wordt elke opdracht (na 
het geven van RETURN) direct uitgevoerd. Zo op het oog lijkt het 
slechts leuk speelgoed! Laten we nu eens een LOGO-programma 
bekijken; liever spreken we van een procedure: 


TO VIERKANT 
REPEAT 4 (FORWARD 75 RIGHT 90) 
END 


Het LOGO-systeem weet dat tussen TO en END een procedure (een 
stuk programma) gedefinieerd wordt. De procedure tussen TO en 
END (in ons voorbeeld VIERKANT genoemd) maken we zelf. Zouden 
we bovenstaande procedure intoetsen, dan zal op het beeldscherm 
een vierkant met zijden ter lengte 75 getekend worden. De program- 
meurs onder u kennen vast de opdracht REPEAT waarmee een her- 
halingsstructuur geprogrammeerd kan worden. 


Als we de procedure VIERKANT in het LOGO-systeem ingevoerd 
hebben, kunnen we hier ook gebruik van maken. LOGO onthoudt 
alle procedures die wij zelf invoeren. We kunnen dan ook later 
nieuwe procedures ontwikkelen waarin we gebruik maken van eer- 
der ingevoerde procedures (zoals VIERKANT). Nu volgt een proce- 
dure met een parameter: 


TO VIERKANT :ZIJDE 
REPEAT 4 (FORWARD:ZIJDE RIGHT 90) 
END 


Toetsen we nu VIERKANT 150 <RETURNP> in, dan zal LOGO een 
vierkant met zijden van 150 tekenen. Niets verhindert u een proce- 
dure in te bedden in een volgende procedure en die weer in een 
volgende en die weer .….……. ‚ enzovoorts. Alleen de beschikbare 
geheugenruimte is hierbij een remmende factor. Bekijk het volgende 
LOGO-programma. 


LOGO-programma nr. 1 


TO VIERKANTPATROON 
REPEAT 8 (FORWARD 20 LEFT 45 VIERKANT 75) 
END 


Dit eenvoudige drieregelige LOGO-programma tekent een patroon 

van acht vierkanten (zie p.87). We zullen dit programma straks in 
BASIC vertalen. Het zal blijken dat hiervoor enige wiskundige en 
programmeertechnische kennis nodig is. Het LOGO-programma kan 
door een leerling van de lagere school gemaakt worden, terwijl het 
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BASIC-programma dat deze acht vierkanten tekent slechts door 
scholieren van de hoogste klassen van het voortgezet onderwijs kan 
worden gemaakt. Het wordt nog moeilijker als we intikken: 


TO STROOK:AANTAL 
REPEAT :AANTAL (FORWARD 100 VIERKANTPATROON) 
END 


Tikken we vervolgens in STROOK 5 <RETURN> dan maken we 

heel eenvoudig een fraai ornament. Probeer dit maar eens in BASIC 
of Pascal. LOGO wordt pas echt interessant als we van de mogelijk- 
heid gebruik maken dat een procedure zichzelf aanroept (recursivi- 
teit). De turtle-graphies berusten op het principe dat we een pro- 
cedure parameters geven en dat deze procedure zichzelf steeds met 
andere parameterwaarden aanroept. Bekijk de volgende procedure: 


TO VEELHOEK:ZIJDE : HOEK 

FORWARD:ZIJDE LEFT: HOEK 

VEELHOEK:ZIJDE:HOEK —- hier roept de procedure 
END zichzelf aan 


Als we nu VEELHOEK 200 144 intikken, zal het LOGO-systeem een 
regelmatige vijfpuntige ster tekenen. We zullen echter op de STOP- 
toets moeten drukken om het tekenen te stoppen. Brengen we in 
deze procedure een stopmechanisme aan en laten we de procedure 
steeds de waarde van ZIJDE veranderen, dan krijgen we de wel- 
haast bekendste LOGO-procedure: 


LOGO-programma nr.2 


TO TURTLE:ZIJDE:GROEI:HOEK 
FORWARD:ZIJDE LEFT :HOEK 
MAKE:ZIJDE : ZIJDE+: GROEI 
IF:ZIJDE>200 THEN STOP 
TURTLE:ZIJDE:GROEI: HOEK 
END 


Als er een dubbele-punt voor een LOGO-woord staat, weet het 
LOGO-systeem dat het om een naam van een variabele gaat en niet 
om de naam van een procedure of een LOGO-opdracht. In de boven- 
staande LOGO-procedure zien we hoe de procedure zichzelf steeds 
met een nieuwe waarde voor ZIJDE aanroept. Ook dit programma 
vertalen we in BASIC. Hiermee kunt u elke denkbare rechtlijnige 
turtle-grafiek maken. U hoeft alleen de waarden van de invoerpara- 
meters (ZIJDE, GROEI en HOEK) te veranderen. 
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Het derde LOGO-programma tekent een reeks vierkanten in spiraal- 
vorm. Het zijn bekende figuren in het 'land van de computer- 
graphics'. 


LOGO-programma nr. 3 


TO VIERKANTSPIRAAL:ZIJDE : HOEK 
IF:ZIJDE>150 THEN STOP 
VIERKANT :ZIJDE LEFT :HOEK 
VIERKANTSPIRAAL:ZIJDE+4: HOEK 
END 


Ook dit LOGO-programma vertalen we in BASIC. Met deze drie pro- 
gramma's hebt u niet alleen met LOGO maar ook met de 'turtle- 
graphics! kennis gemaakt. Wilt u meer weten van deze graphics- 
wereld, dan wijzen wij op het boek Turtle Geometry: The Computer 
as a Medium for Exploring Mathematics van Harold Abelson en 
Andrea di Sessa, uitgegeven door Cambridge, MA:MIT Press, 1981. 
Liefhebbers van wiskunde, informatica en LOGO vinden in dit boek 
een schat aan informatie. 


Natuurlijk zijn we in LOGO niet gebonden aan ‘rechte lijnen'. In 
LOGO hebben we de beschikking over alle wiskundige functies en 
bewerkingen. Alle programma's uit dit boek kunnen in LOGO 
geschreven worden. Zij zouden dan vast korter, overzichtelijker 
en begrijpelijker geweest zijn. In LOGO kunnen we heel eenvoudig 
met lijsten en tabellen werken. In LOGO hoeven we de variabelen 
niet te declareren. Achter een naam in LOGO kan gewoon een getal, 
een string, een n-dimensionaal veld en nog veel meer schuil gaan. 
Ook kan het de naam van een procedure zijn. We hoeven ons niet 
te bekommeren om het reserveren van geheugenruimte (DIM- 
opdracht in BASIC), noch om het aangeven van het type (real, 
integer, character, boolean, enz.); het LOGO-systeem regelt dit 
zelf. 


Wie in LOGO programmeert leert automatisch netjes programmeren. 
Begrippen als procedure, recursie, herhaling, toekenning worden 
op een natuurlijke manier aangeleerd. Met een paar opdrachten is 
een kind al in staat ingewikkelde figuren te tekenen. Met LOGO 
degraderen we de computer niet tot een programmeerbare zakreken- 
machine, maar halen we er alles uit wat erin zit. 
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Vertaling LOGO-programma nr.1 in BASIC 


De kern van elk LOGO-graphies-programma wordt gevormd door 
twee opdrachten: 


FORWARD (resp. BACK) :ZIJDE 
LEFT (resp. RIGHT) :HOEK 


Hoe vertalen we deze opdrachten in BASIC? Bekijk hiertoe de vol- 
gende figuur. 


HOME 


DX=S"COS(W) 
DY=S"SIN(W) 


De PEN staat in punt Pi(x1,y:) en moet in de richting van hoek W 
over een afstand S verplaatst worden. In P:(x2,y2) aangekomen 
moet de PEN (eigenlijk de turtle) een hoek van DW° linksom maken. 
De volgende BASIC-opdrachten voeren dit uit: 


INPUT "coordinaten start '";X1,Y1 
INPUT "“beginrichting(gr) '';W 
PI=4&XATN(1) 
H=.5:RD=PI/180:WI=WERD 
X2=INT(XI+STCOSCWIJ +H) 
Y2=INT(Y1=-STESINCWI)+H) 

LINE CENXCXI),YII=CENK CHI, Y 2) 
xi=X2:yl=Y2 

WaW+DW 

IF W>360 THEN W=W-360 

WI=W*RD 


Deze opdrachten komt u in alle volgende BASIC-programma's tegen. 
Meer theorie is niet nodig! De programma's moeten met bovenstaan- 
de informatie gelezen kunnen worden. 
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Experimenteer met programma 31 LOGO-1 vierkantpatroon. Met 
W=90°, S1=20, DW=45, S2=75 en N=8 krijgen we de bovenstaande 
tekening. 


100 REM programma 31 vierkantpatroon 
(LOGO=1) 

110 ENPUT "coordinaten start "sx1,Y1 

120 INPUT "“beginrichting(gr) ";W 


130 INPUT "verplaatsing 091 
140 INPUT "draaihoek Linksom "';DW 
150 INPUT "zijde vierkant VEGË 


160 INPUT "aantal vierkanen "';N 

170 PI=4XATN(I) 

180 H=.5:RD=PI/180:WI=WERD 

190 CLS 

200 COLOR 1,15,4 : SCREEN 2 

210 DEF FNX(XD=INT(374X/1.4+.5) 

220 FOR J=1 TO N 

230 X2=INT(X1+STCOSCWI)+H) 

240 Y2=INTCY1-STSINCWI)+H) 

250 LINE (FNXCX1),YI)=-CFNXCX2) ,Y2) 

260 x1=X2:Y1=Y2 

270 AX=X1:AY=Y1 

280 W=W+DW 

290 IF W>360 THEN W=W-360 

300 WI=W*RD 

310 FOR K=0 TO 2 

320 X2=INTCXT+S2COSCWI+EK/24PIJ+H) 
330 Y2=INTCYI=S2*SINCWI+K/24PIJ+H) 
340 LINE CEFNX(X1),Y1)=CFNX(X2) ,Y2) 
350 xi=X2:Y1=Y2 

360 NEXT K 

370 LINE CENXCTD , VANS CENKCAKD „AND 

380 X1=AX:Y1=AY 

390 NEXT J 

400 A$S=INPUTS$S(1) 

410 COLOR 15,4,4: END 
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W=90, S1=30, DW=48, S2=70 en N=15 geeft de onderstaande tekening. 
Als DW een deler is van 360, is N gelijk aan 360/DW. Kiest u voor 
DW een willekeurige waarde, neem dan voor N een groot getal, bij- 
voorbeeld 100, en breek het tekenen met de stoptoets af. De 'kin- 
deren! doen dat in LOGO ook op deze manier. 


ZE 
ET DN Ad 
q o. ee LS 


Vertaling LOGO-programma nr.2 in BASIC 


Het programma lijkt sterk op het vorige en behoeft daarom geen 
commentaar. 


100 REM programma 32 turtle-grafiek 
(LOGO=2) 

110 INPUT "coordinaten start '";X1,Y1 

120 INPUT "beginrichting(gr) '";W 


130 INPUT "verplaatsing WES 
140 INPUT "draaihoek Linksom ';DW 
150 INPUT "toename zijde SDS 


160 PI=4XATN(1) 
170 H=.5:RD=PI/180:WI=WERD 
180 CLS 
190 COLOR 1,15,4 7 SCREEN 2 
200 DEF FNX(X)=INT(37+X/1.4+.5) 
210 X2=INT(XT+S*COSCWI) +H) 
220 Y2=INTCYI=S*SINCWI)+H) 
250 IF X2<0 OR KES255 OR 
Y2<0 or Y2>191 
THEN AS=INPUTS$(1) :STOP 
240 LINE CFNXCXI),YI)=(FNX(X2) ,Y2) 
250 x1=X2:Yl=Y2 
260 W=W+DW 
270 IF W>360 THEN W=W=360 
280 WI=W&RD:S=S+DS 
290 GOTO 210 
300 AS=INPUT$(1) 
310 COLOR 15,4,4: END 
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Het is handig de ingevoerde waarden na het tekenen op het scherm 
of op een printer te laten afdrukken. In de volgende tabel zien we 
de diverse waarden voor de daaronderstaande tekeningen. 


EEE 


1 
96 
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Vertaling LOGO-programma nr.3 in BASIC 


100 


110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 


250 
260 
270 
280 
290 
300 
310 
320 
330 
1000 
1010 
1020 
1030 


1040 


REM programma 33 vierkantspiraal 
(LOGO=-3) 

INPUT "coordinaten start ";X1,Y1 

INPUT "beginrichting(gr) '";W 


INPUT "verplaatsing mn 
INPUT "draaihoek Linksom ';DW 
INPUT "toename zijde “DS 


PI=4XATN(1) 
H=.5:RD=PI/180: WI=WERD 
CLS 
COLOR 1,154 1 SCREEN 2 
DEF FNX(X)=INT(37+X/1.4+.5) 
AX=X1:AY=Y1 
FOR FASE=0O TO PI STEP PI/2 
GOSUB 1000 ‘punt (X2,Y2) 
IF VLAG=| 
THEN AS=INPUT$S(1): 

COLOR 15,4,4:STOP 
LINECFNX(XI) „YI =-CENXCX2) „Y2) 
xi=X2:Y1=Y2 

NEXT FASE 
LINE (FNX(X1),Y1)=(FNXCAX),AY) 
X1=AX:Y1=AY 
WEW+DW:IF W>360 THEN W=W-360 
WI=W&RD:ES=S+DS 
GOTO 220 
END 
REM subroutine (X2,Y2) 
X2=INT(X1+S*COSCWI+FASE)+H) 
Y2=INT(Y1=-S*SINCW1I+FASE)+H) 
IF X2<0 OR X2>255 OR 
Y2<0 OR Y2>191 
THEN VLAG=1 
RETURN 
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Voor de boven het programma staande tekening kozen we X1=128, 
Y1=96, W=90, S=5, DW=15 en DS=3,5. 


Vergelijk deze drie BASIC-programma's eens met de overeenkom- 
stige LOGO-programma's. In BASIC moeten we zelf alle graphic- 
software schrijven, terwijl deze in LOGO als machine-routines aan- 
wezig is. 


Het 4e en 5e LOGO-programma 


Tot slot van dit hoofdstuk geven we nog twee BASIC-programma's 
waarmee de turtle behalve rechte wegen ook kromme wegen kan 
bewandelen. Dit voegt een nieuw element aan de turtle-graphics 
toe. De LOGO-structuur zullen we stapsgewijs verfijnen. Als uit- 
gangspunt nemen we de onderstaande figuur. Deze bestaat uit een 
lijnstuk L, met daaraan vast een cirkel met straal R. We noemen 
deze figuur voor het gemak even 'de steelpan’. 


Als de turtle de steelpan tekent, begint hij in A. Hij legt vervol- 
gens een afstand L in positieve x-richting af en komt in B. Hier 
draait hij 90° met de klok mee en legt daarna de omtrek van de cir- 
kel af tot hij weer in B uitkomt. Daar draait hij zich 90° met de 
klok mee en kijkt dan rechtvooruit naar het punt A. Met deze steel- 
pan kunnen we twee soorten tekeningen maken. Laten we deze 
mogelijkheden bekijken : 


LOGO-programma 4 


TO CIRKELFIGUUR-1I:LIJNSTUK:STRAAL:HOEK 
STEELPAN :LIJNSTUK :STRAAL 

LEFT : HOEK 
CIRKELFIGUUR-1:LIJNSTUK:STRAAL:HOEK 
END 
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LOGO-programma 5 


TO CIRKELFIGUUR-2:LIJNSTUK:STRAAL: HOEK 
REPEAT 4 (STEELPAN:LIJNSTUK:STRAAL LEFT 90) 
LEFT : HOEK 
CIRKELFIGUUR-2:LIJNSTUK:STRAAL: HOEK 

END 


Deze twee LOGO-programma's kunnen nog niet uitgevoerd worden, 
want de procedure STEELPAN is nog niet in het LOGO-systeem 
gedefinieerd. We gaan de structuur verfijnen: 


TO STEELPAN:LIJNSTUK:STRAAL 
FORWARD:LIJNSTUK RIGHT 90 
CIRKEL:STRAAL 

RIGHT 90 

END 


Ook in deze procedure zien we een nog niet gedefinieerde proce- 
dure, CIRKEL, opduiken. Om de turtle een cirkel te kunnen laten 
maken zullen we deze procedure CIRKEL moeten maken. Veel LOGO- 
versies bevatten reeds zo'n voorgedefinieerde CIRKEL-procedure. 
MSX-BASIC beschikt ook over een CIRCLE-functie. We zien deze in 
regel 290 van het volgende programma. Bedenk dat de X-coördinaat 
van het middelpunt met de FNX-functie moet worden aangeroepen, 
anders krijgen we een ellips. 
correctiefactor 
CIRCLECFNXCUD „VD ,R‚1,0,2*PI,1.4 * voor de horizon- 
ef | k a tale straal 
coördinaten straal |} begin bind tekenhoek 
middelpunt tekenhoek 
afdrukkleur 
(1=zwart) 


Als we echt LOGO in MSX-BASIC willen nadoen, moeten de cirkels 
getekend worden vanuit het uiteinde van de steel (zie punt B in 
tekening op p.92). De MSX-CIRCLE-functie begint met tekenen aan 
de andere kant van de cirkel. 


Als we het 'echt' willen doen moeten we de regels 290 t/m 320 ver- 
vangen door: 
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AX=X2 : AY=Y2 

REM CIRKEL M(U,V) EN STRAAL R 

UZINT(AX+R*COS CWI) +H) 

VZINT(AY=R&SINCWID+H) 

X1=AX : Y1=AY 

FOR WW=C(W1=-PI) TO (W1+PI) STEP PI/32 
X2=INT CUERKCOSCWW) +H) 
Y2=INT(V=-R&SINCWW)+H) 
LINE (CFNX(X1), YT) =(ENX(X2) ,Y2) 
xi=X2:Yi=Y2 

NEXT WW 

WEW+DW : IF W>=360 THEN W=W-360 

WI=WERD 7 X1=AX : YI=AY 


Vertaling van LOGO-programma nr.4 in BASIC 


100 REM programma 34 cirkeliguur=1 
(LOGO=4) 

110 CLS 

120 INPUT "coordinaten start ";X1,Y1 

130 INPUT "beginrichting(gr) '';W 

140 INPUT "lengte steel in 

150 INPUT "straal cirkel WeR 

160 INPUT "draaihoek Linksom '';DW 

170 PI=4KATN(1) 

180 H=.5:RD=PI/180:WI=WERD 

190 CLS 

200 COLOR 1,15,4 : SCREEN 2 

210 DEF FNX(X)=INT(37+X/1.4+.5) 

220 xx=X1:YY=Y1 

230 REM begin herhaal=lus 

240 X2=INT(XT+L#COS CWI) +H) 

250 Y2=INT(Y1=L&SINCWI)+H) 

260 LINE GENWCRTD VAN CENK CHR VI 

270 U=INT(X2+R*COS CWI) +H) 

280 V=INT(Y2=R*SINCWI)+H) 

290 _ CIRCLECFNX(UI,V),R„1,0,2*PI,1.4 

300 WEW+DW:IF W>=360 THEN W=W=360 

310 WI=W&RD 

320 Xl=X2:Y1=Y2 

330 IF X1<>XX OR YTI<>YY THEN 240 

340 REM einde herhaal=lus 

350 AS$S=INPUTS$(1) 

s60 COLOR Tbr 7 END 


De bovenstaande tekening krijgen we door X1=128, Y1=135, W=90, 
L=90, R=20 en DW=170 te kiezen. 
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Vertaling LOGO-programma nr.5 in BASIC 


100 REM programma 35 cirkeliguur=-2 
(LOGO=5) 

VU IGS 

120 INPUT "coordinaten start ";X1,Y1 

130 INPUT "beginrichting(gr) ";W 

140 INPUT "lengte steel viel 

150 INPUT "straal cirkel eN 

160 INPUT "draaihoek Linksom '';DW 

170 PI=4KATN(1) 

180 H=.5:RD=PI/180:WI=WERD 

120: ELS 

200 COLOR 1,15,4 : SCREEN 2 

210 DEF FNX(X)=INT(37+X/1.4+.5) 

220 REM begin eindloze Lus 

230 FOR J=1 TO 4 

240 X2=INTCXI+L*COS CWI) +H) 

250 Y2=INTCYT=L&SINCWI)+H) 

260 GENE GENK CRT NTI CENK CHAIN 2D 

270 REM cirkel met middelpunt 
(U,V) en straal R 

280 U=INT(X2+RKCOS CWI) +H) 

290 V=INT(Y2=R&SINCWI)+H) 

300 GIRCLECFNX CUI VI Rs TO 2PI, 1-4 

310 W=W+90 : IF W>=360 THEN W=W-360 

320 WI=W*RD 

330 x1=X2:Y1=Y2 

340 NEXT J 

350 W=W+DW 

360 WI=W&RD 

370 GOTO 230 

380 REM 'einde' eindloze Lus 

390 END 


Voor deze tekening geldt: X1=128, Y1=96, W=0, L=50, R=10 en 


DW=45. 
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Met andere waarden voor deze variabelen kunt u de mooiste turtle- 


plaatjes maken. 


7 Educatieve toepassings- 
programma’s 


De voorgaande 35 grafische programma's zijn geen toepassingspro- 
gramma's. Ze waren bedoeld om een overzicht te geven van de 
mogelijkheden van graphies met hoog oplossend vermogen. In dit 
hoofdstuk zullen we vijf praktijkgerichte grafische programma's 
presenteren. Deze programma's zijn: 


1 Tekenen van een landkaart 

2. Maken van een histogram (een stavengrafiek ) 

3. __Demonstratieprogramma voor de breking van lichtstralen 
: Demonstratieprogramma voor de 'speldenworp' van Buffon 
5 Prooi-roofdierpopulaties 


1. Tekenen van een landkaart 


Educatieve programma's over topografie gebruiken dikwijls land- 
kaarten van een land of van een werelddeel die door de computer 
getekend worden. We zullen laten zien hoe de computer zo'n land- 
kaart kan tekenen. We gaan de landkaart (althans de grensomtrek ) 
van Zwitserland tekenen. Waarom Zwitserland? Wel, in de eerste 
plaats omdat de auteur dit gekozen heeft en in de tweede plaats 
omdat Zwitserland geen eilanden heeft. Het tekenen van de contou- 
ren van Nederland gaat in principe net zo, alleen kosten al die 
eilanden een hoop extra werk, vandaar! 


Hoe tekenen we nu de omtrek van een bepaald land? Heel eenvoudig! 
We pakken gewoon een atlas, een stukje overtrekpapier, een pot- 
lood, een lineaal en millimeterpapier. We kiezen in de atlas een land 
of werelddeel uit en trekken het op overtrekpapier over. Daarna 
leggen we het overgetrokken plaatje op millimeterpapier, waar we 
van te voren een coördinatenstelsel op getekend hebben (een x- 

en een y-as). We bepalen nu van een (groot) aantal punten op de 
omtrek van het land (of werelddeel) de coördinaten (x,y) en schrij- 
ven deze op. Als we ons op het beeldscherm ook een coördinaten- 
stelsel voorstellen en als we hierin de punten, waarvan we de coör- 
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dinaten in een programma opnemen, met elkaar laten verbinden, dan 
ontstaat een 'landkaart' op het scherm. 


In het onderstaande programma, dat 'de kaart' van Zwitserland 
tekent, is een 'echte' kaart gebruikt met een schaal van 

1 : 2.000.000. Om een enigszins natuurgetrouwe weergave te ver- 
krijgen zijn de coördinaten van 90 grenspunten berekend. De coör- 
dinaten van deze punten, die opgegeven zijn in millimeters ten 
opzichte van de oorsprong van het gekozen coördinatenstelsel, zijn 
in DATA-regels opgenomen. Als het programma het coördinaten- 
paar (0,0) leest, weet het dat de 'kaart' af is. 


De x-coördinaten liggen tussen 6 en 177. Om deze naar het hoge- 
resolutiebereik 0-255 en 0-191 te transformeren vermenigvuldigen 
wij de x- en y-coördinaat met 1,35. Hierdoor blijft er nog ruimte 
over voor bijvoorbeeld een kader rond de kaart. De werking van 
het programma zal hiermee duidelijk zijn. 


We hebben op deze manier ook een kaart van Europa en zelfs een 
wereldkaart getekend. Het probleem met de eilanden hebben we als 
volgt opgelost. Als het programma in een DATA-regel negatieve x- 
en y-coördinaten leest, weet het programma dat dit punt niet met 
het vorige verbonden moet worden en dat dit punt dus het begin 
is van een apart stukje 'land'. 


ZWITSERLAND 
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100 REM programma 36 kaart 
van zwitserland 
1D: CLS 
120 COLOR 1,15,4 : SCREEN 2 
130 DEF FNX(X)=INT(37+X/1.4+.5) 
140 K=1.35:H=.5:U=18:V=191 
150 READ X,Y 
160 X1=INTCU4K*X+H) 
170 YIZINT(V=K«Y+H) 
180 REM begin herhaal Lus 


190 X2=INT(U4K*X+H) 

200 Y2=INT(V=K*Y+H) 

210 LINE CENX(X1), YI) -CFNX(X2) ,Y2) 
220 X1=X2:Y1=Y2 

230 READ X,Y 

240 IF X<> O0 THEN 190 


250 REM einde herhaal Lus 

260 AS=INPUT$S(1) 

270 COLOR 15,4,4 3 END 

280, 

290. DATA -69,108,71,107, 70, 104,75 ,1,04 
300 DATA 75,104,76,106,80,107,81,104 
310 DATA 86,105,91,108,94,107,101,106 
320 DATA 100,108,105,108,106,110,101 
330 DATA 110,98, 712,102,117, 108, 118 
360 DATA Wie, 112 1411511811 01 28 
550. DATA 110713910251 655 105, 146,95 
360 DATA 142,86,144,78,154,77,154,77 
370 DATA 154,72,163,67,168,68,173,76 
380: DATA VIT Tar TUN Srl 15e 
390 DATA 1/1,51,164,56,164,50,165, 43 
400 DATA 166,34,162,34,157,42,143,36 
410 DATA 139,48,136,45,133,48,132,40 
420 DATA 12de leortostezel kelken 
450 DATA 114,20,116,25,100,35 102,45 
440 DATA 94,42,88,35,90,28,79,15,75 
650. DATA 106,19, 600,14 ese la 
460 DATA 37,29,39,36,37,40,39,43,29 
470 DATA 45,16,38,18,33,13,29,6,28,6 
480. DATA Sez11,34,135,40-10,455 12755 
690 DATA 25,6326,73,50,73,485 94 rbe 
200 DATA P4546,102,:54, 102,55, 9961 
510 DATA 98,63,102,69,108 

520 DATA 0,0 


Zie de appendix voor een uitbreiding van dit programma. 
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2. Maken van een histogram 


Vaak willen we een aantal waarnemingen grafisch in een stavendia- 
gram of histogram weergeven. Hiertoe dient het volgende programma. 
Dit programma tekent een horizontale en een verticale as. De verti- 
cale as wordt in stukjes van acht schermpuntjes verdeeld. Dit komt 
overeen met 5% van de hele verticale as. Elk tweede streepje op 

deze as geeft de volgende 10% aan en is iets breder getekend. Hier- 
mee kan de lengte van de staven redelijk geschat worden. 


Het programma kan maximaal 100 waarnemingen verwerken. Uit este- 
tische overwegingen moet u echter niet veel meer dan 40 waarnemin- 
gen invoeren, omdat de staafjes anders meer op lijntjes dan op balk- 
jes gaan lijken. 


De waarnemingen worden zo 'geschaald' dat de grootste waarneming 
precies met 160, de lengte van de verticale schaalverdeling, over- 
eenkomt. 


Het zou mooi zijn als we bij de staven, de assen en de schaalverde- 
ling tekst en getallen konden zetten, maar dat gaat heel lastig, dus 
doen we het nu maar niet. Denkt u erom dat de getallen die u 
intoetst de lengte van de staven aangeven. Als elke staaf een 
bepaalde klasse met waarnemingen voorstelt, dan voert u dus steeds 
het aantal waarnemingen (de frequentie) van een klasse in. Het pro- 
gramma kan uitgebreid worden door er een stuk voor te zetten dat 
ruwe gegevens inleest, die vervolgens netjes in klassen verdeeld 
worden. De frequentie van de waarnemingen in de klassen wordt 
vervolgens gebruikt om de stavengrafiek te tekenen. 
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100 REM programma 37 histogram 

110 CLS 

120 PRINT "histogram tekenen! 

130 PRINT "======vememme ke 

140 PRINT:PRINT 

150 INPUT "hoeveel gegevens(<40)'';N 

160 DIM ACN) 

170 MX==10000000000! : PRINT 

180 FOR J=1 TO N 

190 PRINT "waarde"; J; TAB(12); 

200 INPUT A(J) 

210 IF ACJ)>MX THEN MX=ZA(J) 

220 NEXT J 

250 GES 

240 COLOR 1,15,4 : SCREEN 2 

250 DEF FNXCH)=INTC3A4K/1.4t. 5) 

260 REM horizontale as 

270 LINE (FNX(O),180)=-(FNX(256),180) 

280 REM verticale as 

290 LINE (FNX(10),0)=CFNX(10),180) 

300 REM schaalverdeling 

310 FOR J=10 TO 1 STEP =1 

320 X1=7:Y1=180-J*16:X2=10: Y2=Y1 

330 LINE CENKCRTI,NYTISLENK: CEI 5 Y ZI 

340 X1=4:Y1=Y2+81:X2=10: Y2=Y1 

350 LINE (FNXCX1), YI) =-(FNX(2) ,Y2) 

360 NEXT J 

370 REM staven tekenen (B=breedte) 

380 B=INT(240/N) 

390 FOR J=1 TO N 

400 X1=(J-1)#B+15:Y1=180 

410 X2=X1: 
Y2=INT(180-160%A(J) /MX+H) 

420 LINE CENXCKII YT CENKCK2I 7 Y 2) 

430 Xl=X2:Yl=Y2:X2=X1+B: Y2=Y1 

440 LINE (FNX(X1),Y1)=(FNX(X2) ,Y2) 

450 Xl=X2:Yi=Y2:X2=X1: 22180 

460 LINE CFNXCXI) „YI =CENXCH2) ,Y2) 

470 NEXT 

480 AS$S=INPUTS(1) 

490 COLOR 15,4,4 : END 
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3. Demonstratieprogramma voor de breking van lichtstralen 


Met dit programma willen we laten zien hoe de MSX-computer bij 
natuurkundelessen gebruikt zou kunnen worden. Voordat we het 
programma geven leggen we nog iets uit van de natuurkundige 
beginselen van de breking van licht. 


Als een lichtstraal vanuit de ene stof (bijvoorbeeld lucht) een ande- 
re, optisch dichtere, stof (bijvoorbeeld water) binnenkomt, wordt 
de straal op het scheidingsvlak van beide stoffen naar de verticale 
as toe gebogen. Zie onderstaande figuur. 


[Verticale Werticale 
|as pe 
| 


Medium I 
(lucht) 


| 
| Medium II 
| (water) 


ej en eg zijn de lichtsnelheden in respectievelijk de eerste en de 
tweede stof. De constante n heet de brekingsindex. Bij de overgang 
van lucht naar water geldt een brekingsindex van 1,33. Voor de 
overgang van lucht naar glas gelden andere waarden, enzovoorts. 


Als het licht vanuit een 'dichter' medium overgaat in een 'dunner!' 
medium geldt het omgekeerde: de lichtstralen worden nu van de 
verticale as, die loodrecht op het scheidingsvlak van beide stoffen 
staat, afgebogen. De brekingshoek a kan nooit groter dan 90° zijn. 


Voor dit grensgeval (a = 90°) geldt: 


(want sin90° = 1) 
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Bij de overgang van lucht naar water geldt voor B*: 
sin'=7z = sing=0,7519 = B*=48,75° 


Groter dan 48,75° kan B dus niet worden. Dit betekent dat als het 
licht van water overgaat in lucht met een invalshoek a die groter is 
dan B* (zie rechter figuur), het licht niet meer het water 'uitkomt'. 
Op het scheidingsvlak van water en lucht wordt het licht geheel 
teruggekaatst. Op dit principe berusten de moderne glasvezelkabels, 
waarin informatie in de vorm van licht wordt getransporteerd. 


In het onderstaande demonstratieprogramma kan de brekingsindex 
n ingetoetst worden. De lichtbron bevindt zich in het punt met 
schermcoördinaten (0,190). Het scheidingsvlak ligt horizontaal en 
is de lijn V=96. De invalshoek van een lichtstraal die van medium 2 
(de dichtere stof, onderste helft) in medium 1 (de lichtere stof, 
bovenste helft) overgaat wordt van 0° steeds met stapjes van 3° 
opgehoogd. Op het moment dat deze invalshoek groter wordt dan 6* 
(deze is afhankelijk van de ingetoetste brekingsindex) treedt totale 
reflectie (terugkaatsing) op. Rond de figuur wordt een kader gete- 
kend. De eindpunten van de diverse lichtstralen worden met trigo- 
nometrische functies berekend. De commentaaropdrachten in het 
programma leggen nog eens uit 'wat waar! gebeurt. 


De inverse functie van de sinusfunctie (de boogsinus of aresinus) 
bestaat niet in MSX-BASIC. We lossen dit op door de inverse 
tangensfunctie (ATN in BASIC) te gebruiken. 


Het voordeel van een dergelijke computersimulatie, boven het uit- 
voeren van het natuurkundige experiment, is dat we heel gemakke- 
lijk verschillende brekingsindexen kunnen invoeren zonder steeds 
andere apparatuur op te hoeven stellen en dat de lichtstralen als 
dunne lijnen zichtbaar gemaakt kunnen worden. 
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100 REM programma 38 breking van 
Licht 
110 CLS 
120 INPUT "brekingsindex N';N 
130 PI=Z4&KATN(1) 
140 V=96:H=.S:RD=EPI/180 
150 CLS 
160 COLOR 15,4,4 : SCREEN 2 
170 DEF FNX(X)=INT(37+X/1.4+.5) 
180 LINE (FNX(O),V)=-CFNX(255) „V) 
190 LINE (FNX(O),190)=(FNX(255),190) 
200 LINE (FNX(255),190)=(FNX(255) ,0) 
210 LINE (FNX(255),0)=(FNX(O) ,0) 
220 LINE (FNX(O),0)=(FNX(O),190) 
230 REM stralen in medium 1 en 2 
tekenen;B=beta in graden; 
B1zbeta in radialen; Al= 
alfa ij rad.; SAssin(alfa) 
240 B=B+3;B1=B*RDEX1=0:Y1=2190 
250 X2=INT(9O*TAN(B1)+H) : Y2=V 
260 IF x2>255 
THEN AS=INPUTS$(1): 
COLOR 15,4,4 : STOP 
270 LINE (FNX(X1),Y1)=(FNX(X2) ,Y2) 
280 xi=X2:Y1=Y2 
290 REM alfa en sin(alfa) berekenen 
en totale reflectie nagaan 
300 S=N*SIN(B1) 
310 IF S>1 THEN GOSUB 1000 
ELSE GOSUB 2000 
320 GOTO 240 
330 END 
1000 REM subroutine totale reflectie 
1010 x2=X1+x1:Y2=2190 
1020 rif x2<255 
THEN LINE (FNX(X1),Y1)= 
CENX(X2) „Y2) 
ELSE X2=255:Y2=2INT(V+(255=-X1) 
/TAN(B1) +H) : 
LINE (FNX(X1),Y1)= 
CENX(X2) „Y2) 
1030 RETURN 
2000 REM subroutine breking 
2010 AT=ZATN(S/SAR(1-S*S)) 
2020 X2=INT(X1+9ORTANCAI)+H):Y2=0 
2030 If x2<255 
THEN LINE (FNX(X1),Y1)= 
CENX(X2) „Y2) 
ELSE X2=255:Y2=2INT(V-(255=-X1) 
/TANCAT) +H) : 
LINE (EFNX(X1),Y1)= 
CENX(X2) „y2) 
2040 RETURN 
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4, De speldenworp van Buffon (1773) 


Stelt u zich voor dat we in een plat vlak een aantal evenwijdige 
lijnen met een onderlinge afstand a tekenen. 


We werpen nu, zonder echt te 'mikken', een speld met een lengte c 
die kleiner dan of gelijk aan a is (c<a) op het vlak met de evenwij- 
dige lijnen. Hoe groot is nu de kans dat een speld een van de lijnen 
treft, dat wil zeggen snijdt of raakt? 


speld snijdt een lijn speld raakt een lijn 


Dit probleem kwam in 1773 op bij de Franse wiskundige Buffon na 
het zien van de Amerikaanse vlag met de 'stars en stripes'. Buffon 
heeft berekend dat deze kans gelijk is aan 


Omdat deze formule de constante n bevat, heeft men deze 'spelden- 
worp' vaak gebruikt om de waarde van n door simulatie te bepalen. 
We werpen hiertoe vaak, bijvoorbeeld een lucifer, op een papier 
waarop een aantal evenwijdige lijnen (met een onderlinge afstand 
die groter dan of gelijk aan de lengte van de lucifer is). Als de 
lucifer in n worpen k keer een lijn treft, dan geldt dat 


2e rn k 

an ongeveer gelijk is aan Rs 

Beide uitdrukkingen geven namelijk een indruk van de kans dat 
de lucifer een lijn treft. 

We kunnen als volgt de waarde van n benaderen: 


heks an 
an n ‘ak * 
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Op deze manier heeft de astronoom Rudolf Wolf in 1850 met 5000 luci- 
ferworpen voor n de waarde 3,1596 gevonden (n = 3,141592654....). 


Met een computer kunnen we gemakkelijk duizenden worpen simule- 
ren. Deze programma's vinden we in bijna elk informaticaleerboek. 
Dit zijn echter programma's die vrij lang draaien, maar waaraan 
niets te 'zien' is. Als we bijvoorbeeld in zo'n programma voor n de 
waarde 10.000 zouden intikken, zien we eerst een tijd niets en ver- 
volgens verschijnt de mededeling dat van de 10.000 keer 6.349 keer 
een lijn is geraakt, hetgeen voor n de waarde 3,1501024 oplevert. 


We gaan een programma maken dat ook dergelijke berekeningen uit- 
voert, maar dat bovendien elke speld die geworpen wordt laat zien. 
We zien het experiment als een film aan ons voorbijgaan. Bekijk 
hiertoe de onderstaande tekening. 


speld AB is 60 lang 


De manier waarop een speld valt kan met twee toevalsgetallen 
bepaald worden. De twee toevalsgetallen bepalen respectievelijk de 
afstand YM van de speld tot de daarboven liggende lijn en de hoek 
W die de speld met de lijnen maakt: 


Os YMsa ; YM is de afstand van het midden van de speld 
tot de daarboven liggende horizontale lijn 


OsWs180° ; Wis de hoek (in positieve zin, dat wil zeggen 
bij draaiing tegen de klok in) die de speld met 
de positieve x-richting maakt 


Hierna volgt een illustratie van hoe het eruit kan zien en het pro- 
gramma. 
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Ds PZ 


NP AN 5-7; 
ZE OPSREN 


REM programma 39 speldenworp van 


buf fon 
CLS 
PRINT "speldenworp van Buffon" 
PRINT "„smsengentaesnsmmnssen ne 


PRINT:PRINT 

INPUT "hoeveel worpen';N 

M=0O:H=.5:PI=4KATN(I) 

CLS 

COLOR 1,15,4 : SCREEN 2 

DEF FNX(X)=INT(37+X/1.4+.5) 

REM Lijnen tekenen 

FOR Y1=10 TO 190 STEP 60 
LINE (FNX(O),Y1)=-CFNX(256) ,Y1) 

NEXT Y1 

REM n maal gooien en tekenen 

FOR J =1 TO N 
XM=INT (226*RND(1)+30+H) 
YM=ZINT(6OERNDCI) +70+H) 
WEPIERND(1) 
DX=3OKCOSCWI:DY=3ORSINCW) 
X1=INTCXM=-DX+H) 
Y1=INTCYM+DY+H) 
X2=INT(XM+DX+H) 
Y2=INTC(YM=-DY+H) 
LINE CPNRENII ZY TI SCENE RI pT 2 
IF Y1>130 OR Y2<=70 

THEN M=M+1 

NEXT J 

AS=INPUT$S(1) 

COLOR 15,4,4 z: SCREEN O 

CLS 

PRINT "aantal worpen "iN 

PRINT "aantal doorsnijdingen';M 

PRINT "benadering voor pi"; 

INT(2*N/Mx*10000+.5)/10000 
END 


Kijk in de appendix hoe we de tekst en de waarden uit de regels 


400, 


410 en 420 onder de tekening op het graphic-scherm kunnen 


afdrukken. 
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Uit de waarden voor YM en W berekenen we de coördinaten van de 
uiteinden A en B van de speld. Het programma tekent hiermee de 
speld op het beeldscherm. 


Een speld treft een lijn als de y-coördinaat van A groter dan of gelijk 
aan 130 is of als de y-coördinaat van B kleiner dan of gelijk aan 70 
is. De x-coördinaten van A en B doen er in het geheel niet toe. 


Het programma tekent als 'speelveld' vier evenwijdige lijnen met een 
onderlinge afstand van 60. De spelden zijn trouwens ook 60 lang. 
Als alle spelden geworpen zijn kunnen door een toets in te drukken 
de waarden voor n, k en mn worden afgelezen. 


5. Prooi-roofdierpopulaties 


Het laatste educatieve programma is een (deterministische) simulatie 
van een ecologisch systeem. Het is een demonstratieprogramma voor 
een biologieles. 


Het ecologische systeem bevat gras, hazen en vossen. Tussen deze 
drie ecologische componenten gelden de volgende betrekkingen: 


1. De hazen eten gras en de vossen eten hazen. 

2. Als er meer gras groeit, neemt ook het aantal hazen toe. 
Deze hazen eten echter van het gras en verminderen zo hun 
eigen groei. 

3. Als er meer hazen komen, neemt ook het aantal vossen toe. 
Omdat vossen hazen eten, verminderen zij zelf hun groei, net 
zoals bij de hazen en het gras. 


Als we het aantal hazen op een bepaald tijdstip t aangeven met h(t) 
en het aantal vossen met v(t), kunnen we, volgens Lotka en 
Volterra (1920), voor het aantal hazen en vossen op tijdstip t+1 de 
volgende vergelijkingen opstellen : 


h(t+1) = h(t) +a.h(t) - b.h(t). v(t) vergelijking (1) 
v(t+1) = v(t) t+te.v(t).h(t) - dev(t) vergelijking (2) 


De toename van het aantal hazen tussen de tijdstippen t en t+1 is 
evenredig met het aantal hazen op tijdstip t‚ dus 


h(t*I)-h(t) =a.h(t) 


toename hazen Ea hazen op tijdstip t 


groeifactor 
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De afname van het aantal hazen is evenredig met het aantal aanwe- 
zige hazen (natuurlijk verloop) en met het aantal vossen, want die 
eten hazen, dus 


h(t+I)-h(t) =a.h(t) - b.h(t).v(t) 
groei van afname 
het aantal hazen ! aantal hazen 
toename 
aantal hazen 


Dit is vergelijking (1). We nemen in het programma aan dat er altijd 
genoeg gras voor de hazen voorhanden is. 


De toename van het aantal vossen is evenredig met het aantal aan- 
wezige vossen en met het aantal hazen (prooi). De afname van het 
aantal vossen is alleen evenredig met het aantal, omdat in dit 
systeem de vossen zelf geen prooidieren zijn. Voor de vossen krij- 
gen we dus: 


v(t+1)-v(t) = e.v(t).h(t) = d.v(t) 


groei van afname 
het aantal vossen fi aantal vossen 


toename 
aantal vossen 


De vergelijkingen zijn als zogeheten differentievergelijkingen opge- 
steld. Het gras speelt, zoals gezegd, eigenlijk geen rol; er is altijd 
genoeg om alle hazen te voeden. 


Bij het draaien van het simulatieprogramma hebben we de volgende 
waarden gekozen: 

X (= h(0)) = 200; Y (= v(0)) = 20; a=0,3; b=0,01; C=0,002 en 
d=0,5. 


De grafiek die het programma tekent laat heel mooi de groei en de 
afname van de hazen- en vossenpopulaties zien. Het aantal hazen 
groeit eerst, bereikt een maximum en neemt vervolgens af. De vos- 
sen groeien ook, maar later, bereiken later een maximum en nemen 
dan ook af, waarna de cyclus zich herhaalt. In de biologie noemen 
we dit een dynamisch evenwicht. 


U hoeft de waarden voor de coëfficiënten a, b, c en d maar iets te 
veranderen of het systeem kan ontregeld worden. Hierbij neemt òf 
het aantal hazen enorm toe òf alle hazen en vossen sterven snel uit. 
Deze eenvoudige simulatie toont aan hoe desastreus een kleine 
ingreep in een bestaand ecologisch systeem dat in een dynamisch 
evenwicht is kan zijn. 
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vossen 


REM progamma 40 roofdier-prooi- 
dier ecosysteem 
GES 
INPUT'pop. prooidieren(200) '';X 
INPUT'pop. roofdieren(20)'';Y 
INPUT'"groei prooidieren(.3)'"';A 
INPUT"afname prooidieren(.01)';B 
INPUT"groei roofdieren(.002)'';C 
INPUT"afname roofdieren(.5)'';D 
ELS 
COLOR 1,15,4 : SCREEN 2 
DEF FNX(X)=INT(37+X/1.4+.5) 
K=.3:H=.5:V=191 
FOR X1=0 TO 255 STEP 8 
XP=EX+(ARX=-BAXHY) 
YREY+H(CHXHY=-DHY) 
REM pop.prooidieren tekenen 
Y1ZINT(V=-K&X+H) 
X2=X1+8:Y2=INT(V=-K*XP+H) 
IF V2<0 Or v2>191 
THEN AS=INPUTS(1): STOP 
LINE CENXCKTD VII CENX CX2), VA) 
REM pop.roofdieren 
Y1=ZINT(V=K*Y+H) 
X2=X1+8:Y2=INT(V=K*YR+H) 
IF Nes OR- 422191 
THEN AS=INPUT$S(1): STOP 
ENE CÈNNCRTD pV TOT CENN END SLD 
X=XP:Y=YR 
NEXT X1 
AS=INPUTS$S(1) 
GOLOR 15544 # END 


hazen 
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Appendix 


In deze appendix laten we zien hoe een aantal programma's uitge- 
breid kunnen worden om leuke effecten te bereiken. Ook geven we 
nog een aantal tekeningen die u met de programma's kunt maken. 
Tot slot volgt dan nog een betere versie van het bolprogramma 25. 


Programma 4 


In de tekening, die programma 4 oplevert, zien we dat tussen een 
zeshoek en de daarbinnenliggende zeshoek zes driehoeken ontstaan. 
We gaan deze driehoeken om en om '‘inkleuren',. We doen dit niet 
alleen voor de driehoeken tussen de buitenste zeshoek en de eerste 
ingeschreven zeshoek, maar ook bij de driehoeken tussen de eerste 
ingeschreven zeshoek en de tweede ingeschreven zeshoek, voor de 
driehoeken tussen de tweede ingeschreven zeshoek en de derde 
ingeschreven zeshoek, enzovoorts. Er ontstaan dan drie 'zwarte'! 
spiraalvlakken. 


Als we een vlak, in dit geval een driehoek, dat geheel omsloten 
wordt door lijnen willen inkleuren, moeten we in dat vlak een punt 
lokaliseren. De schermcoördinaten van dat punt gebruiken we dan 
in de PAINT-opdracht om het vlak waarin dat punt ligt te kleuren 
(zwart, of een andere kleur). Omdat er een groot aantal driehoeken 
gekleurd moeten worden, moeten we ook een groot aantal coördinaten 
van punten in die driehoeken berekenen. Deze punten moeten op 
een systematische manier in de FOR-NEXT-lus (regels 230-350) 
bepaald worden. Als we in het X-Y-vlak een driehoek ABC tekenen, 
waarvan de coördinaten van de hoekpunten (Xa,Ya); (XB,Yg) en 
(Xc,YC) zijn (zie onderstaande tekening), dan weten we uit de 
meetkunde dat de coördinaten van het zwaartepunt (Z) van de drie- 
hoek gelijk zijn aan 


XA *XB * XC Ya *YB *Yc 
3 E 3 
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Het zwaartepunt van een driehoek is het snijpunt van de drie lijnen 
die elk een hoekpunt met het midden van de daartegenoverliggende 
zijde verbindt. In de tekening is Z het zwaartepunt en er geldt dus 
dat 


Xa + XB * XC MEt: Maak fo 
rh ne tene: sch 
Y C 


Xz XC 


= (K4tXBtXC0)/3 
= (YatYBtYC)/3 


We kiezen in ons programma steeds het zwaartepunt van de driehoek 
als inwendig punt. Om de coördinaten hiervan de kunnen berekenen 
moeten we dus de drie x- en de drie y-coördinaten van de hoekpun- 
ten weten. We hebben in de onderstaande tekening, die door pro- 
gramma 4 getekend wordt, één driehoekje als voorbeeld genomen. 
Van dit driehoekje zijn de hoekpunten A en B twee hoekpunten van 
de ingeschreven zeshoek, terwijl hoekpunt C een hoekpunt is van 
de daarvoor getekende (in ons voorbeeld de buitenste) zeshoek. 

We moeten dus steeds van twee zeshoeken de coördinaten van alle 
zes hoekpunten ter beschikking hebben. In het onderstaande pro- 
gramma hebben we, om de coördinaten van de hoekpunten van de 
zojuist getekende zeshoek te kunnen onthouden voordat de nieuwe 
berekend worden, twee arrays extra opgenomen, en wel XO(6) en 
YO(6) (van XOud en YOud). In regel 325 maken we XO(J) en YO(J) 
gelijk aan de coördinaten van de zojuist getekende zeshoek (X(J) 

en Y(J)), vlak voordat X(J) en Y(J) gelijk worden gemaakt aan de 
coördinaten van de hoekpunten van de volgende zeshoek (MX(J) 

en MY(J)) in regel 330). 
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X(K+1) 
X(K) 
XO(K+1) 


We zien vervolgens in de regels 263 en 264 hoe we de schermcoördi- 
naten van het zwaartepunt berekenen uit de twee hoekpunten van 
de nieuwe zeshoek (X(K), X(K+1) en Y(K) en Y(K+1) en één hoek- 
punt van de daarvoor getekende zeshoek (XO(K+1) en YO(K+1)). 
In de FOR-NEXT-lus van regel 262 nemen we de stapgrootte 2 om 
steeds één driehoek over te slaan. Natuurlijk hoeven we dit inkleu- 
ren pas te doen als de tweede zeshoek getekend is, vandaar de 

IF N=1 THEN 270 in regel 261. Het inkleuren gebeurt tenslotte in 
regel 265 met PAINT(FNX(X),Y),1. Denk erom dat ook hier FNX(X) 
in plaats van X genomen wordt anders wordt een verkeerd vlak 
gekleurd. 


Het resultaat staat op p.113. Daarna volgt het aangepaste program- 
ma 4. 
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REM programma & ingeschreven 
zeshoeken 

CLS 

COLOR 1,15,4 : SCREEN 2 

DEF FNX(X)=INT(37+X/1.4+,5) 

DIM X(6),Y(6) ,MX(6) ,MY(6) ,XO(6), 
Yo(6) 

PI=4&xATN(1) 

U=128:V=96:R=100:H=.5 

W=60#PI/180 

FOR J=0 TO 6 

Wiz J*w 

X(JDZINTCUIR*COS CWI +H) 

YOJDZINTCV=RESINCWI)D +H) 

NEXT J 
FOR N=1 TO 20 
FOR J=0 TO 5 
EME CENNENEI IN oe FAALII 
CENXOXCIHIDD YH) 

NEXT J 

IF N=s1 THEN 270 

FOR K=1 TO 5 STEP 2 

XEINTCOKEKDENCKHIDHKOCKHTD DD 
13) 

YZINTCCYEKDSYEKSID+YOCKHIDD) 
13) 

PAINTCENX (XD, YI, 1,1 

NEXT K 

FOR K=0 TO 5 

MXAK)DZINTC(OKEKDEXCKS IDD /2+H) 

MY(K)ZINTC(Y(KD EYCK) /2+H) 

NEXT K 

MXC6DEMX(O) 7: MYC6)=MY CO) 

FOR J=0 TO 6 
XOCIJDEKCIHIEYOCIHD ETI) 
XCJISMKCI) 2 YCI)EMY CI) 

NEXT J 

NEXT N 
AS=ZINPUTS(1) 
COLOR 15,4,4 : END 
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Programma 19 


We kunnen van de vliegekop, die programma 19 tekent, een 'echte' 
vliegekop met ogen maken. Hiertoe tekenen we twee ogen (cirkels) 
en kleuren het vlak dat door de cirkels en de daarbuitenliggende 
lijnen begrensd wordt. Het resultaat en de benodigde extra 
opdrachten in het programma staan hieronder. 


221 CIRCLE (FNX(148),96),10,1,0,2*PI, 
1.4 

222 CIRCLE (FNX(108),96),10,1,0,2*PI, 
ek 

223 PAINT (FNX(163),100),1,1 

224 PAINT CFNX(93),100),1,1 


Neem gerust een andere kleur door in de PAINT-opdracht de 1 te 
vervangen door een getal tussen 0 en 15. 
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Programma 21 


Allereerst geven we een paar voorbeelden van wat je zoal met pro- 
gramma 21 kunt tekenen. Je zou dit wellicht '‘computerkunst' mogen 


noemen. 


A=1/B=-100/C=-200 


A=3/B=100/C=200 
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Tekst op het grafische scherm 


Het zou leuk zijn op het grafische scherm de waarden voor A,‚ B en 
C af te drukken naast de tekening. Bevalt zo'n tekening, dan kun- 
nen we direct de waarden van A‚ B en C noteren. In de regels 231- 
238 van het onderstaande programma zien we wat we hiervoor moe- 
ten doen. In regel 231 openen we een bestand (met nummer 1) voor 
uitvoer naar het grafische scherm (GRP :). We kunnen niet van de 
PRINT-opdracht gebruik maken, maar we moeten met de 

DRAW "BM...'"-opdracht de grafische cursor naar een bepaalde 
schermpositie brengen (BM betekent Blank Move; alleen de grafische 
cursor verplaatst zich). In regel 232 gaat de cursor naar scherm- 
positie (10,100). Aangezien de tekeningen pas vanaf de x-coördinaat 
37 (zie FNX-functie in regel 170) getekend worden, hebben we links 
nog een veld van 0-36 bij 0-192 om tekst af te drukken. Rechts 
hebben we ook nog zo'n veld! 


Regel 233 zorgt ervoor dat de tekst 'A: ' gevolgd door de waarde 
voor A (hoogstens vier posities) afgedrukt wordt. De volgende 
regels doen hetzelfde voor B en C. In regel 238 sluiten we het 
tekstuitvoerbestand. 


Denk aan de '#1' achter de PRINT-opdrachten in de regels 233, 235 
en 237. Via deze omweg kunnen we dus toch op een grafisch scherm 
iets afdrukken. 
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REM programma 21 symmetrische 
figuren 
ELS 
PI=4&KATN(I) 
uU=128:V=92:H=. 5: RD=PI/180:K=100 
INPUT "EOEtS: AB en Cutn’zs ABC 
CLS 
COLOR 1,15,4 : SCREEN 2 
DEF FNX(X)=INT(37+X/1.4+.5) 
FOR W=0 TO 360 
T=WERDER=EKSIN(CHT) 
X2=INT(U+RKCOS (AT) +H) 
Y2=INT(V=RXSIN(B*T) +H) 
IF wW=0 
THEN X1=X2:Y1=Y2 
ELSE LINE(FNX(X1),Y1)= 
CENXCX2) „Y2): 
xi=X2:Y1=Y2 
NEXT W 
OPEN "GRP:'" AS #1 
DRAW "BM10,100"' 
PRINT#1, USING "A: HHHH";A 
DRAW "BM10,108"' 
PRINT#1, USING "B: #H###";B 
DRAW "BM10,116"' 
PRINT#1, USING "C: H##H#R!";C 
CLOSE#1 
AS=INPUTS$S(1) 
COLOR 15,4,4 : END 
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Programma 25 


De auteur Marcel Sutter heeft onlangs in het tijdschrift Mikro+Klein- 
computer een verbeterde versie van het 'bolprogramma' 25 gepubli- 
ceerd. Wij laten dit nu zien. Het programma tekent bollen in allerlei 
standen en met verborgen lijnen. De draaihoeken om de x-, y- en 
z-as alsmede de 'afstand' tussen de breedte- en lengtelijnen kunnen 
worden ingevoerd. Eerst volgt het programma, daarna enkele voor- 
beelden. 


100 


110 
120 


130 
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160 
170 
180 
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220 
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260 
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380 


REM bol met hidden lines en 
draaiing 

CLS 

PRINT “bol met Lengte en breedte" 

e Lgnen” 

PRINT:PRINT 

INPUT "draaihoek om x=as'"';A 

INPUT "draaihoek om y=as'";B 

INPUT "draaihoek om z=as"';C 

INPUT "afstand Lijnen(10=-45)"';D 

PI=4&XATN(1) 

U=128:V=96:R=I0OERD=EPI/180:H=.5 

S1=SINCA&RD) :S2=SIN(B*RD): 

S3=SIN(C*RD) 

C1=COS(A*RD) :C2=COS(B&RD) : 

C3=COS(C*RD) 


Ne 


REM rotatiematrix berekenen 
AX=C24CTzAY=-C2HSIEAZES2 
BX=CI#S3+ST4S2HCI 
BY=C1#C3-STS24S3:BZ=-STC2 
CXESTKSI-CTKSAHCF 
CY=STKCIHCIKS2#SI:E CZECH HCR 
REM hoge resolutie inschakelen 
CLS 
COLOR 15,15,4 © SCREEN 2 
DEF FNXCX)=INTC37/4+X/1.4+.5) 
REM bolomtrek tekenen 
CIRCLECFNX(U)D ,V),R,1,0,2*PI,1.4 
REM Lengtelijnen tekenen 
FOR L=0 TO 180-D STEP D 

F1=0 

FOR P=0 TO 360 STEP 5 

GOSUB 1000 'XX,YY,ZZ berekenen 
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390 IF YY>0 
THEN F2=0:F1=0 
ELSE XB=INT(U+XX+H) : YB=INT(V= 
ZZ+H):F2=1:IF F1=0 THEN 
X1=XB:Y1=YB:F1=1 ELSE 
LINECFNX(X1) YI) CFNX (XB) 
„YB):X1=XB:Y1=YB: FI=F2 
400 NEXT P 
410 NEXT L 
420 REM breedtelijnen tekenen 
430 FOR P==90+D TO 90-D STEP D 
440 F1=0 
450 FOR L=0 TO 360 STEP 5 
460 GOSUB 1000 'XX,YY,ZZ berekenen 
470 IF YY>0 
THEN F2=0:F1=0 
ELSE XB=INT(U+XX+H) : YB=INT(V= 
ZZ+H):F2=1:IF F1=0 THEN 
X1=XB:Y1=YB:F1=1 ELSE 
LINECFNX(X1) „Y1)=CFNX (XB) 
„YB):X1=XB:Y1=YBiFI=F2 
480 NEXT L 
490 NEXT P 
500 AS=INPUTS$S(1) 
510 COLOR 15,4,4 2 END 
1000 REM XX,YY,ZZ berekenen 
1010 XERKCOS(PARD)KCOS(LARD) 
1020 YaRKCOS(PARD)*SIN(LARD) 
1030 ZERKSIN(P*RD) 
1040 REM NVE MAAR CRA NVZ 
draaien 
1050 XX=AXX+AY*Y+AZ*Z 
1060 YY=BX*X+BY*Y+BZ*Z 
1070 ZZ=ECK*XACYKY+CZ*Z 
1080 RETURN 


119 


120 Appendix 


draaihoek om x-as: 0 
draaihoek om y-as: 40 
draaihoek om z-as: 60 
afstand lijnen : 30 


draaihoek om x-as: -50 
draaihoek om y-as: 20 
draaihoek om z-as: 60 
afstand lijnen Er T 


draaihoek om x-as: 90 
draaihoek om y-as: 30 
draaihoek om z-as: 40 
afstand lijnen : 30 
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Programma 36 


We gaan de kaart van Zwitserland wat verfraaien. Eerst geven we 
het gewijzigde programma. 


100 REM programma 36 gekleurde kaart 
van zwitserland met vlag 

TLD ELS 

180 COLOR 6,15,15 £ SCREEN 2 

130 DEF FNX(X)=INT(374+X/1.4+.5) 

140 K=1.35:H=.5:U=18: V=191 

150 READ X,Y 

160 X1=INTCU*K*X+H) 

170 Y1=ZINT(V-K*Y+H) 

180 REM begin herhaal Lus 


190 X2=INTCU4+K*X+H) 

200 Y2=INT(V=-K*Y+H) 

210 LINE (EFNX(X1),Y1)=(FNX(X2) ,Y2) 
220 xi=X2:Y1=Y2 

230 READ X,Y 

240 IF X<> O0 THEN 190 


250 REM einde herhaal Lus 

260 AS$S=INPUTS(1) 

261 IF A$="p"' THEN PAINT(FNX(128) ,96) 
„6,6:AS=INPUT$S(1) 

262 IF A$="'v!"' THEN U=128:V=96:K=1: 

COLOR15,11,11: GOTO 150 

263 IF A$='"'k"' THEN PAINT(FNX(128),96) 
„15,15:AS=INPUT$(1) 

270 COLOR 15,4,4 : END 

1 


290. DATA. rader ed ej we a eee dn veces ee 


520 DATA 0,0 

530 DATA =40,10,-10,10,-10,40,10,40 
540 DATA 10,10,40,10,40,-10,10,-10 
550 DATA 10,-40,-10,=40,=10,=10,=40 
560 DATA -10,-40,10,0,0 
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Wat zijn nu de veranderingen? 


regel 120: COLOR 6,15,15; afdrukkleur rood, achtergrond en 
kader wit 

regel 261: IF A$="p" THEN PAINT(FNX(128),96),6,6; 
als we nu een kleine p intoetsen, wordt Zwitser- 
land rood gekleurd 

regel 262: IF A$="v" THEN U=128:V=96:K=1: 
COLOR 15,11,11:GOTO 150; 
als we nu een kleine v intoetsen, wordt een wit 
omrand kruis (de vlag van Zwitserland) op het 
rode vlak getekend en het kader boven en onder 
de tekening wordt groen. Door GOTO 150 worden 
de extra data-regels 530-560 gelezen. Dit zijn de 
coördinaten van de hoekpunten van het kruis. 

regel 263: IF A$="k" THEN PAINT(FNX(128,96),15,15; 
als we nu een kleine k intoetsen, wordt het kruis 
wit gemaakt 


Het effect is nu als op de onderstaande tekening: wat daar zwart is, 
is op het scherm wit; wat wit is, is op het scherm rood en boven en 
onder de tekening zien we op het scherm nog twee stroken groen. 
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Programma 39 


Dit educatieve programma vraagt om het afdrukken van de benader- 
de waarde van pi op het grafische scherm, onder de tekening met 
de 'gevallen' spelden. We willen eigenlijk dit zien: 


Ea 
aantal raa k = A1 
benader ming vwoom pis an DD 


De regels 361-369 zorgen ervoor dat de tekst en de waarden uit de 
PRINT-opdrachten in de regels 400, 410 en 420 van het programma 
op p.106 nu op het grafische scherm worden afgedrukt. 
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REM programma 39 speldenworp van 
buf fon 
CLS 
PRINT "speldenworp van Buffon'' 
PRINT "=—vvaanvvaasaavvsvvz……. ze 
PRINT:PRINT 
INPUT "hoeveel worpen'';N 
M=0O:H=.5:PI=4kATN(I) 
CLS 
COLOR 1,15,4 : SCREEN 2 
DEF FNX(X)=INT(37+4+X/1.4+.5) 
REM Lijnen tekenen 
FOR Y1=10 TO 190 STEP 60 
LINE (FNX(O), YI) -(FNX(256) ,Y1) 
NEXT Y1 
REM n maal gooien en tekenen 
FOR J =1 TO N 
XM=INT(226*RND(1)+30+H) 
YM=INT(6OERND(1)+70+H) 
WEPI&RND (1) 
DX=30OKCOS(W)I :DY=ZORSINCW) 
X1ZINT(XM=DX+H) 
Y1=INTCYM+DY+H) 
X2=INTCXM+DX+H) 
Y2=INT(YM=-DY+H) 
GENE SENT NA IS GENK CED VB 
IF Y1>130 OR Y2<=70 
THEN M=M+1 
NEXT J 
OREN “‘orpe" AS #1 
DRAW "BM37,160"' 
PRINT#1,'"aantal worpen:'';N 
DRAW "BM37,168'' 
PRINT#1,"aantal raak :'';M 
DRAW "BM37,176" 
PRINT#1,'"benadering voor pi:''; 
PRINT#1, USING " #.HH#H#';24N/M 
CLOSER1 
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Het zal duidelijk zijn dat de in dit boek gegeven programma's onbe- 
perkt uitgebreid kunnen worden. Ze kunnen vast ook verbeterd 
worden; de invoerwaarden zouden bijvoorbeeld gecontroleerd kun- 
nen worden, zodat het programma niet ergens halverwege afbreekt 
door een foutieve schermcoördinaat. Wij hebben dat allemaal niet 
gedaan, omdat we de programma's klein wilden houden, waardoor 
alle aandacht op de tekentechniek gericht blijft. 
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Over de inhoud van dit boek 


Dit boek bevat 40 grafische program- 
ma's voor het programmeren met hoge 
resolutie in MSX BASIC. Dit betekent dat 
het scherm van uw MSX-computer 
verdeeld wordt in 256 puntjes horizontaal 
en 192 puntjes verticaal. Op dit scherm 
kunnen fraaie tekeningen gemaakt 
worden. 

Bij bijna elk programma wordt de 
nodige tekst en uitleg gegeven. Daarnaast 
bevat het boek veel illustraties van dat- 
gene wat de programma’s op het scherm 
tekenen. Alle programma's zijn getest op 
een MSX-computer en draaien foutloos. 
De programmatekst is direct vanuit het 
geheugen van de computer op de printer 
afgedrukt en zo in het boek opgenomen. 

De fraaie figuren berusten dikwijls op 
enige wiskundige vergelijkingen. In het 
boek wordt hier ook aandacht aan 
besteed. Dit vormt voor de niet-wiskun- 
dige lezer echter geenszins een belemme- 
ring het boek te gebruiken. Het alleen 


De illustratie op de voorkant is 
gemaakt met het onderstaande 
programma. 


100 REM programma 32 turtle-grafiek 
(LOGO-2) 

110 INPUT "coordinaten start ";X1,Yi 

120 INPUT "beginrichting(gr) ";W 


130 INPUT "verplaatsing "38 
140 INPUT "draaihoek linksom "4 DW 
150 INPUT "tosaname zijde "‚DS 


160 PI=4#ATN(1) 
170 He, SiRDePI/180rWisWeRD 
180 CLS 
190 COLOR 1,15,4 1 SCREEN 2 
200 DEF FNX(X)=INT(37+4X/1.4+.5) 
210 X2=INT (X1+8#COS CWI) +H) 
220 Y2=INT(Y1-S#SIN(W1) +H) 
230 IF X2<0 OR X2>256& OR 

Y2<o OR Y2>192 

THEN A$=INPUT$(1) 1 STOP 

240 LINE (FNX(X1),Y1)=(FNX(X2) ,Y2) 
250 Xi=X2iYi=Y2 
260 WeW+DW 
270 IF W>360 THEN WeW-360 
280 Wie=WeRD:5=5+DS 
290 GOTO 210 
300 AS$=INPUTS(1) 
310 COLOR 15,4,4: END 


maar intikken en draaien van de program. 
ma's geeft al zulke fraaie beelden dat het 
bestuderen van de erachterliggende 
theorie geenszins noodzakelijk is. 

Naast het tekenen van allerlei een- 
voudige en meer ingewikkelde functies 
(scholieren en leraren zullen nu hun oren 
spitsen) worden de beginselen van het 
driedimensionaal tekenen uitgelegd. 

Vijf programma’s laten zien wat voor 
figuren je met de taal LOGO kunt tekenen. 
Hieronder zijn ook de beroemde turtle- 
graphics. Ook bevat het boek vijf 
educatieve toepassingsprogramma's, 
waarin graphics de hoofdrol spelen. 

In een appendix wordt een aantal 
programma's uitgebreid om de lezer te 
laten zien dat de mogelijkheden met 
MSX-graphics haast onbegrensd zijn. 

Er zijn voorbeelden van het werken met 
kleuren (COLOR en PAINT) en van het 
afdrukken van tekst en getallen op het 
grafische scherm. 


